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On 6 Values in the Resistance of Metais 


By M. BLACKMAN 
Physics Department, Imperial College, London 


MS. received 13th April 1951 
ABSTRACT. The @ value to be expected on the Bloch theory of the temperature 


dependence of the resistance of metals is shown to be much larger than that observed 
experimentally. 


with temperature, MacDonald and Mendelssohn (1950) have determined the 

6 values, which we shall denote by @, for the alkali metals using the Griineisen 
(1933) function. For sodium, in particular, the 7° law was found to hold at 
low temperatures as predicted by Bloch (1929). The 6, value found is 202, 
while that derived from specific heat data, termed @), is 150. 

It is suggested that this 9) value should be increased to 202 on the basis of 
Simon’s (1926) theory of ‘thermal quantum jumps’. ‘The difference between 
the experimental 6, value for lithium, 270, and that deduced by Simon, namely 
512, is not considered. It should be pointed out that the elastic constants 
of sodium determined by Quimby and Siegel (1938) and by Bender (1937) 
lead to a 0 value (termed 6,,) of about 140. ‘This should equal the 0, value at 
low temperatures, and the agreement between 6)=150 and 6,,~140 definitely 
disposes of the 6)=202 value. In what follows we take 0, =06, as we shall be 
concerned with low temperatures. 

An important point which needs emphasis is that no agreement is to be expected 
between the 0) and 0, values. ‘The Bloch (1928) paper, and indeed all detailed 
accounts of Bloch’s theory, make it perfectly clear that the only lattice vibrations 
which play a part in the resistance are the longitudinal waves. ‘This is in marked 
contrast to what happens in specific heat theory, into which all vibrations enter. 
In fact, the 0, value is determined in the main by the transverse vibrations. 

It is surprising that, in experimental papers such as those of Griineisen and 
others, it should be tacitly assumed that 0, and @, should be identical. Even 
more surprising is that the same tacit assumption is made in theoretical pape:s, 
including, in some measure, Bloch’s original paper. 

On the Bloch theory we should expect 0g=0,, where 4, pertains to the 
longitudinal vibrations only. For an isotropic elastic medium this should be 


defined as i i 3 \18 
int ys Cy (=) , 


[ a recent experimental work on the variation of the resistance of metals 
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where A is the volume per atom, c;,, the ‘longitudinal’ velocity, # Planck’s constant, 
and k Boltzmann’s constant. For media with cubic symmetry pc,? =c,, where 
p is the density and c,, the ‘longitudinal’ elastic constant. This 6; value is 
appreciably larger than the 6) value. The latter, as mentioned above, can be 
calculated as 6, from the elastic constants, where 


Pee 5 oN 18 
oh vas ones (=) , 


Here Z is the mean wave velocity, which can be calculated from the elastic 
constants (see Born 1923). As examples we consider tungsten which is 
isotropic within experimental error, and aluminium which is nearly isotropic. 
In this latter case an average ‘longitudinal’ velocity is calculated, using the 
Born—-von Karman (1913) method. The results of the calculation are: 
Al: 0;,=707, 6)=415; W: 6,=636, 0p =384. 

The Bloch theory assumes inter alia elastic isotropy and an isotropic Fermi 
distribution in momentum space. The latter assumption should be reasonably 
good for sodium, lithium and copper. ‘The former does not hold at all well, 
though what is meant by ‘elastic anisotropy’ needs careful consideration. At any 
rate, we have to calculate, in these cases, an average c,, value. 

A reasonable method of calculating an average c;, value would be 


1 (S, .S,)? dQ 
= = | ee aya es (1) 
where S, is a unit vector in the direction of the displacement associated with a 
wave in a crystal, and S, a unit vector in the direction of propagation. The 
summation 7 is taken over the three types of waves, and dQ is an infinitesimal 
solid angle. ‘This form would be correct for the isotropic case, and also for 
the (purely hypothetical) special case of equal transverse and longitudinal 
velocities. Another method of averaging which is suggested by the Bee TS 
low temperature law, is obtained by averaging c;° instead of c;%. ‘These 
averages lead to 0, values denoted by 0,°, 0, respectively: 

Both methods of averaging were cated ehfouehs for the three metals, sodium, 
lithium and copper, by solving numerically the cubic equation for the corres- 
ponding elastic media. The elastic constants used for sodium were those 
measured by Bender (1937), those for lithium and copper were taken from the 
work of Fuchs (1936). For symmetry reasons only an eighth of the total solid 
angle in (1) needs to be considered. ‘This was divided into 90 equal sets. Each 
of these fixes a direction in space, and this, the three associated velocities. From 
these velocities the directions of the amplitude, associated with a wave travelling 
in the direction considered, are calculated. 


Substance 6,2) 6,5) nialih 
Na 282 260 202 
Li 540 530 260 
Cu 540 505 330 


The disagreement is obvious, and is independent of the method of averaging. 
The detailed calculation shows that the main contribution to 6; comes from only 
one of the three solutions of the elastic equations. This is mainly longitudinal 
in character and more isotropic (as the detailed calculations show) than the other 
two solutions which are mainly transverse in character. Further, the elastic 


it ia si i 
lot 
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anisotropy (defined as the ratio of the maximum and the minimum value of the 
velocity for a particular type of solution), is smaller for copper than for sodium 
_ or lithium, yet the agreement between 6; and 0p isno better. When one remembers 
that the resistance of metals at low temperatures is proportional to 6,~, the dis- 
crepancy shown above means a discrepancy of about five in the value of the 
resistance. 

It appears from the above, that while there seems to be qualitative agreement 
between the Bloch theory and the experimental data, there are considerable 
quantitative discrepancies. 

It is a rather curious empirical fact, established mainly by Griineisen (see, for 
instance, Griineisen 1933), that there is a remarkable correspondence between 
the 0, and 6, value for a large number of metals. It is true that the assumptions 
of the Bloch theory do not apply strictly to most of these cases; why the deviations 
from these assumptions should lead to 6g =0) is a matter which is in great need 
of clarification. It has been one of the purposes of this paper to emphasize that 
this need exists. 
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Lattice Stability and Oriented Overgrowth 
By M. SMOLLETT anp M. BLACKMAN 


Imperial College, London 
MS. received 29th Fanuary 1951 


ABSTRACT. The elastic constants of a two-dimensional ionic lattice, strained to fit on 
to the cube face of a substrate of the rocksalt type, are worked out. From these constants: 
the stability conditions of the strained layer are deduced. ‘These conditions are studied 
numerically in a number of selected cases where oriented overgrowth is observed in experi- 
ment. The calculations lead to conclusions which are in disagreement with a suggestion 
of Frank and van der Merwe as to the general mechanism of oriented overgrowth. 


Sl StNTRODUCTION 
N a series of papers Frank and van der Merwe (1949) have given a 
| mathematical discussion of a theory of oriented overgrowth. In this theory 
oriented overgrowth is associated with an initial layer uniformly strained to fit 
the substrate; where this initial process is not possible the orientation is random. 
The energy of a linear chain with quasi-elastic forces in a periodic potential 
field representing the substrate is compared in two cases: when the chain is strained 
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to fit the substrate, and when the chain contains dislocations. ‘The authors find 
that the strained chain has the lower energy when the degree of misfit (defined 
by the relative difference in spacing) is less than 9°4. Up to a larger critical 
misfit (about 14°% in the above example) the monolayer can be deposited metastably 
in exact fit on the substrate at sufficiently low temperatures. Since according to 
Frank and van der Merwe the dislocated lattice is mobile on the surface, complete 
oriented overgrowth of one crystal on another is expected only if the first monolayer 
can be formed on the surface under ‘ sub-critical’ conditions. 

The above theory is claimed to be ‘in general agreement’ with experiment. 
The critical misfits of 99% or 14°% are not strictly applicable to experimental 
cases, and it is assumed that a larger critical misfit can occur, for instance, in 
oxides. A particular case which is specifically quoted is that of zinc oxide (Cabrera 
and Mott 1949), where the misfit is 20%. 

Such misfits involve remarkably high strains in the initial layer formed on 
the substrate. As it is known that strains of this order would be close to the 
theoretical elastic limit to which a crystal can be strained (see Zwicky 1923), it 
was decided to study the stability of a strained two-dimensional layer. In many 
cases the overgrowth consists of an alkali halide of the NaCl type, on a substrate 
of similar type. A two-dimensional ionic lattice, in arrangement similar to the 
(100) face of rocksalt, fitting on the (100) face of a selected alkali halide of the 
NaCl type, has therefore been chosen as a suitable model. ‘To study the stability, 
the elastic constants of the two-dimensional layer are worked out in selected 
cases (Tables 2(a), 2(b)) where the layer is strained to fit the substrate. The 
effect of the substrate is taken into account. Jn all these cases oriented overgrowth 
is observed in practice. It will be shown that the strain leads to instability in all 
these cases, and hence the Frank—van der Merwe theory cannot account for these 
experimental results. 


§2. THE STABILITY OF A TWO-DIMENSIONAL LATTICE 

The equilibrium of two-dimensional lattices is considered in this work from 
the point of view of elastic stability with respect to small vibrations. The stability 
conditions which we shall use follow from the equations of motion of a two- 
dimensional ‘regular’ elastic medium, which can be written in the form 

[pox? — 77(Cy18q7 + C4aSy7)] Ug — (Ce + Cqa)SasyU, =0 | 

= 77(C19 + C44) 5,U 2 + poo? = 77(C745,7-F C4485) U,=0. i 
U,, U, are here the components of the amplitude U of a wave in the medium, 
the complete wave form being given by Uexp {—ir(xs,+ys,)} exp (—iw?). 
The density of the medium is p, and its elastic constants are given by ¢1,, C19, C44y 
S,, S, are the components of a unit vector normal to the wave front, r=2z/A 
where A is the wavelength, and w is the angular frequency. 

The conditions of stability may be deduced from (1), and these are (Blackman 
1938) that ¢,>0, ¢1;—¢y.>0, ¢y4>0 for stability. If any of these expressions 
vanish or are negative, w becomes imaginary and the medium is unstable. The 
stability problem reduces to the calculation of ¢,,, cy. and c4, for an ionic lattice 
of the type described in §1, and then the determination of the vanishing points 
Of Cy;—Cyg) £41 OF Cyy. Alternatively we can calculate for a specific case of a 
strained lattice, whether these quantities are positive, zero or negative, and hence 
deduce the stability. 
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$3. THE ELASTIC CONSTANTS 

A method due to Born (1923) is used for the evaluation of the elastic constants. 
This proceeds by formulating the vibrational equation for the lattice and then 
specializing the equation for long wavelengths (and low frequencies). The 
limiting equations can then be identified with the equations of motion of the 
elastic continuum, and the elastic constants can then be calculated in terms of 
lattice sums. This procedure can be carried out for a uniformly strained lattice 
in the same way as for a lattice at equilibrium, and the equations of motion are 
formally the same. The only difference is that the Cauchy relation cy, =c4, no 
longer holds in the case we are interested in, but we still obtain the equations (1) 
as the limiting equations. 

The lattices we are dealing with contain, strictly speaking, two particles per 
cell. When the masses of the two particles are equal, we can work with one 
particle per cell if the repulsive forces between any pair of particles are taken to 
be independent of the type of particle. The equation of motion of the positively 
charged particle is then identical with that of the negatively charged particle. 
Further, as we are interested only in the long-wave limit of the equations of motion, 
and as these depend only on the mean mass of the particles (or the density), we 
can work with one particle per cell even when the masses are different, provided 
we use the mean mass M. The equations of motion will be correct in the limit 
we need, though not outside this limit. 

We can now write the vibrational equations in a form similar to that given © 
by Born (1923): 

(w®M+A,,)u,+ Ag uy =0, | 
ot aes (2) 
ee () 


As +(w*?®M+ A,,)u 


a 


where w,, uw, are the components of the amplitude of a particle of mass M in the 
reference cell, and 
57 ASSET It aaa o> 40 Fearn 6 CuBEWEST SMa MN I te UMM RDAAY cote (3) 
MM j 

Here the summation m,, m, extends over all cells, fj’ is the appropriate 
derivate of the potential energy ¢ of a pair of particles whose distance is |r, »,| - 
For the particular type of square lattice with which we are dealing 7%, ,,. = x7, +, 
with x,,, = 119) Vm, = M7 Where ry is the lattice spacing and the m values are, 
of course, integral. ‘The unit vector s specifies the direction of propagation of 
a plane wave in the lattice, w is the angular frequency and r = 277/A, where Q is the 
wavelength. 

The coefficients A,, are expanded for small values of 7. In equation (2) 
only the second order terms in s,, s, will be retained. The coefficients of these 
terms can then be directly related to the elastic constants in (1). 


$4.-LATTICE SUMS 


The potential energy of a pair of particles in the lattice is written in the form 
d(7r)= +e7/r+ f/r*’. The sign in the first term depends on whether the two 
particles have like or unlike charges. It is usual to split the lattice sums (3) 
into two parts, a ‘Coulomb’ part and one depending on the repulsive forces. 
These latter sums converge quickly and can be summed directly, while the 
Coulomb terms require a special technique. 
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A typical sum in (2) involving the Coulomb part of the potential has the 
form 
Ale = E (PS /i2™ exp [= i7(%nS2+IYm,Sy)| 


M1 Ms 


ee BerNx? f 
= (- = ot =") XD [27 (NaS yee) eee (4) 
Mm, Me Vin, Me Ving Ma r 
Replacing x,,. by m7, Vm, bY moro, and putting 78,%9= —27hy, 75,79 = —27the, 
we can write the above in the form 
B 2m," —m," : : 
A®=— LY —+— exp [2mi{m,(h, +4) + m(hgt+4)}], oe eee (4a) 


ro My M2 (m,” nF mM, )°!? 
where the factor +1 is replaced by exp {27i(4m,+4m,)} and is incorporated 
into the exponential. Similarly 
Se mm 

Ae) — 73 a a2 m2 exp [27t{m,(h, +4) + mo(ho+4)}]. «20... (46) 
Sums of the type given in (4a) and (4) can be evaluated numerically with the 
help of the Epstein zeta functions (Epstein 1903, 1906) which are defined as 
Oral) _ og, exp {271(m,h, + mohz)} 
 eptind ie (s) = baees (m4? + m2)” =) a ete (5) 
where the prime indicates that the (0,0) term is excluded. It can then be shown 
that the zeta function we need can be expressed in the following form: 


eS 0 0 
id r(3)2| My+y hgt+% ) 


=— : + | dzz'2 XY” exp {—s[m,? + m,?] + 2rim,(h, + 4) + 207m,(h, + 4)} 
1 


M1 Ms 


+{ dzz-? & exp {—al(m, +h, +4)? +(myt+hot+h)}} ...... (5 a) 
1 


My Ms 


Z 


Here again the prime indicates that the (0,0) term is to be omitted; I denotes 
the usual gamma function; s is an integer, and the summation extends from 
—o to +0. Inthe above form the summations are rapidly convergent. 

The zeta function given above is not in the final form required for the sums 
occurring in (4a), (46), which contain factors of the form m,?, m,?,m,m,. These 
sums can be obtained by suitable partial differentiations with respect to h, or hg. 
This is illustrated in (6): 

0” 0 0 


ig _ gg yr My exp [2rr0{my (hy + 5) + mo(hy + 4)}] 
Apa ee eet (s) = —4n? XY =e 


My Mz (m,? + m,?)°? 


If we substitute s=5 in equation (5a) and carry out the differentiation as in 
(6) we find 


coe Nae oO 
—5/2 =| eee ae 
a (5) a3 | aes 


=| dz 9°S! —4%m,? exp [ —72(am,2 + my?) + 2rrimny hy, +4) + 2rrimmg(hy + 3)] 
=a | 


+ | dz x V?L477(m, + hy +4)? exp [ —az{(m, + hy +B)? + (my + hy +4)?}] 
J4 
- | dz 2-9? — 27 exp[—a2{(m, +h, +4)? +(myth,t+h)}].  — ... eee (7) 
1 
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To obtain the terms needed in the approximation we need an expansion of this 
result up to the second order in h, and h,. ‘The terms which will interest us will 
be the second order terms. ‘These are 


h,?[X’87r4m,4 exp [7i(m, + m2) | F(3; m,; me) 

+ X824(m, + 4)4F(3; m, +4; m.+4) 
+E = 2M + YEFGS m +45 ma +4) 
+ L672 F(—4; m, +4; m,+}4)] 

+ he2[S! 8ar4m,2m,2 exp [7i(m, + m,)]F(3; m,; mz) 
+ X824(m, + 4)?(m, +4)? F (33 m, +4; m.+4) 
+3 —4n9l (m+ 4) + (my + SPECS my +45 m+ 3) 
+D2n2F(—4; m,+4; m,+4)], Sern ate © Re (7a) 


where F(n; XE »)=| 2” exp {—7a(x? + y")} dz. 
i 


We have indicated how to reduce a particular lattice sum to the required form. 
‘The integrals which occur in (7a) can all be transformed to error functions, and 
the coefficients of h,?, 4,2 are evaluated numerically without difficulty. In this 
way all the lattice sums can be dealt with. 


§5.-THE ELASTIC CONSTANTS 
Collecting the electrostatic and repulsive terms together and comparing 
equation (2) with (1), the elastic constants are found to be 


Pies ea ioes aailign.g ih pli oobi pet acl 
1p 6s 8 li? iOS &>3 i> 
- "0 "0 —. ..(8) 
a feresn its ato) | 
ea = 104, Hy 7-19 


‘These are valid for a stretched or compressed ionic lattice of the type described 
in §1, where 7 is the actual lattice spacing of the strained lattice. 

The value of 8 has been calculated, where required, from the properties of 
the (three-dimensional) crystal, using the same assumption as above for the 
potential energy of a pair of particles. The equilibrium condition for the 
potential energy ® of a crystal of the rocksalt type can then be expressed as 

a -0-n( = an) fomen=d, 

or 
where WN is the number of cells in the crystal, A = 1-7476, B=6-6288. With the 
known equilibrium distance d, we find the value of B. 

We next consider the effect of the substrate on the values of the elastic 
constants. In the calculation a positive ion in the strained lattice is assumed 
to be directly above a negative ion in the substrate, and a negative ion in the lattice 
to be above a positive ion in the substrate. 

The method used for calculating the contribution of the substrate is similar 
to that given for the lattice. In this case, however, the summations in (2) are 
carried out between the particle in a reference cell in the lattice and those in all 
layers of the substrate. As far as the electrostatic part is concerned, the summation 
is carried out plane by plane, this summation being rapidly convergent. A brief 
account of the method is given in Appendix J. The distance of the lattice 
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from the substrate is found by working out the mutual potential energy and 
making this a minimum (Appendix II). 

It should be noted that the substrate is treated as a fixed potential field. 
This may be justified by the fact that the contributions of the substrate are 
relatively small, as may be seen from Table 1. Only the Coulomb term 
contributions are considered, as the contributions to the repulsive terms are 
negligible in comparison. 

In Table 1 Ac&) and A(c,,—¢2)" are the contributions of the substrate to 
the value of c,, and of (¢,;— C2), these being electrostatic in nature, as explained 
above; cand (¢,;—Cyg)" are the electrostatic parts of the elastic constants 
given in (8). It should be noted that Ac{? is always negative, so that it cannot 
assist the stability of the deposited layer. The contributions are given for 
various values of a where av, is the distance of the strained lattice from the substrate, 
7) being the lattice spacing in the strained state. ‘The values of a were chosen 
to suit experimental cases. It can be seen from the Table that the influence 
of the substrate diminishes rapidly as a becomes large. 


Table 1 
a 0-92 1-00 1-16 1:30 ° 1-60 2-00 
Act)/e(e) 0:22, 018, “S012 0-08, 20-06, 0.00: 
Aer — C42) /(C44 — C10) +0°13,.. --0+10,- -10-055 | .4:0:03, =-0+02,) 120-008 


§6. LATTICE STABILITY 

The effective elastic constants of the deposited layer (i.e. the strained lattice) 
are found by adding the contribution of the base to the ordinary elastic constants. 
given in (8). ‘These new constants are used in estimating the stability. 

Calculations have been carried out for a number of cases where, in experiment, 
an alkali halide of the NaCl type is deposited on the cube face of an alkali halide 
substrate of this type. Whenever (¢4,+ Acq.) or (€y, — C9) + A(C,1 — yg) becomes 
negative, the deposit is unstable in accordance with the criteria given in §2. 
It has been found that the c,, value vanishes first for a compression and the 
(Cy; — yg) value for an extension of the lattice, so that only those two expressions 
have been given in the tables. As will be seen below, in all the cases given in 
Tables 2(a) and 2(4), an oriented deposit is found in practice. The data on 
orientation are taken from van der Merwe (1949), except for the first case in 
Table 2 (a), which was reported by Schulz (1950). 


§7. DISCUSSION 

We have evaluated above the elastic constants of a two-dimensional alkali 
halide layer strained to fit on to a substrate consisting of an alkali halide crystal. 
In each of the cases quoted in Tables 2(a) and 2(4) oriented growth is observed 
experimentally with the (100) face parallel to that of the cubic substrate. 

Our theoretical calculations are designed to examine the dynamics of the 
uniformly strained lattice, treated statically by Frank and van der Merwe. In 
each of the examples investigated in ‘Tables 2(a) and 2(b) instability is found in 
the initial deposited layer if this is strained to fit the substrate. We must 
therefore conclude that the orientation observed in experiment cannot arise from 
an initially uniformly strained layer. It seems reasonable therefore to suppose 
that orientation is associated with an initially ‘distorted’ lattice where by 
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‘distortion’ we mean any type of displacement of a lattice particle from its original 
position leaving the mean lattice spacing more or less constant. 

Table 2 (c) has been included to show that our calculations will give dynamical 
stability in cases where orientation is found, and also in cases where orientation 
is not found. This should emphasize the point that stability of the type we 
have been discussing is a necessary but not a sufficient criterion. It is quite 
possible for a particular arrangement to be dynamically stable and yet to be far 
from the state of minimum free energy. 


Table 2 (a) 

S D C a Cas Now Total Stability io) 
LiF KBr 64 2:90 —12:21 0:00 —12-21 Unstable ae 
NaCl RbI 30 1:30 — 3-06 —0-42 — 3:48 Unstable ae 
NaCl KI 25 1-24 — 1:05 —0-45 — 1:50 Unstable ~ 

Table 2 (6) 

Ss) D E a (C11— C12) A(¢y;— 12) ‘Total Stability © 
KC} LiBr Ais 0:96 —2:81 +0-65 —2:16 Unstable -+ 
KCl AgClal2 0-92 —2-65 +0-71 —1-94 Unstable ae 
KCl NaCl 10 0-93 —2:25 +0-68 —1:57 Unstable =i 

Table 2 (c) 

S D Cc a Chg Gan Total Stability o 
NaCl KBr 17 1:16 1:35 —0:59 0-76 Stable _ 
KCl RbI 17 1:16 1:14 —0:43 0-71 Stable ct 
NaCl Nal 15 1-17, 1:86 —0-56 1-30 Stable se 
NaCl RbCl 17 1-18 1:54 —0:59 0:95 Stable - 


S, substrate crystal; D, deposited growth; C, percentage misfit for a compression; 
E, percentage misfit for an extension. The elastic constants are given in units of 
1072? (dyne/cm.). ©, orientation found in experiment; +, orientation present; —, no 
orientation. 


It should further be pointed out that though Frank and van der Merwe 
(1949) claim agreement between theory and experiment on the basis of van der 
Merwe’s review (1949), other reviews of the evidence on oriented growth 
(Thomson and Cochrane 1938, Thomson 1949) do not show any marked 
preference for the view that a critical degree of misfit is the deciding factor 
between complete and random orientation. 

Nor is it clear how the idea of critical misfit explains the temperature depen- 
dence of the orientation of evaporated layers. For instance, in the case of silver 
evaporated on to rocksalt at low temperatures a random orientation is observed, 
while at about 150° c. silver will give a (100) orientation on the cube face of rocksalt 
(Bruck 1936). Further, in some cases orientation changes with temperature 
(see Thomson and Cochrane 1938); the orientation is definite above and below 
certain temperatures, but it passes through a partially oriented phase in the 
intermediate region. This again is difficult to reconcile with a critical misfit 
criterion. 

The Frank—van der Merwe hypothesis has also been applied (Cabrera and 
Mott 1949) to chemical growth, using the idea that in oriented chemical growth 
the first layer (or layers) has the same spacing as the substrate. This concept is 
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applied particularly to oxide films, but presumably the theory is intended to be 
general. The main experimental evidence adduced is the pseudomorphic 
zinc oxide structure found on zinc by Finch and Quarrell (1933, 1934). However, 
a recent investigation (Raether 1950) has been made of the structure of the thin 
protective zinc oxide layer on the cleavage face of a single crystal of zinc. The 
oxide is found to have the normal structure, though this should bea very favourable 
case in which to observe the pseudomorphic form. The results of Raether 
are in agreement with independent work carried out at the Imperial College by 
Mire DON D ylnicas: 

A large number of investigations has been made on the growth of compounds 
on single crystals using the electron diffraction technique. The general rule is 
that the overgrowth has the normal crystal structure. This point could be 
studied in more detail if crystals with flat surfaces (produced by cleaving, 
electrolytic polishing or thermal etching) were used. With these surfaces 
layers of from five to ten atomic spacings should produce an observable pattern. 
It should then be possible to decide whether structures differing from the normal 
structure play a part in the mechanism of chemical growth or in the general 
process of oriented overgrowth. 
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APPENDIX I 

When calculating the effect of the substrate, the additional contributions 
appear in the form given below. ‘The expression for Ac,, is given as an example: 

Ac, = koh D> m,” exp {7i(mj +ms)} — 3e x, my* exp {at(m, a my) } 

ote M, Ms O (m,? a Ms” a0 b?) 32 279° My Ms D (m? a my” =f 6?) 
Here m, and m, refer to lattice points in the planes of the substrate; b7, is the 
distance of a substrate plane from the particle in the cell which is taken as origin. 
The summation takes place over all m,,m, and 6. To evaluate the above sums 
which contain 6 we use a modified function which is defined by 
00 Lg ye yy explemitmy ty +4) + (y+ I} 

hy sf ss hy 1% 3 My Mo (m2 a ms? 45 eg 
and this zeta function can be expressed in the form 


ieee i 0 0 
2 r(5 “lh +4 hy +4 (e 


Z 


= | -dsa* 2 exp (—nb22) E exp[—m2(m,2 + m,2) + 2mi{m,(hy +4) + mo(hy + 1)}] 


il M, Ms 
oO 
+] dzz*? & exp(—7b?/z) exp[—a{(m, +h, +3) + (mz +h, +43)}].- 
1 MMs 
This form enables us to compute the sums in the expression for Ac,,. The 


second integral cannot be reduced to an error function and must be calculated 
numerically. 
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ATORIEON DX ail 


To evaluate the distance of the stretched lattice from the substrate we form 
the potential energy of a representative particle in the lattice in the potential field 
of the substrate. ‘The minimum value of this energy is found as a function of a. 
‘The form of the potential energy used for the Coulomb terms is (see Born 1923, 
p. 738) 


see yy oXP { — 1b(m,? + m,*)1?} 
YQ mms b (m,? + m,”)"? 


where the summation is over all positive odd values of m,, m, and over all values 
of 6 which varies from plane to plane. The potential energy for the repulsive 
forces is 

us XX (m,? + mo” + 6?) 9? 

To mms b 
where the summation is over all integral values of m, and m, from — oo to + ©. 
The quantity 5 takes the values a,a+1,a+2etc. The value of f in the preceding 
equation is taken as the mean of the values for the substrate and the deposit. 


REFERENCES 

-BLackMaN, M., 1938, Proc. Roy. Soc. A, 164, 62. 

Born, M., 1923, Atomtheorie des festen Zustandes (Berlin : Teubner). 

Bruck, L., 1936, Ann. Phys., Lpz., 26, 233. 

‘CaBRERA, N., and Mort, N. F., 1949, Rep. Prog. Phys., 12, 163 (London: Physical Society). 

EPSTEIN, P., 1903, Math. Ann., 56, 615; 1906, Ibid., 63, 205. 

Fincn, G. I., and QuarRrELL, A. G., 1933, Proc. Rey. Soc. A, 141, 398; 1934, Proc. Phys. 
Soc., 46, 148. 

FRANK, F. C., and vAN DER MErRwE, J. H., 1949, Proc. Roy. Soc. A, 198, 205, 200, 125. 

VAN DER Merwe, J. H., 1949, Disc. Faraday Soc., 5, 201. 

RAETHER, H., 1950, 7. Phys. Radium, 2, 11. 

Scuuz, L. G., 1950, Phys. Rev., 78, 638. 

‘THomsoNn, G. P., 1948, Proc. Phys. Soc. A, 61, 403. 

“THomson, G. P., and CocHRANE, W., 1938, Theory and Practice of Electron Diffraction 
(London : Macmillan). 

TIWICKY, E.,, 1923, Phys:.Z., 24, 131. 


Some New Bitter Patterns on a Single Crystal of Nickel 
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ABSTRACT. Some new Bitter figures have been obtained on a single crystal of nickel 
whose surface is a (111) plane when magnetized along a [110] direction, and a curve obtained 
of domain boundary spacing as a function of the effective magnetic field. Pctterns of a 
complex nature, probably due to closure domains, were also obtained on a (112) plane. 


$i INTRODUCTION 
HE present work was undertaken to extend our knowledge of the domain 
| magnetization processes undergone by a single crystal of nickel when a 
magnetic field of varying intensity was applied to it. Much work has been 
done using the technique originated by Bitter (1932); most of it, however, has been 
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devoted to the magnetic properties and domain structures of crystals of silicon— 
iron, notably by Bozorth, Williams and Shockley (1949), and by Bates and Neale 
(1950). Nickel is of particular interest and difficulty because it possesses eight 
crystallographic directions—the [111] or trigonal axes, which are easy directions 
of magnetization, while iron and silicon—iron possess only six—the [100] or 
tetragonal axes. 

Kaya and Sekiya (1935) published a paper containing some beautiful photo- 
graphs showing powder patterns they had obtained on (100), (110) and (111) 
planes in nickel. Their method of preparation of the surface was, however, 
likely to result in strain, and as was found in the present work, this affects the 
nature of the powder patterns obtained very considerably. Indeed, on account 
of the large magnetostriction coefficient of nickel, the greatest care must be taken 
throughout the preparation and manipulation of the specimen to avoid distortion. 


§2. PREPARATION OF THE SPECIMEN 


The crystal was cut into the form of a slab having approximate dimensions 
15mm.x5mm.x3mm. Its main surface was approximately parallel to a 
(111) plane, and its sides roughly parallel to a [110] direction forming (112) 
planes. 

It was mounted in a specially designed holder and the orientations of the 
required crystallographic planes and directions relative to the surfaces cut were 
found with an x-ray goniometer. The crystal surfaces were then carefully 
ground on a high-speed emery wheel so that they corresponded, within three 
degrees, to the required atomic planes. During this part of the preparation 
care was taken that the specimen never became unduly warm, in order to keep: 
distortion and disturbance of the czystal to a minimum. 

The surfaces were then polished mechanically on emery papers down to 0000. 
This removed the layer of disturbed metal produced during the cutting and 
grinding processes but left in its turn a thin polycrystalline Beilby layer which was 
covered with very fine scratches. ‘To remove both and expose the true crystal 
surface, the specimen was polished electrolytically (Bates and Mee 1950) in a 
solution consisting of 100cm*® of orthophosphoric acid, 25cm? of sulphuric acid. 
and 75cm? of distilled water. This electrolyte was best used at a temperature 
of 40-50°c. in an inverted-bottle type of cell. The specimen formed the anode 
and was surrounded by acylindrical iron cathode. Acurrent density of 1 amp/cm? 
of anode surface was employed, and a polishing time of from three to four minutes. 
was required to remove the scratches. 

To remove strains the specimen was annealed zm vacuo at 900° c. for about 
two hours and then cooled slowly at the rate of 120°c/hour. The importance 
of doing this is shown by the unusual ‘maze’ pattern, obtained on the surface 
before annealing, shown in Figure 1 (a) (Plate *). 

The crystal specimen was placed in a horizontal field between the pole pieces 
of a small electromagnet of special design, giving but little stray field in a vertical 
direction. ‘The pole-pieces were of similar shape to those used by Bozorth,, 
Williams and Shockley (1949) in their work on silicon-iron. The figures were 
made manifest with a colloidal suspension of magnetite in soap solution, 


developed by Elmore (1938). 


* For Plates see end of issue. 
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§3. EXPERIMENTAL PROCEDURE 

To remove the effect of stray magnetic fields encountered during the pre- 
paration of the specimen it was first saturated by passing a heavy current through 
the electromagnet coils and then carefully demagnetized by reversals. On 
placing colloid on the crystal, no pattern was observed while it was in a demag- 
netized state, but on applying magnetic fields of gradually increasing intensity 
powder patterns were obtained. Preliminary experiments on a single crystal 
disc whose equatorial plane was a (111) surface showed that patterns were only 
obtained when a field was applied along a [110] direction (face diagonal). 

The patterns obtained were straight lines perpendicular to the direction of 
the applied field. Examples of these patterns are given in Figures 1 (d), (c) and (d). 
The spacing of these lines was sufficiently regular to enable some quantitative 
measurements of the line spacing to be made. The spacing varied in a regular 
manner as the applied field was altered. Before each measurement of a spacing 
the crystal was carefully demagnetized by reversals and then put into the cyclic 
state corresponding to the required field value. 

Graph (a) of Figure 2 shows the variation of line spacing with strength of 
the effective magnetic field, measured by a magnetic potentiometer. It will 
be seen that it bears a marked resemblance to the Néel (1944) case for a single 
crystal of iron, in the form of a slab bounded by (100) and (110) planes with a 
magnetic field applied along a [110] direction. No useful patterns of any 
regularity were observed when a magnetic field was applied in a vertical direction. 
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Figure 2. (a) Graph of line spacing against applied magnetic field. (b) J, H curve for nickel crystal 
strip, (111) plane, [110] direction of magnetization. 


Colloid was also applied to a polished (112) plane on the side of the above 
crystal and patterns obtained. ‘These were rather more complex in nature 
(Figure 1(e)). As can be seen from some of the accompanying photographs, at 
relatively low values of the applied magnetic field the lines appear very broad, 
with side branches parallel to the field. As the field is increased the lines become 
more definite and also move closer together, following the lines on the (111) 
plane. A peculiar branching of the powder lines can be seen at a crystal 
boundary on one of the photographs, Figure 1(f), and it is thought that this may 
be due to surface closure domains. If this is correct, closure domains will be 
formed on both (111) and (112) planes; that these closure domains are likely 
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from the energetic point of view can be shown by energy considerations given 
below. 

The magnetization curve for the crystal was obtained experimentally. The 
effective field acting on the crystal was measured directly by a very small magnetic 
potentiometer in series with a sensitive Tinsley galvanometer, a standard mutual 
inductance being used for calibration purposes. ‘The intensity of magnetization 
was found by using a search coil wound directly on the specimen itself. As can 
be seen by comparing graph (b) of Figure 2 obtained in this way with that of 
Honda et al. (1935), the measured intensity of magnetization is somewhat lower 
than that generally accepted. This may perhaps be explained by assuming that 
the most stable domain configuration corresponding to any particular value of 
the magnetic field takes time to form. 


§4. DISCUSSION OF RESULTS 

Nickel has eight directions of easy magnetization—the [111] body diagonals. 
of the face-centred cubic lattice. ‘The portion of the magnetization curve over 
which the patterns appear suggests that it is a two-phase region’(Néel 1944). ‘This. 
will extend from the main knee of the magnetization curve towards saturation.. 
Now the (111) plane has no easy direction of magnetization in it: the nearest 
makes an angle of approximately 20° with the surface and an angle of 35° with 
the applied magnetic field. ‘Thus if the crystal is divided into domains, as the 
powder patterns suggest, these will be in the form of sheets perpendicular to the 
direction of the applied field, i.e. the [110] direction. If these were magnetized. 
along the two easy [111] directions nearest to the direction of the applied field, 
they would, if they extended up to the crystal surface, have associated with them. 
a component of magnetization perpendicular to the surface, which would give 
rise to free poles on this surface. 

Such free poles would have associated with them such a large magnetostatic 
energy that the arrangement would be unstable. This magnetostatic energy 
can be eliminated in two ways. First, the directions of magnetization of the: 
fundamental domains can be rotated from the favoured [111] directions into two: 
of the three [110] diagonal directions lying in the plane of the slab. This would. 
result in small increases in the magneto-crystalline anisotropy energy as the [110] 
directions are less favoured energetically than the [111] directions in the 
boundary energy (approximately a three-fold increase) and in the magnetic field 
energy. Alternatively, the flux giving rise to the surface free poles could be 
closed through small prismatic closure domains magnetized in a direction parallel 
to the applied field. The latter suggestion is more in agreement with observation, 
for reversal of the field direction does not seem to affect the position of the patterns. 
The rotation of the magnetization vector involved in traversing the boundary 
between a fundamental domain and a closure domain is small and the wall energy 
of such a domain is also small. 

As can be seen from the photographs, the patterns formed on nickel are rather 
thin and diffuse compared with those found on cobalt or iron. For this there 
may be two reasons. First, the magneto-crystalline anisotropy constant of 
nickel is small (—5-1 x 10*ergs/cm® as compared with 4-2 x 10° ergs/cm*. and 
4-1 x 10% ergs/cm* for iron and cobalt respectively). Consequently a domain 
boundary is considerably thicker and the gradient of H2 (where H, is the stray 
surface field) is less and the attraction on the colloid particles is much reduced. 
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Secondly, the intensity of magnetization of nickel is less than that of iron or 
cobalt (500 against 1,700 and 1,400 £.M.u/cm'), and, as it is the component of the 
magnetization in a vertical direction at the domain interface which gives rise to 
the stray field H,, this results in a further reduction of the forces attracting the 
colloid to the domain boundaries. 

Unfortunately, because of the weakness of the stray fields existing at the 
surface of a nickel crystal, none of the techniques so far used to ascertain the 
directions of magnetization of the surface domains, e.g. the scratch, colloid 
striation or magnetic probe techniques, gave conclusive results. Hence at the 
moment one cannot safely decide which of the two domain arrangements suggested 
above is the more correct. 
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On Stueckelberg’s Treatment of the Vector Meson Field 
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ABSTRACT. Stueckelberg’s treatment of the vector meson field is investigated, and is 
found to involve difficulties similar to those of Fermi’s treatment of the radiation field. 
It is shown that these difficulties may be removed by the same procedure as has already been 
applied to the radiation field in an earlier paper, which involves the use of an indefinite 
metric and a modification of the supplementary condition. 


§1. INTRODUCTION 
T has been shown by Stueckelberg (1938) that vector mesons may be treated 
] in very much the same way as the radiation field. Stueckelberg’s treatment 
seems to have some advantages, but it involves the same difficulties as 
appear in Fermi’s treatment of the radiation field, according to which one 
cannot find a normalized state satisfying the supplementary condition, and 


+ VAS '. * State Scholar of the Government of India. 


696 Suraj N. Gupta 


therefore one cannot define the state of vacuum in an entirely satisfactory way 
(Ma 1949, Belinfante 1949). Recently the present writer (Gupta 1950) has 
shown that the difficulties in the quantization of the radiation field may be 
removed by using an indefinite metric and modifying the supplementary condition, * 
and here it is proposed to show that the same procedure may also be applied to 
Stueckelberg’s form of the meson theory. It must, however, be admitted that 
the difficulties to be treated in this paper are by no means as serious as those of 
the radiation field. For, instead of Stueckelberg’s treatment one can always 
use the usual form of the theory given by Kemmer and others (Kemmer ite 
which does not involve any such difficulty. 

We shall throughout follow the same treatment as has been applied to the 
radiation field in the earlier paper (Gupta 1950), which will here be referred 
to as I. Moreover, we shall consider only the neutral meson field, for the charged 
field may be treated in a similar way. 


§2. FIELD EQUATIONS FOR THE VECTOR MESONS 
According to Stueckelberg (see Pauli 1941) vector mesons may be described 
by means of two auxiliary fields A(x) and B(x) satisfying the field equations 


CEP —K*) Ane) = 0, aC" 2) B(x) = 0 ae eee (1) 
and subject to the supplementary condition 
0A, a 
(= +B) ¥=0. a ee (2) 


‘The field equations (1) can be derived from the Lagrangian density 


= 1 { (0A, 24 2 0B 2 B2 . 
L=~31(S) +A, +(=) eel, tees (3) 


which gives for the Hamiltonian density 


04,04, BOB 


= Ab Tdi th GDIMGIS TE a eee 2 (4) 

One also obtains for the commutation relations in the usual way: 
[A (e),4(4) = Aw aw) 9” Te eee 1)) 
L B(x), Be =A (0 oe eee (6) 


where the function A(x —x’) has been defined by Schwinger (1949). 
If we now define a field operator U,, as 
1 0B 
U, =A, seals K On, 65h b0 (7) 


the Sy esl condition (2) may be written as 


0U,, 


oe bo | a RT er seh oo (8) 
and the commutation relations for U,, will be given by 
eas I Gee ; 
U(x), U,(x')] =2 (3, ae ina) A(x =x!)torten aft RoR (9) 


In the present treatment we shall retain all the above equations, except that 
the supplementary condition (2) or (8) will be replaced by a weaker condition. 


On Stueckelberg’s Treatment of the Vector Meson Field 697 


§3. QUANTIZATION OF THE MESON FIELD 
Following I, we regard A, and B as self-adjoint quantities, and expand 
them as: 
A,(x) =X(2ky) 1? {a,(k) exp [i(k . r—Rot)] + a,'(k) exp [—i(k.r—Aot)] | 
k 


eet tO 
B(x) ee ek) exp [7(k . r—pot)] + '(k) exp [—7(k . r—Apt)] ) 


with Ree kee alt wait (EL be Hrinesneces (11) 
where a dagger denotes an adjoint. Substituting the Fourier expansions (10) 
in (5) and (6), we get 

last) a,\(k)=8tlo(k)aeb i(k) ald porreaes (12) 
Also, substituting the above Fourier expansions in (4), using the commutation 


relations (12), and omitting the zero-point energy, we obtain for the Hamiltonian 
of the meson field: 


JHdv = Eko {ay'(k)a(k) + a,'(kk)aa(k) + a5'(k)ag(kk) — a'(kk)ao(Kk) + 5"(kt)b(k) }. 


We now choose our zg axis along the vector k, so that the operators a,(k), 
a,(k) and their adjoints correspond to the transverse mesons. Moreover, for the 
present purpose it is more convenient to replace the operators a,(k) and a)(k) 
by a,(k) and c(k), which are given by 


a(t) = "2 a(e) ~ * ay(k), (Kk) = — aa(k) + “ Pa a (14) 

It may be verified from (11), (12) and (14) that 
: a (k)ay(k) — c(h) ck) = a5 (Is)a(k) — ay! (s)ag(k), eee (15) 
[a,(k),, a,'(k)|=1,~ lek), c'(k)]=—1, [a,(k), cl(k)J=0. -....: (16) 


Thus, the operators a,(k) and c(k) have properties similar to those of a;(k) and 
4o(k), and the Hamiltonian (13) may be written as 


[ Ho = ZRo(a,i(k)a,(ke) + a"(k)as(Ie) + a,"(ke)ar (Ie) — c¥(ke)c(lk) + 5"(k)b(K)). 


Besides the transverse mesons corresponding to a,(k) and a,(k), and the 
longitudinal mesons corresponding to a,(k), we have two more types of mesons, 
which we shall refer to as the b- and c-mesons. Since according to (16) the 
c-mesons have a negative commutator, we use an indefinite metric for these 
mesons as explained in I, while we treat the other types of mesons in the usual 
way. Thus, all the operators a,(k), a@,(k), a,(k), 5(k), and c(k) appear as 
absorption operators, while their adjoints are all emission operators. Moreover, 
now all the mesons have positive energy, though the states containing an odd 
number of c-mesons have a negative normalization. 


$4. THE SUPPLEMENTARY CONDITION 
It can be easily seen that according to the present treatment the supplementary 
condition (2) or (8) cannot be satisfied by any state of the meson field. Therefore, 
following I, we choose our supplementary condition as 
dU?) 
OK 
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Tey ae re (18) 
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where U,'+) is the positive frequency part of U,. It follows from (18) and its 
adjoint equation that the expectation value of aU (0x, Vanishes, so that our 
supplementary condition agrees with the usual one in the reveal limit. 

Substituting the Fourier expansions (10) in (18), and noting that we have 
chosen our z axis along k, we obtain: 


(i * a(k) —1%2 ay () + 5() e=0; ee (19) 
which, on account of (14), may be written as j 
{b(k)=c(k)}V-= Oe oh iy Ee (20) 


Since b(k) and c(k) are both absorption operators, (20) is evidently satisfied by 
any state not containing any b- and c-mesons. In general, as shown in I, a 
solution of (20) will consist of a normalized state containing no b- and c-mesons, 
and a superposition of a series of ‘redundant states’ with arbitrary coefficients, 
zero normalization, and zero energy. However, these redundant states do not 
give rise to any physical effect.* ‘Thus, for all practical purposes we may ignore 
such states, and regard the b- and c-mesons as entirely absent in a non-interacting 
meson field. 

Now the state of vacuum may be defined as that containing no transverse and 
longitudinal mesons, because the b- and c-mesons will be automatically excluded 
by the supplementary condition as explained above. Moreover, since the use 
of an indefinite metric for the c-mesons ensures that all the components of 4,‘+ 
contain the absorption operators, the state of vacuum will satisfy the relations 


ALLOY =0, BOK 5 =0) 30 OL =O “eee (21) 
It must be observed that if one does not use an indefinite metric for the 
c-mesons, the operator c(k) occurring in (20) will have to be treated as an emission 
operator, and it will not be possible to find a normalized state satisfying (20). 
But we have seen that no such difficulty arises in the present treatment. 


§5. TREATMENT OF INTERACTION 

For the treatment of interaction, it is very convenient to use the interaction 
representation, which has already been used by Miyamoto (1948) and Matthews 
(1949) for Stueckelberg’s formulation. We shall, therefore, add only a few 
remarks to indicate the points of difference. 

It should first be observed that in the interaction representation our 
supplementary condition retains the form (18), which may easily be seen to 
satisfy the usual requirements of consistency. This situation is different from 
that arising in the treatment of the electromagnetic field. For Bleuler (1950) 
has extended the theory, given in I, to the interaction representation, and he has 
shown that the appropriate supplementary condition for the electromagnetic field 
in the interaction Re is given by 


0A, Hox 

a al DO(x—x')j, (x")do, ‘| ¥te1- 0. dae (22) 
which is obtained by replacing A,,(”) and D(x — x’), occurring in the usual supple- 
mentary condition (Schwinger 1948), by their positive frequency parts A,'+)(x) 
and D(x x’) ). The physical significance of this difference in the two situations 
is as follows: The ‘charge term’ in the supplementary condition (22) allows 


* It should be noted that a redundant state will always remain redundant, because the normali- 
zation is conserved in course of time. 
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virtual longitudinal and scalar photons to appear in the presence of interaction. 
But, since there is no corresponding term in the supplementary condition for 
the meson field, the b- and c-mesons cannot appear even in virtual states. ‘Thus, 
these mesons play no part at all in the present theory, which is therefore 
completely identical with the usual theory, involving only transverse and 
longitudinal mesons. 

For the invariant perturbation theory based on the interaction representation, 
we require, in addition to the commutators (5) and (6), the vacuum expectation 
values of the anticommutators. Using the vacuum condition (21), we find by 
following Schwinger (1949): | 


({A,(2), Al2)}20=3- Ax), (B(x), BG} Y= Aa"), | 
({U,(x) U,(x’)}) =(8 wulll ihe ) ame —2") | cece (23) 
oe TH ; PN REO OX: ; 
It is to be noted that the use of an indefinite metric does not make any difference 
to the calculation of the vacuum expectation values, for the indefinite metric 
operator 7, defined in I, contributes a factor (— 1)” to a state containing m c-mesons. 
Therefore, for the state of vacuum 7 is just equal to unity. 


§6. LORENTZ-INVARIANCE 

In conclusion, we discuss the Lorentz-invariance of the present formulation. 
It will be observed that in I a preferred direction was introduced in our 
treatment, viz. the direction of the time-like anti-Hermitian component of the 
electromagnetic potential, and the Lorentz-invariance of the theory was ensured 
by the fact that the results of physical interest were independent of this direction. 
But, in the present case we do not have to introduce any preferred direction. 
For, one can split the vector field A,(x) into two parts, A,(«) and A,(x), 


such that 
1 0°A,(x) cat 1 0°A,(x) 
2 Ox,0%, ” A; (x) = A,(x) — oD) Ane ae sees s (24) 


1 — 
Ae) RTGS 
Following the interaction representation and using (10) and (14), one can also 
find separately the contributions of 4,(«) and A,@(x) to the Hamiltonian of 
the meson field. In this way, one can see that A,“(x) just corresponds to the 
c-mesons, while 4,(x) corresponds to the transverse and longitudinal mesons. 
Since one can separate the c-mesons from the transverse and longitudinal ones 
in a Lorentz-invariant way, the use of an indefinite metric for the c-mesons does 
not conflict with the requirements of relativity. 
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ABSTRACT. Using properties of states of quantized angular momentum, general 
formulae are obtained for the angular distribution in nuclear reactions and scattering 
processes in which the bombarding particles and target nuclei are polarized. A simple 
symmetrical formula for use in specific calculations is given. ‘The following cases are 
considered in detail : (a) polarized bombarding particles only, (6) target nuclei polarized in a 
magnetic field and bombarded by polarized particles. Finally, general predictions are made 
concerning the polarization of particles resulting from any reaction or scattering process 
and the results applied to the D—-D reaction which is suggested as a source of fast polarized 
neutrons and protons. 


§1. INTRODUCTION 
HE spin dependence of nuclear reactions and scattering processes might be 
effectively investigated by observing the effect of polarization of the target 
nuclei and bombarding particles on the angular distribution of the outgoing 
particles. 

Several methods have been suggested for the production of a beam of polarized 
elementary particles, the most successful being the polarization of thermal neutrons 
by scattering in magnetized iron (Bloch 1936, Hughes et al. 1948). The polari- 
zation of fast neutrons has been discussed by various authors (Schwinger 1946, 
1948, Lepore 1950), whilst Wolfenstein (1949) has proposed methods for the 
production of polarized protons. 

Partial alignment of the spins of the target nuclei can be effected at low temper- 
atures by direct coupling of the nuclear magnetic moments with a large external 
magnetic field, or, in the case of a paramagnetic substance, by using a smaller field 
merely to align the moment due to the orbital electrons. In this latter case, the 
magnetic field at the nucleus due to the electron shells is then usually sufficient to 
produce appreciable nuclear polarization (Rose 1949). 

The form of the angular distribution in the unpolarized case has been discussed 
by Eisner and Sachs (1947) using the theorem of vector addition of angular 
momentum and by Yang (1948) who considers the properties of states of quantized 
angular momentum under rotation and inversion. 

Wolfenstein has extended the former method and investigated the case of 
polarized bombarding particles, but, like the preceding authors, confines himself 
to obtaining general limitations on the angular distribution rather than obtaining 
formulae which may be used in specific cases. 

In this present paper use is made of the above methods in order to obtain a 
general expression for the angular distribution when both the bombarding 
particles and the target nuclei are polarized. This expression is then reduced toa 
simpler symmetrical form, use of which considerably facilitates calculations in 
specific cases. 

The formulae are then used to obtain more detailed results in two particular 
cases: (a) bombarding particles of spin } polarized in an arbitrary direction and 
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target nuclei unpolarized; (b) target nuclei polarized in a magnetic field and 
bombarded by polarized particles of spin 4, the axis of polarization being parallel 
to the direction of incidence. 

Finally, the general predictions made by Wolfenstein concerning the polari- 
zation of particles of spin 4 resulting from any reaction or scattering process are 
obtained in a form suitable for calculation in actual cases. As an example, the 
results are applied to the D—D reaction which is suggested as a source of fast 
polarized neutrons and protons. 


§2. POLARIZATION 
A spin function representing a state of polarization of a nucleus of spin s, will 
be of the form 
Pau a. where? $s |.U a1) Pile ee (1) 
m m 


m is the component of spin along the axis of quantization (z axis) and |U,,|? 
represents the probability of the nucleus being in the spin state 7". 

If the axis of polarization is defined by the spherical polar angles ©, ® and if 
the probability of the nucleus having a ‘spin orientation’ m’ referred to this axis 
is proportional to «,,, then since a spin function %” quantized with respect to 
this axis is related to %" by 


Pia ED (Dawe OOWM ere ek. Fea. (2) 


where the matrix elements Ds,,,, form an irreducible representation of order s of 
the group of real rotations of space, it follows that U,,, is given by 


by ( Sew (On) Dyn / ( Et i i OV Sioe 5 (3) 


€,/ 18 a Statistical variable indicating that for a given ‘spin orientation’ m’ there 
will be a random distribution of spins about the axis of polarization. It possesses 
the property (€n€,,) =Sin'm, Where ¢ ) indicates a statistical average. 

In order to describe the magnitude and direction of the polarization, it is 
convenient to introduce a polarization vector P defined by 

Reig ® (ret nn)\/ (St) rei (4) 
m m 
where r is a unit vector in the direction of polarization and P (referred to as the 
polarization factor) can be regarded as the ratio of the mean spin of the nucleus 
along the axis of polarization to its actual spin. 

Thus in the case of polarization of nuclei by a magnetic field H, «,, has the form 
exp (ym), where y=yH/ksT and » =nuclear magnetic moment, k = Boltzmann’s 
constant, 7’=temperature of the substance polarized. 

It then follows that (Rose 1949) 


P=[(2s+1)coth3(2s+1)y—cothdy]. ——...... (5) 
In most practical cases, y<1 and P then reduces to the form P=4(s+1)y. 


§3. THE ANGULAR DISTRIBUTION 
(i) The Nuclear Process 
For a nuclear process represented symbolically by A+ P->B+Q in which 
the target nuclei A and the bombarding particles P are polarized, the angular 
distribution of the outgoing particles in the centre-of-mass system can be written 


1(0$) =| © (2L+1%YU,V FRU OB)|, aes (6) 
1d ap 
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the direction of incidence being along the x axis. Here a, p, b and q represent the 
= components of the spins $4, sp, Sg and Sg of the nuclei A, P, B and Q, and Land / 
are the incoming and outgoing angular momenta respectively. yj"(r@¢) is a 
normalized wave function representing the relative orbital motions of B and Q in 
terms of the spherical polar variables 7, 0, 6 measured in the centre-of-mass system. 

No summation is taken over m, since, by conservation of the z component of 
angular momentum m=a+p—b-—q. 

FirrM is related to the probability of the nuclear process proceeding from 
an initial state f,={,? 47" to a final state be= Pe Pie 

Equation (6) has a similar form to that given by Yang (1948) except that since 
it is concerned with polarized incident particles, their spins are not orientated at 
random so that factors U,, and V,, dependent on the nature of the polarization have 
to be introduced and the summation over p and a kept inside the square modulus. 

Since ¢/"(r0d) involves 6 and ¢ in the form of a normalized spherical harmonic 
Y7"(@¢), it follows that, as in the unpolarized case, the maximum degree spherical 
harmonic that may occur in the angular distribution is 2/,,,,, subject to any other 
limitation that may be imposed. 

In Appendix I, by using the theorem of vector addition of angular momenta, 
it is shown that, in terms of the well-known transformation coefficients 


Foam = & Prat (Sas 43 Sp, P|S)(S,a+p; LM|J)(s, 6 +951 m|J)(sp, 05 sq, 418), 
Se 


where p7/” is a function of the energy of the incident particles and is related to the 
probability of the reaction proceeding through different eigen-states of the total 
angular momentum J of the whole system. S and s are the total initial and final 
spins. 

It is also shown in Appendix I from parity conservation arguments that 


—a—pL—M _ apLM , 
| ae cas te Pais Dn a Pe ny Qe O'S (7) 


only one of the signs holding for a particular process. 


(ii) Rotational Properties of I(6¢) 

Following the procedure used by Yang and obtaining a transformation 
equation for Fj?! ™.""(r0h) under an arbitrary rotation of axes through the Eulerian 
angles yo, 9%, ¢o it is then easy to show that 

HOS )=B| ALI yV yD xaPobs) Diy xB) Dla xeBbo) 
q Ma‘ap 
x Fram Pi" (rO$) P, 


where 6’ and ¢’ are measured in the orientated system of axes and are functions of 
6, d; Xo» A; do: 


If 6=¢=0, then 6’=6), ¢’ =7—¢y and, since #/"(rd) involves 6 and ¢ as 
Y7;"(6¢), it will only take a non-zero value when m=0, i.e. when M=b+q—a—p. 

By pone (2L + 1)Meih ReREo+a= ayy) =f, (8) 
and 6)=0, d5=7-—4¢, a final expression for the angular distribution is obtained: 

T(6$)== |X UyVypDi'y4q—a—v(9, 7 — $)Dy'p(9, 7 — $) Dita 9, 7 — 9) fon P- 


bq Lla’ap’p (9) 


Xo has been omitted since it only occurs as exp[—1(b+q)xo] and therefore dis- 
appears since the summation over b and q is taken outside the square modulus. 
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The maximum degree spherical harmonic that can occur in (9) is 
2(Linax +54 +5p) or 2(Lyyax +54 +Sp—4) whichever is even. This results from 
the form of the rotation group elements and standard parity conservation argu- 
ments. If the particles A and P are unpolarized, the maximum degree is 21), 
since, in this case, U,, and V,, are replaced by «,/(2sp+1)¥? and €q/ (2s, +1)? 
so that the summations over p’ and a’ come outside the square modulus and the 
corresponding rotation group elements disappear because of the orthogonality 


mm’ mm” — Pm'm”s 


property - SD! p* 5 


Using (8) and (7) an expression for f#%" is obtained and if this is used in (9), 
then two applications of the theorem given in Appendix II enable J(@¢) to be 
written in the following symmetrical form: 

CG) =D Uae DI pe (Osa) Gz, Paks TNO lenlt. wear (11) 


bq Jap 


where Hig CELL a, 4; sp, P| S)(S,a+p; L,0|J)p™*” 
and GL = 2 $7(700)(s5b5 sq gl s\(s,b+ 9; 1,0|S) piss, 


the factor p//5’ being written as p”” x p,,7. 

Consideration of (11) shows that the complexity of the angular distribution 
is further limited by the maximum value taken by the angular momentum J of 
any compound state, the maximum degree being 2/,,,,. In the particular case of 
the reaction proceeding through one compound state only, this state having 
angular momentum J =}, the angular distribution is spherically symmetrical. 
This arises because G, _,= + Gj, for a given J value (following from a similar 
argument to that resulting in (7)), so that on summing over 6 and g, all rotation 
group elements disappear by (10). 

If the reaction results in the emission of a gamma-ray photon, care must be 
taken in interpreting the wave function of the outgoing particle. Goertzel (1946) 
has shown that a photon may be treated as a particle of spin 1 having total 
angular momentum j (z component m) and z component of spin p, if the 
wave function ¥/"(r¢u) is replaced by A,,,(ré¢u). Here «= +1 corresponds 
to (A*+7A"),,, and ~=0 to A%,, the Cartesian indices referring to the 
' three components of the vector potential of a spherical wave decomposition of 
the electromagnetic field corresponding to a photon of total angular momentum 
quantum numbers 7 and m. 

Using this interpretation of the wave function of the out-going particle and an 
analysis similar to that already given, it is found that the same limits on the com- 
plexity of the angular distribution are effective. The limit 2/,,, arising in the 
particle case, however, is now replaced by 2). 

Thus, collecting results, the following general statement can be made about 
the complexity of the angular distribution of the outgoing particles in a nuclear 
reaction or scattering process in which the initial particles are polarized. 

If x is the maximum degree spherical harmonic occurring in the angular 
distribution, then the following inequalities hold: «<2(Lyaxt+s54+5p) if 
2(Linx +S, t+5p) is even, x <2(Lyaxt5,+Sp—4) if 2(Lmaxtsa+sp) is odd, 
# <QWJ max, X <Wnax (OF Zax in the case of gamma-emission). «=0 if only one 
compound state is involved, this state having J =}. 
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These inequalities should be compared with those concerning unpolarized 
particles. The latter can in fact be obtained from the above expressions by 
putting s, =sp=0. 

Calculations of the angular distribution in specific cases can be made using any 
of the general forms that have been obtained for J(6¢). In practice, however, the 
symmetrical formula (11) is found to simplify such calculations considerably. 

So far, then, general forms for the angular distribution have been obtained and 
limits on its complexity imposed. In subsequent calculations use is made of these 
forms and detailed results obtained for cases in which the bombarding particles 
have spin 3. It should be stressed that these results do not apply in general 
(i.e. for bombarding particles having spin greater than 4). 


§4. POLARIZATION OF THE BOMBARDING PARTICLES ONLY 
Suppose that only the bombarding particles P are polarized and that the axis 
of polarization is defined by the spherical polar angles © and ®. Then by (3) 


Uy = (2ey(%y-)"Dyry(—7, ~ ©, 0)) / (2). 


Since the nuclei A are unpolarized, V,, can be replaced by e€,/(2s, + 1)¥?. 
Then substituting for U,, and V,, in (9), there results 


1(,4)=Iuy(8)— Get sin Osin (=) E(RYOY*RyO)p eee (12) 
where R3?(0) = = fai Do, Sigh (C)yo. At eee (13) 


I,n(9) is the unpolarized angular distribution and P» is the polarization factor for 
the particles P. 
The following properties of Rj? are used in the deduction of (12): 


Eee) Sr) eit 1 Sage [(Rei(O)) *(Reg *) +¢.c.] =0. 
aog 
They result from (7),.and the fact that 


Dy le) — ( = 1 eA). sVehsielsle (14) 
Now if n is a unit vector defined by the vector equation 
kKjAk=nasin’, © 9° 4 Sree (15) 


where ky is a unit vector in the direction of incidence and k a unit vector in the 
direction of the outgoing particle, then (12) can be written in the vector form 


© 1(66)=Lyn(8) + 0- Pp E gr Vr(O)V7*7"@), wee (16) 
LLM 
where P, is the polarization vector for the particles P. g/4., is real and is defined 


M ait 4ar a-M—4 


pe oY eee , 4, L)«/ ¢b-+q—M-—H, —3, Ll’ 
Sub Isa #1) py {QL + 12k’ + ype ee aaa a Care er) 


In (16) use has been made of the property Di_ y,(0) = {4n(2L + 1)}12 Y¥(6), where 
Y¥(@) is a normalized spherical harmonic. 

Since n is perpendicular to the axis of incidence, it follows at once that 
polarization parallel to this axis produces no change in the angular distribution. 
Further, a reaction induced only by S waves still maintains a spherically sym- 
metrical angular distribution when the bombarding particles are polarized. 
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The maximum effect is to be expected when the polarization is perpendicular 
to the axis of incidence, for example along the y axis, in which case an asymmetry 
of the form 


T(Od)e= 1 ,(0); ees c0s ox Sine VE (0) VU (0) 9 o.. oac. (17) 
LLM 


is obtained as predicted by Wolfenstein who considers this particular case only. 
The results of this section may of course be restricted by any of the other 
limitations given at the end of §3. 


§5. BOMBARDING PARTICLES AND TARGET NUCLEI POLARIZED 
In this section the case of polarization about the axis of incidence is considered, 
the polarization of the target nuclei.being accomplished by means of an external 
magnetic field. ‘The polarization factor for these nuclei is then given by (5). 
The z axis is defined by © =0, so that by (3) 


Uy = € (tty)? / ( ay) | Vo se (By)? / ( =Bu Een (18) 


Also for polarization in a magnetic field, 8,,=exp (ya’) (see § 2). 
Substituting (18) into (9) gives 


1(0) = { 25 exp (ya)I,,(9) } / { 2s exp (ya) } Te See (19) 


where J, (0) is the angular distribution to be expected when the particles A and P 
have ‘spin orientations’ denoted by a and p respectively, and is given by 

I ,5(9) = 2s lees Dita fay) Dt, (0) DE(8) ie ae coy eee (19a) 
The rotation group elements have been written as functions of @ only, since, by 
the nature of the summation, it is obvious that ¢ is not involved and that the 
angular distribution is therefore symmetrical about the z axis. 

I,,,(8) is found to be a polynomial in powers of cos 6, only even powers occurring 
if L is confined to either all odd or all even values. Further, by (7) and (14), it 
possesses the property J,,(0)=_,_,(0). Using this property, (19) may be 
written 


1(0) = alcosk (ya) + Pp sinh (ya)]J,;(8) }/ { exp (ya) } fee. (20) 


where P, is the polarization factor for the particles P. 
In most practical cases y <1, so that it is of interest to consider the form of (20) 
for such low values of y. In Appendix III it is shown that to order y?, the angular 


distribution may be written 
2(Lmax +1) 


1(9)=Iy(9)+y % (apPet+yby)Cos"6, seuss (21) 
n=0 


i.e. the maximum power of cos@ involved is 2(L,,.x +1) independent of the spin 
of the nucleus A. As further powers of y are considered, then higher powers of 
cos @ will enter the expression up to the maximum permitted by the general 
limitations given in §3. 

It is at once obvious that if y=0, then the angular distribution is unchanged. 
This bears out the result of the preceding section that polarization of bombarding 
particles of spin 4+ along the axis of incidence has no effect on the angular distri- 
bution. This fact also implies that ifs, =4, alld, =0. 

If the bombarding particles are unpolarized, then Pp =0 and the change in the 
angular distribution is of the second order in y. In particular for a reaction 
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induced only by S waves, the angular distribution, which would normally be 
spherically symmetrical, has the asymmetrical form J(@)«1+Ay*(1—3 cos? 8), | 


ae a= {x (S49=54) 7,,0=0} /{31,,0=0}. 


For the case of a resonance reaction (one compound state only) induced by 
S waves and resulting in the emission of a gamma-ray photon, the following 
simple expression for J,,,(0) is obtained by using the concept introduced in §3 for 
the ‘wave function’ of the photon: 
Tn(8) £|(Sps Pi Sa @lF)P2{[ Darn, APPL 15 5p, A= L/P 
©* S41G; E13 sy AFT |S) PT 
Here J is the angular momentum of the resonance level and / is the multipole 
order of the transition resulting in the emission of the gamma-ray. 

Calculations for particular reactions using this expression for J,,(9) indicate that 
with y=4 and P,»=4, asymmetries of up to 25% might be expected. The 
observation of such effects however would in many cases be masked by gamma- 
rays from accompanying cascade transitions. 


§6. PRODUCTION OF POLARIZED PARTICLES 
(1) Nature of the Polarization 
Methods suggested by the authors referred to in §1 for the polarization of 
fast nuclear particles consist mainly of scattering processes in which there is strong 
spin-orbit coupling. Any nuclear reaction or scattering process in which such 
coupling plays a part might be expected to result in polarization of the outgoing 
particles, and in this section general results relating to the nature of this polarization 
are given. 
In the process represented symbolically by A+ P—B-+Q, consider the polari- 
zation of the particles Q (assumed to be of spin 4). If the initial particles A and P 
are unpolarized, the final state of the process can be represented by 


> TOT TOAD gag 2k + eat PREM HIE OIB). ooo (22) 
The polarization of the particles Q is defined by 
pt) P= (be Slit ot. Seseateunanes (23) 
where o is a vector representing the three components of the Pauli spin operator 
for the particles Q. (It is easily shown that this definition of P is consistent with 
that given in (4) when applied to particles of spin 4.) 
Substitution of (22) into (23) and use of the well-known properties of o gives 


; 2 
HOP (08) +2PLO8)= Ee S (BL + yeetReeNNeB)) 


x (saz+ IA eee gaan (t6)), 
Ll I 


1 
(2s, +1)(2sp +1) apo 


EAE + IEE Sime" (79) r) ; 


1(6) P.(64) = 5 (15 2 + 1) i Fentou(r64) P 
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where /(@) is the angular distribution of the particles Q in the centre-of-mass 
system. 

From (7) it follows at once that P,(0¢) =0. 

By investigating the behaviour of the quantity [P,(0¢6)+7¢P,(0¢)] under an 
arbitrary rotation of axes and adapting the procedure of § 3 to this case, there results 


HA). (08) +7P 08) = aay & (SDLRIO))” (ED* (OREO)) 


apb \q 
Be dee (24) 
where R;7(@) is defined in (13). 
Using the properties of R#?(@) given in §4, (24) reduces to 
1(0).P,(06) +iP,(04)] =iexp (id) D kM, VE@)YH(B), ce. (25) 
LUM 
where the coefficients ki}, are real and are given by 
— 21 z 4a 
ee aes ne ee Ce eee tepsinl pele 
SS Gt Cee CEANeL aye ee 


If a vector n is defined as in (15), then (25) can be written in the vector form 
TO) Risen ie HEU (O)V PL (0) ac gnraree ty. (26) 
LUM 


i.e. if the outgoing particles are polarized, the axis of polarization is perpendicular 
to the plane containing their direction of motion and the axis of incidence. Also, 
for P(@) to have a non-zero value, P waves at least must be effective in the reaction. 


(ii) Example: Application to the D-D Reaction 


Analysis of the D—D reaction (Konopinski and Teller 1948, Beiduk, Pruett and 
Konopinski 1950) shows that in order to explain the observed angular distribution 
up to deuteron energies of about 0-5 Mev., it is necessary to introduce strong 
spin-orbit coupling and P wave effect. As indicated in the previous section, these 
are the necessary conditions for the production of polarized particles and it might 
therefore be expected that the outgoing neutrons and protons are polarized. 

Assuming that only S and P waves are effective in the reaction (26) gives 
I(0)P(8) = n[(R8, — Rio!) YG(8) Y7(8) + (Ao, — RG) YY(8)Y4(0)]. If in the reaction 
the transitions taking place are 1S,>15S,, 3P—>°P, *P,—>1P, then in f}2” the 
following have non-zero values: $3}, pti» Pill» Pils Plot 

To conform to the notation used by Beiduk et al., the following substitutions 


for matrix elements are made: 


201 %ol’ =| Povo! o4| Bul? =5] pio 


Vt mm = 2%,,(L,0; 1, M|J)\(L,M—m; 1,m|J). 
J 


oles? =sleiyPs 


The energy dependence of the reaction is taken to enter only through the 
approach cross sections oy, (for definition see the references given above). 

In the above substitutions, M represents the initial z component of total spin 
and m the final ¢ component, so that matrix elements having Mm are those 
arising from spin-orbit coupling. 
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Then, with these substitutions, 


koi — Rig =0, Ay — RG = — 144/20, Im iors + Byy/+/2) + yovint$ 


1 |2 O [2 
and 1(8) =|29[?o {1 ee [ (en) cos? 6 
9% | xo| 
eS (“eh Seat) ane a} ; 
2| xo 


If J(@) is written in the form J(#) «1+ A(E)cos?6, where E represents 
the deuteron energy, then the expression for the polarization is easily shown 
to be P(E,#)=nBp(E, 6), where p(E, 0)=A(E) sin @ cos 6/(1+ A(E) cos? 6) and 
B= —4iIm[yio(yin + Bul V2) + Movin I/[4 lio? + 2 bol? — (Baal? — 2b’ P — 2byin* PA. 

The coefficient B is energy independent but cannot be evaluated without a 
knowledge of the relative phases of the various matrix elements.* If there is no 
spin-orbit coupling y}, =y;;' =f, =0 and it takes the value zero; otherwise to be 
consistent with experimental results for the angular distribution, its value must 
lie within the approximate range 1>B>-—1. 

The maximum value taken by p(£, 6) for a given value of E is readily shown to be 
Pmax(£) = A(E)/2(A(E) + 1)¥? at 6=cos1[(A(E) +2)-17]. Using values of A(E) 
obtained from fitting a curve of the form J(@)<1+ A(E)cos*@ to the observed 
angular distribution of outgoing protons for deuteron energies up to 0-5 Mev. 
(Huntoon ez al. 1940) we see that p,,,,(#) varies continuously from zero, at 
zero energy, to 0:54 at 0-5 Mev. At this latter energy, then, the polarization may 
take any value in the approximate range 0-50%. 


§7. SUMMARY 

From the foregoing work we see that polarization of the initial particles in a. 
nuclear reaction or scattering process leads to an increase in the complexity of the 
angular distribution of the outgoing particles. This increase, however, is limited 
by a series of inequalities referring to the maximum degree spherical harmonic that 
can occur in the angular distribution. 

If the bombarding particles have spin 4 and if only they are polarized, then the 
azimuthal variation of the angular distribution can be given exactly ($4). In 
particular, polarization parallel to the axis of incidence produces no effect on the 
angular distribution; also, a process induced only by S waves will result in a 
spherically symmetrical angular distribution for any direction of polarization. 

The target nuclei may be polarized by means of an external magnetic field,. 
and if they are bombarded by polarized particles of spin 4, the axis of polarization 
in both cases being parallel to the direction of incidence, then the resulting angular 
distribution is symmetrical about this axis and is in the form of a polynomial in 
powers of cus @, the maximum permissible power being 2[Linax +5,(+4)]. If the 
field is ‘small’, then the maximum power that may occur is 2[L,,,. + 1] whatever 
the spin of the target nucleus A (§5). 

The increase in complexity of the angular distribution is thus dependent on the 
form of the polarization. It is also related to the spin dependence of the nuclear 


* Note added in proof. Beiduk et al. (1950) have attempted to calculate these matrix elements: 
assuming that the spin orbit coupling may be attributed to tensor forces only, and to their 
approximation £,;=0 and y;;+ =—y?,, so that B is identically zero. However, their calculated 
results for the angular distribution show considerable disagreement with experiment and indicate 
that other spin orbit forces are probably acting. 
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forces that are effective in the process considered and should therefore give some 
information about the nature of these forces. (Bhatia (1950) has discussed in 
detail the information that could be obtained from an experimental investigation 
of the scattering of fast polarized neutrons by protons.) The multipolarity of 
transitions resulting in gamma-emission might also be deduced from this depen- 
dence of the angular distribution on polarization, but normally, experimental 
procedure would be complicated by accompanying cascade transitions. 

General considerations (§ 6) show that a nuclear reaction or scattering process 
resulting in polarized outgoing particles must necessarily be induced by at least 
P waves and must involve strong spin-orbit coupling, the polarization then being 
perpendicular to the plane containing the axis of incidence and the direction of 
motion of the outgoing particle. The D—D reaction satisfies the above conditions 
and it seems probable that the outgoing protons and neutrons are polarized 
although the exact magnitude of the polarization cannot be predicted without 
performing detailed calculations assuming some theory of nuclear forces. 
Measurement of the polarization, of course, will give information about the nature 
of such forces. 


AP EIN DL xed 
An expression for F¥2%™ can be obtained by applying the theorem of vector 
addition of angular momenta to an initial wave function of the form Wry, Be 
The product #3? can be expressed in terms of eigenfunctions J§*? of the 


total spin S, and the product 44°? in terms of eigenfunctions ” °s4*?*™ of the 
total angular momentumJ. ‘The indices L and S indicate the mode of formation 
of the eigenfunctions. Then 


W=X(s4,4; Sp, p|S)\(S,at+p; L, M[J)» Spgter™, 
SI 


where the coefficients of ” %s4*?*™ are the well-known transformation coefficients. 

The effect of the nuclear reaction is to transform the states % Sy¥4+?+™ of the 

incoming particles to Lpp/S? & °4+?t™ where ’%s7+?+™ is one of the states of the 
ls 


outgoing particles. An interpretation of the factors p47 has been given in the 


text. 
ow the state’ °4*”*™ can be written asa linear superposition of the individual 


states ¥7"(r, 9, >), a and Pro where m= M+a+p—b—gq, so that 
te= ¥ ph (6q,43 sp plS\(S,a+p; L, M|JV(s,b+45 1m) 
lsSJb m 
: x (Sp, b; SQ» qs) foes l (7, 0, $), 
ie. 


FalM — D pits’ (54,4; Sp, p|S)(S,a+p; L, M[J)(s,b+4; 1, m|J)(sp, 8; 59, 4|5). 
SsJ 


balm 
Also (S, M; s,m|J) possesses the property that 
(S, —M; s, —m|J)=(—1)§t*-7(S, M; s, mJ) 
so that Frp=f im = & (—1)"pins"(Sar 4s sp, p|S)(S,a+p; L, MV) 
a x(s,b+q; Lm|J)(sp, db; Sq, 9/5), 
where n=5,+5ptSp+Sg+L+l—2J. 
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But in a nuclear reaction in which the parities of the initial and final nuclei are 
fixed, it readily follows from parity conservation considerations that the integer 
quantity (L+/—2/) cannot have both odd and even values, and since 
5S, +5p+5p + Sq is fixed for the reaction, 7 can take only odd or even values but not 
both. Then, Fj? (5.4 = + F{ZkM, only one of the signs holding for a particular 
reaction. 


AIP REIN Dlx i 


Consider two states pe and pe quantized with respect to a set of axes 0x’, 
Oy’, 02’, obtained from the set 0x, Oy, 0z, by rotation through the Eulerian angles 


x, 9, >. 
The product state fF can be written as a superposition of eigenfunctions 


7+” of the total angular momentum J: 
pe a - 284, as Sp, [J )pot? meet ae Sp, p’ |.J) Dj, +2’, u(X 6, phys 


where yj’ is a state quantized in the initial set of axes. But the product state 
pe pr" can also be written as 


Ye Pp = 2Dy a X9$) Dyn xOP) PS bP, 


and then ff can be encased as a superposition of eigenfunctions of the total 
angular momentum J to give 


Pe Ve, = es Dil x9P)Dy'y(X94)(Sa,43. Sv, P|T V5"? 
Comparing the two expressions for pf 4?" and equating coefficients of }' gives 


(54,45 Sp, P"|J)Da +9, (x98) = ue (XO) Dy'y(x94)(8,, M—P5 sp, P|), 


which is the required theorem. 


Ded, BIND TEXS SU 


To order 7” (20) can be written 


18)= aq (1 PAE) Br) + y Pedal (0) + Dell ul) 


{1/(2s, + 1)}XJ,,(@) is just an expression for the unpolarized oe distribution 


I,,(6). Consider Bala’): By (19a) the summation over a involves terms of 
the form adi, )DiA(@) which can be written = (Der ‘(89))- a8 aq Dot vad)» 


since the matrices [D;4 (@)] are unitary and also real. This expression represents 
the elements (a,a,) of the transform by the matrix D/‘ (0) of the diagonal matrix 
Baa’ 

But D;\,(0) represents a rotation about the y axis through an angle 0, so that 
in terms of rotation operators the above expression becomes, using an obvious 
notation for Cartesian components of angular momentum operators and matrix 
elements, 
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{s,a,|¢~@a,6° | 5,44) = (4a, |a, cos 0 +a,8in6|s,a,), 
=acos@ if a,=a,=a; 
=4[(s, —a)(s, +a+1)]"? sin@ ifa,=a; a,=a+1; orifa,=a+1; a=a 
=0 otherwise. 


6 therefore enters XJ,,(@) through terms of the form 
a 


Ditpegaee oO) DS peat a—p(9)D} (0)D},,(8) cos 6 


$D1 


and DP 9 -a—p,(0)DEy4.¢-0-1-p,(9)Dbp,(6)Dt,,(8) sin 8. 


$P1 
By inspection it can be seen that both terms may be expressed in powers of 
cos @ and that the maximum power occurring is 2(L,,,.+ 1). 
The term Xa*I,,(0) can be treated by the same method and is also found to 
a * 


involve powers of cos@ up to 2(L,,,x +1). 
Collecting results then leads to equation (21). 
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The Scattering Phase Shifts of Two-Nucleon Systems 
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ABSTRACT. Energy expansions for simple functions of the scattering phase shifts are 
given for all relative angular momenta for central forces in neutron—proton and proton— 
proton scattering and for non-centrel static forces in proton-proton scattering including 
the coupling of P and F states. The differential scattering cross section for this latter case 
is given in an appendix. 


INTRODUCTION 

HE recent papers of Blatt and Jackson (Blatt and Jackson 1949, Jackson 
and Blatt 1950) have demonstrated the importance for the analysis of 
nucleon scattering experiments of the expansion of certain functions of 
the phase shifts as power series in the energy. Such expansions have been derived 
by a number of authors using various techniques (Landau and Smorodinsky 
1944, Schwinger 1947, Hatcher, Arfken and Breit 1949, Peierls and Barker 
1949, Chew and Goldberger 1949, Christian 1949, Bethe 1949). In this paper, 
the techniques used by Bethe (1949) are extended to the cases of general orbital 
angular momentum states and of tensor forces for both neutron—proton and 

proton—proton systems. 


§1. NEUTRON-PROTON CENTRAL INTERACTIONS 
The radial wave function multiplied by 7 satisfies the differential equation 


sos [ e— SY - P(r) |y=0 Bre (1.1) 


if 


with the boundary conditions (0, k) =0 and 

y(7, k) ~ A,[sin (kr — l/2) + tan 8 cos (kr —Ir/2)].....-- (1.2) 
In the above P,(r) represents a central, short range nuclear potential for states 
with angular momentum /h. The remainder of the notation is standard, 8 being 
the corresponding phase shift. The subscript / will in future be dropped where 
no confusion arises. 

Expressions for the phase shifts are readily derived by converting (1.1), along 
with the boundary conditions, into an integral equation. Let vp and v, be 
regular and irregular solutions of (1.1) in the absence of the nuclear potential, 
oe 

Up =hkrj,(kr) ~sin (kr —In/2), vp= —krn,(kr) ~cos(kr—In/2). ...... (1.3) 
The Wronskian of vg and v; is W=—1. ‘Then, by variation of parameters, 
it is readily established that 

1 7 CO i 
y(r) =Aeg(r) + aps) | a()POD() ds+ex(r) | en(s)PON) as} (1A) 


~ Asin (kr —/7/2)+cos (br —In/2) 5, | Only ds. nro? (155) 
0 


* Now at Atomic Energy Project, Chalk River, Ontario, Canada. 
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Comparison of (1.2) and (1.5) gives 
kWA =cot 8 i optt)P(y(ndr. scan (1.6) 
0 


The normalization of y is still arbitrary. A formula for the phase shift 
independent of the normalization may be readily derived. Substituting the 
value of A from (1.6) into (1.4), multiplying through by yP and integrating 
throughout the range of 7 yields 


J Pens?) ar— |" arP@yr) | G(r.) PC) ds 


ial maa RW niece J) ets (1.7) 


G(r, s) is Green’s function given by 
RWG(r, 5) =Vp(r)0x(5), Sac; 
=0,(r)v,(s); s>r. 


This is the Schwinger variational form. Other expressions for the phase shifts 
may be deduced by specializing the normalization. A convenient choice, 
particularly if y is determined by applying the iteration procedure to (1.4), would 
be A=T. In this case by (1.7) 
av]. UpPy dr. ye (1.9) 
The expansions for the phase shifts may now be obtained using relationships 
of the type (1.7) and (1.9) directly. Alternatively, Bethe’s (1949) methods may 
be generalized to the case of general orbital angular momentum. We derive 
below the expansion for k?’+1cotd, using Bethe’s techniques. In this instance 
we choose A= cot6 so that the solution of (1.1) with P=0 is 


v=coto[vpg(7,k)+tandv,(7,k)]  — ...... (1.10) 


and for large 7, y is identical with v. 

Let y=y(r,k) and yo=y(7, ky) be the two solutions of (1.1) corresponding 
to the wave numbers k=k and k=ky. The corresponding solutions (1.10) are 
then labelled v and v). Then, following Bethe (1949), we have 


[229 — %y0']p, = (H—Re?)| (Vo 2%) dr nee (1.11) 


ro 


~tand = 


where 7, implies the limit y>0. Substituting for v from (1.10), equation (1.11) 
reduces to 
— RAR *cot 8 + kok" cot by 


= (k? — k,?) | (Y¥q — VUp) 47 — [Uy —Vyp?'],,- se eee cll Pay 


Noting that s [Wyo —VypPz | =(R2—Ro2)Vip, wee we Gi3) 

it iS easy ‘to see that the square bracket in (1.12) just cancels the singularities 

arising in the integral. We can then write 

R2z+1 cot 8 = Rot cot 65 + R'hy (Rk? — Ry?) ie {(vUp — Yo) — Singularities at.r =0} dr. 
ok er (1.14) 
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We now expand 


Y(r, k) =R'y(r, k), Vins k\=R OR) meee (1.15) 
into the power series 
VY =YV p+ (R?— Ry) VY +(R? Ro? )P Vote ces | ene ses (1.16) 
V =Vo+ (R= Ry?) V, (Re Re? )PPVagt.... ea ae (1.17) 
Substituting into (1.14), we have 
R2+1 cot 8 = Bo(Ry) + (A? — Ro2)B1(Ro) + (KR? — Ro?)*Po(Ro) +--2. ee eee (1.18) 


in which B,) =k,"*! cots, 


9 : me er eer (1.19) 
B,= | (VoVin-1— YoY n-1 Singularities at r=0) dr 
10) 


Bo is thus determined by the asymptotic form of yy which must be obtained, 
perhaps numerically from its differential equation, perhaps by iteration in a 
formula like (1.5). 

From (1.10) and (1.17), we have. 


Ve= Boker Sono hoot oe a eran ra (1.20) 

Vi=[ BV) | =the ent Bol ae gt |, + [ae | 
pe FP nolo Sa Po liye: peel CB |e een eee 
Pe Bibs (1.21) 


i aeds f 
Vo=3| qe |, =Paltdkr tno +B Lae], 


a? 1p @ 
vi ____ spl a3 Asien yy) 
+ 2Bo baer 5 ex | nag) E aye" vs | on 


0 


eige (1.22) 
etc. The coefficients Y,, in the expansion (1.16) may also be defined as 
lelerd: 
Y= VG Ra) Mo ay Ee: (| Fe ee (1.23) 
or alternatively through the set of differential equations 
aY Wi+1 
—! + se — -P() | Vi (in ace (1.24) 
dav Wi+1 
a + | Ae oe - P¢*) | Vue. YS Aes (1.25) 


as may readily be established by substituting (1.16) into (1.1) and equating the 
coefficients of each power of (k?—,?) to zero. Keeping in mind that Y, ~V,, 
one can also write the solution of (1.25) as a solution (perhaps by iteration) of 
an integral equation, namely 


Yn=Brko pot (L/W) EE i Vo(8) P(s) Y,(s) ds 


ee, f vno{ P(s)¥,(8) — Yn-a(8)} 45+ Opp I. v10(8) Yna(s) as| wae (1.26) 


The most useful energy about which to develop these expansions is zero, 
in which case kj =0. From the series expansion for j,(kr) and ,(kr) we deduce 
that in this limit 

Vo=PoByr'1+[((2l+1)Br'y*  —,..... (1.27) 
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where 1/B,=1.3.5. ...(2/+1). Also Y4 is that solution of 
GY fare [Kl Eli BP Y =0h a eee (1.28) 
which is asymptotically equal to Vy. The higher order functions V,, and Y,, are 
obtained by taking similar limits. 
Since 8, has the appropriate dimensions, we can write By = —a,*"1, where 
a, is a generalized scattering length for zero energy. Then 8, may be written 
B,=4a,?'r,, where r,is the ‘range’ of forces in the state with angular momentum 
Ih. Thus for ky =0, (1.18) becomes 
aP'k?+1 cotd = —1/a,+4r,k?—713PR44+7,°O Ro—. 1. ween (1,29) 


In this notation (1.19) may be rewritten 


ro=2| [(1 + r/a,)? — Y,2] dr, eels) 
T,=2 ii [Bray era (21 Ay at Vee dr, gp UO, ie ie eo crane (1.31) 
Po= { | desta — Syn? + 4 + [up)x? — da0+ “ YoY, | COR Cae eae (1.32) 

0 0 j 

? 1 

aoe ip | abet? + Rant +E a? — Qu? a ac YoY, | dx, 
Phe ee ee eRe (1733) 
P,= ( [Pt B Px?! +4/2(21+ 3) + $B Pur tat? + Qu .x3/(27— 1)(21 + 1)(274+ 3) 

0 

—4%/2(214+1)+(a?/r?) YoY] dx, I>1, ei arate (1:34) 


where 14, =7;/4;. 


Before leaving the neutron—proton system, it is desirable to consider the 
normalization dA4=1. Then 
y ~sin (kr —Ur/2)+tand") cos(kr—Im/2) (1.35) 
and tan") is given by (1.9). We see that as k>0 


yo= —lim k?“ tan 3 > | " Br P(y(r)dr, 1... (1.36) 
k->0 0 


where y,(r) is the solution of (1.31) with the asymptotic form corresponding to 
(1.36), viz. 
Ven fy (2 Br 7 tees (ho?) 


Alternatively, yp is the solution of 


yo(7) = Britt —ritt ia s~'P(s)yo(s) ds/(21+ 1) 


= rt |" st4P(s)yo(s) ds|(21+ 1) pefannt less) 
0 
The phase shift for general k may be expanded about k, =0 in the form 
kl tand\) = —(yot+ yk? +yokt+....). hon i iGlsso) 
Then we | eee ay xy) Pdn  ...+.- (1.40). 
0 
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where x,, and y,, are the coefficients in the series k-’1y = yo + y,k* + yok*+.... and 

KN = ty + Re + gk o5 Le, Hy = (—) tere ee (32a ee 

from (1.3). 
§2. PROTON-PROTON CENTRAL INTERACTION 


The radial wave function times ¢ satisfies the differential equation 
ze 5 + (A? —Ul+1)/r—2ak/r—P(r)y=0 wea (250) 


and the boundary conditions 
¥(0, k) =0, y(r, k) ~ A, [sin {kr — l7/2 +5 (r)} 
+tan 8) cos {kr —I7/2+5(r)}], «0... (2:2) 


where «=1/2kr, R=h?/Me?, 2«k/r is the Coulomb potential, P(r) is the short 
range nuclear potential, 5(r) the Coulombian phase shift 
8 O(r) = —« log 2kr, yOruargT(ta+l+1),  . 20. (2:3) 


5 is the nuclear phase shift and A, is fixed by the normalization of y. 

If we denote by F', and F; the regular and irregular solutions of equation (2.1) 
in the absence of the nuclear potential, with the asymptotic behaviour indicated 
below, then (7, k) is pipseucee the solution of the integral equation 

y(1, k) =A, F2(7, eee Te A) | F(s, k) P(s)y(s, k) ds 
Se BGR) i Fils, Py kd) 9 eee (2.4) 
a 0 ; 


~ A,sin {kr — [7/2 +8 (r)} + cos {hr —bn/243()} 


x { Fa(s,k)P(s)(s,k)ds, Rr>1. cae ene (2.5) 
0 ‘ 
It follows in the same way as in the n—p case that 
cote” (= | 
Ave aia Fark) P(r) yr ara (2.6) 


and for arbitrary normalization 
| ~ Pe)p*{rk) dr | y(r, k) P(r) dr | G(r, s) P(s)y(r, 8) ds 
( 0 0 
| [Fal POM ® ar | 


To derive an energy expansion for the phase shifts we again use Bethe’s methods. 
We choose for 4, 


cot 8? =kW ee) 


A, = C,cot8™, hope Nee (2.8) 


in which C, is the Coulomb barrier penetration factor, dhol: tO be defined 
explicitly. es let F(r, k) be the complete solution of the auxiliary equation 
such that F(r, k) is identical with y(r, k) outside the range of nuclear forces, i.e. 


F(r, k) = C, cot 8[Fa(r, &) + tan 8 Fi(r, h)). 
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It then follows as in the n—p case that 
LFF,’ — FyF’],, =(# — he?) le oe FRGdan (2.10) 


o The solutions Fp and F, have been iventicated’ in some detail by Yost, Wheeler 


and Breit (1936). Using their results and substituting for F from (2.9) into 
(2.10) we have 


R241C2 cot 8 + ory a: [(21+ 1)(q—plna) +p] =g,(R, Ro) + R22 1C,O” cot 8 
+ kg?! *1[ (21+ 1)(Go — Po In. %) + Po]/(2/ + 1)? 
By Cole A oted (1) a a a i ae ene: oe (2.11) 


in which g,(k, ko) is a known function of the energy: 
81(R, Ro) = —h'hg' (21+ 1)?{(21— 1) a, R'a_y 44 Rg 
+ a,_ yk 1a_,,5\RK'42(21—3) + 2... ayy hk (1 —2Da, ky}, 


Same (2.12) 
(20+ 1)! C,=2[02 +02 ]¥2[(Z— 1)? + o2]¥2.... [1 +2] ¥2(27rx) N2(e?™ — 1) 2, 
A Ps (2:13} 
Se eer Cer ee yd be eters (2.14) 
7 l 3 21+1 i hoe 
g=p {3 si +22) 2 s-1+Re 135) 
; d 2"(ia+n—1)...(ta—/) 
a) Be) ae coe TT 1, ene (2.15) 


and a, are coefficients defined by Yost, Wheeler and Breit (1936) by certain 
recurrence formulae. In (2.11) 


PAR, ky) = if [FF —yyo—singularities atr=O]dr. ...... (2.16) 
0 


The singularities referred to in (2.16) arise in a manner similar to those in the 
n—p considerations; that these singularities are subtracted out may be readily 
established. 


Explicitly, for /=0, formulae (2.11) and (2.12) become 


§ 2 
! al Ree OF ; +h se) | = 5 al ae sag) i+ hla) | + (R? — Rg?) p(s Ao) 


exp (27a) — exp ( (Ore 
ees (2:17) 
£o(R, Ro) =9, and ho(«) is defined by 
Orel (te) 
ho(«) =Re T(—2) —Ine.\ Sane adhe 
For /=1, l=1, this formula is 
nL oe EeEGs)COUOU Un ahp9?(1 + & 9”) iS, a) Ra 
R al Spee exp (277%) — 1 gael) |- 2 oe exp (27a) — 1 ay hy( 4%) | 
+ = (Rk? — ky?) + 9R?y?(R? — Ro”)py(R, Ro)y ww we ee (2.19) 


£3(k, Ro) = ae (Rk? — ky”), and h,(«) is given by 
** hy(&) = —1+(1+a?){Re I” ( —ta)/T'( ~ia) —Ina —3/2}.° 2... (2.20) 
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The expression for g,(k, k9) has also been derived, and is 
2580(k, Ry) = 13(R? -- Rp?)/576R? + 3(k4 — Ry')/96R. ...... (2.21) 


We now perform the expansion of (2.11) in the special case ky=0. For this 
purpose we denote 


lim Rao C2 Coto = 4+") ee eee (2.22) 
k>0 
and also write 
lim RR?"[(21+1)(q—pln«)+p]=h,(R). sae (2.23) 


; k>0 
In particular 


ho(R) =lim hy(«) =0, h(R) =m hie) = —3/38R2 7 (2:24) 
k>0 k+>0 
Similarly, let 
lim g,(k, ko) =2,(k) =yiR?+yoktt....yjR wees (2:25) 
k>0 
where the coefficients y,, are defined explicitly for 7=0, 1, 2 by (2.17), (2.19) and 
(2.21) respectively. Choosing also an expansion for k'F(r,k) and k'y(r, k) as 
RUF(r, k) =Xo(r) + 2X, (r) + BAX, (r)+.... | 


vta-20) 
k'y(r, k) =Vo(r) + R2y,(r) + R4y.(r) +...-, 
we have, on substituting into (2.11) 
RAC 2 cots + k?'+1/(21 4 1)?[(27+ 1)(q—plna) +p] 
=a +h R)/(21+ 1PR + (7, + By)R? + (v2 + Bo)R* + 
(yi B Re Bak? Bea Wee (2a20) 
in which he A (XE Kul Shay ey al? ei eee (2.28) 
The coefficients X,, are now defined in more detail. We recall that 
REG, R)=C Pkt coral (7, k)iC,k 1 PRC ai Rye (2.29) 


_In Appendix A the power series development for F(r,k)/C,k'+1 and 
k'C,F\(r, k)is outlined. Using these developments together with the development 
for C?k?'*! cot d given by (2.27) we get 


Xr) =a72-1S%(r) + NgWM(r), see (2.30) 


ee eae oy Oe te “t (ae ory 
x {(214+-1\(q—pina)+ Pl) pS Di) NO eee (2.31) 


Xmer9,09) + {n+ Be~ (om pr ipl(2l+ Nla-plna) +21) 


x S00) + (ret Bam (aa or pl@l+ Ng—pina) +1) 


x SyO(r) + No(r) ...2+- (2.32) 
etc. 


The coefficients y,, may generally be written down as solutions of either the 
set of differential equations 


= 5 1+ 1) /72 4+ Zakir + Pr) ly, = 9,1, WSL sneer (2.33) 
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or the integral equation 
1 (oe) 
y =A OSM Ze. 8 (2) 
Yn a So (r)+ RWS, N) So (7) I. No (s)P(s)y p(s) ds 


os ams: mye) J So?(S)LP(5)¥n(S) —In-a(s)] ds 


1 ; 
a YY) (0) 
+ EMS, Ny 8) [No OS HOCH VE Seer. (2.34) 


in which 


k=0 


A, =7,+ Bu - Fae [(20+1)(q—p Ina) +p] } Sa (2.35) 


§3. NON-CENTRAL PROTON-PROTON INTERACTION 

We derive the contribution to the triplet differential scattering cross section 
taking into account the P state and the coupling of the F state to the P state, 
having assumed a tensorial (axial dipole dependence) interaction. The relation 
of the phase shifts to the interaction potential in this instance is also derived 
and the expansions in the energy are also indicated for quantities closely related 
to the phase shifts. 

In the non-relativistic Schrédinger equation for the wave function of the 
two protons in the centre-of-mass system, viz. 


[AS GM iE —V )\ib=O, ee iets (3a) 
we choose an interaction of the form 
V=JS(r)+ K(r)o 624+ L(ir)Ditelr we aes (3.2) 
in which D,, is an axial dipole coupling 
Dip So. xoa?) x.) I Dio ws ais (3.3) 


The notation used is standard; J(r), K(r) and L(r) are radial functions, « the 
Pauli spin operator for particle (1), etc., x) a unit vector in the direction of the 
line joining the two particles, e?/r the Coulomb potential. 

The states of the system separate, in the usual way, into two classes: the 
singlet states for / even, triplet states for /odd. In the singlet state, J?=/?, J,=L, 
and the complete spin, angular momentum and radial dependence of 1% for a 
given j and m is 


LU EN A IEG Cay Aan (oR air ae Nee ee (3.4) 
in which 1R,‘”(r) is the solution of the equation 
[d?/dr? + k? —I(1+1)/7? —2ak/r—1 P(r) PR O(r)=0 «00... (325) 


and, for an interaction of the form (3.2), 
IPO) =(MIMYIG)—SKGS he ws (3.6) 


In the singlet state, then, the definition of the phase shifts, the energy expansions 
and the differential scattering cross section are identical in form to the case when 
the potential (3.2) is purely central. 

In the triplet state, pertaining to the possibilities 7 =/ and j=/+ 1, for J odd, 
the wave function 35,” may be generally written as 


3h m= 21 /rPRiO (1720 ueeene (Sta) 
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in which the summation is /=j for j odd, /=j+1 forj even. The radial wave 
functions *R,'” are now solutions of the equations 
Vea: (d2/dr? —j(7 + 1)/r? 2ak/r+ kh? PR 
Tia A Ras eho usa ey On (3.8) 
and the coupled set 
d* (j-1)j 2ck 


f=j—1: 55 P(r) 
, Ais)) , G+ Ir" & J steate 
| 1 eae 9 ) |u,+ PEs Ce ei y ebonck (3.9) 
fjords a GHING+2) Dak pe Pr) 
alas 2) , G+ Dr" Bs 7 | 
aco Or ) e+ ar oy SS bar O(r)ju,=0.0 a8. (3.10) 


In the above we have set 
u,= 3R,'79, 0; = RO), P(r) =(M/#)[J (7) + K()], 
OPS (MRPALG = = eee eae (3.11) 

We note here that for /=7, j odd, the radial wave function is uncoupled ; 
the definition of the nuclear phase shift and the expansion in energy may again 
be performed as outlined in § 2. 

As for the coupled set (3.10) and (3.11), we note first that the solutions u; 
and wv, are to be regular at the origin, i.e. u,(0,k) =v,(0,k) =0, and outside the 
range of nuclear forces they are equivalent to the solution of (3.10) and (3.11) 
with P(r) = Q(r)=0, having the asymptotic behaviour 

u;~A,(cos6,“)* sin (kr — 47-1) + 8! V(r) +8), 

‘uv; ~B,(cos6,”) sin (kr —4(7 + 1)7 + 8 +) (r) + 8,”). 
Furthermore, we see that there are four linearly independent functions, regular 
at the origin, satisfying these conditions. Call these u,™” and uv,” and the 
corresponding nuclear phase shifts 6;““%” and 6,”%”. The triplet wave functions, 
defined by (3.7), will be labelled 35 ,%™ and °%,”™, viz. 

ti tae = (1/r)[u Vie. g-1)m die Ue A 2. 
etc. | 

By a slight generalization of the phase shift concept, we rewrite (3.13) and 
its counterpart with « replaced by y, so that 


Bh 0% ™ —~ (1/7 cos d,*) sin [kr —(j — 1)7/2 4+ 8-(r) +87] 


7 OZ Rico Zee 2 eee (3.14) 
3h ™ ~(1/rcosd,7) sin [kr —(7 + 1)a7/2 + 87 +Y(r) + 8,7] 
187 Me ay eZ Oe a (3.15) 


w,~ and w,Y measure respectively the extent of the admixture of the state with 
orbital angular momentum number /=7+1 to that state with J=j—1 and 
vice versa. We are here following the methods given for neutron—proton tensor 
forces by Rohrlich and Eisenstein (1949). By substituting (3.12) into (3.13) we 
find that (3.14) and (3.15) may be arrived at by setting 
Oj 5 FHS 4 fF YD _ AI (3.16) 
5,” ee Oey =O," aie Ftd) = yD) +7 
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and setting 


a 
G) fe oe cos 6,*[cos (8%; +yFY—yPD—a)Pt (3.18) 
j 
7 BS Y £f54+1) _ (3-1) yp 
wi = Fcos Oy ee TH COSO, | te Pee (3:19) 
iy 
It also follows from the orthogonality conditions on the wave functions *;%™ 


3,/, Y,™m 
and *f,%™ that Cea = eee ee. ose (3.20) 


The total triplet wave function for a fixed 7 and m, 7 even, is then a linear 
combination of the functions °4,%™ and *3,%™. Call this 34,"@). If we denote 
by *%, the triplet wave function for fixed 7 and m, } odd, then the total triplet 
wave function for a fixed polarization, *y,,,”, is 

Prot” = 2 ses et Sine na aes deel Wee wi aaras (3.21) 
jo 


jeven 
The total unpolarized wave, ee is obtained by suitably averaging over m. 
The scattered wave may then be written in the form 
Bh oat ~rocat, ae (1 /2ir) eilkr—Inj2—alog2kr) { p,m 8 Z (om exp (yh Dee 28,) 
m, l odd 

+ [ogy exp UY? +.28744%) — Yaga exp Uy + 28,117) ] PZ). 40" 

+ [ype 1” exp Uy + 28, 1”) — yey. 1% Exp Uy + 28,_1°)] 224 O™ 

SiR) jar oyna cll) ur y exp Zyros (3,22) 
In (3.22), *bcat,¢ 18 the scattered wave when purely Coulomb interactions are 
considered. The constants «,;”, 8;”, y,;” are fixed, in the usual way, by demanding 
that asymptotically *y4o, —*%tot,- be an outgoing wave. The constants *7,, are chosen 
(for unpolarized waves) |*/,, |? =1/4 for m=1,0, —1. The phase shift 5,corresponds 
to the radial wave function /=j, 7 odd (see equation (3.8)). The differential 
cross section in the triplet state is then 


dS, = Boat D Ph aee 1 dain We BI deca (3.23) 


The cross section has been derived taking into account the states °P,, °P,, 
8P, and the 9F, coupling. The formulae are quoted in Appendix B. 

We turn now toa further inspection of 5,”, w;*Y. For this purpose we convert 
the differential equation (3.9) and (3.10), together with the boundary conditions 
on u,, v; into an equivalent integral equation, namely, 


uy AyFy,j 1 (UR)Fu, ea] Fa, wbsls) as (1/R) Fr, a] Fr,» xfsl9) ds 


where Fy.) Fy 5-15 Proj, F, j41 are the Coulombian functions and are the 
independent solutions of the auxiliary equations (3.9) and (3.10) with P=Q=0; 


hye [ Pe) - 3? Os ©) | — TYEE 200, he (3.26) 


Lwe y= [ Pi) sy 219) | v,- we OC = ieee (3.27) 
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Comparing now the asymptotic forms of (3.24) and (3.25) with (3.12), we have — 


A,=— cots ip anata) — Mee (3.28) 


B,= —cot8,"| Dil oda he eee (3.29) 
0 


The phase shifts may be conveniently defined as follows. In (3.28), choose 
A,=1. Then 


Cy. ?k2/*3 tan 8,¢= — | "Cy Fe, glk IC, 1d (s)] ds, oe. (3.30) 
0 


where C,_, is the normalization factor for the Coulomb wave functions. A 
similar definition for tand,” may be arrived at by choosing in (3.29) B,Y=1. 
These choices of the normalization constants would be convenient if equations 
(3.24) and (3.25) are solved, for instance, using the iteration procedure. 

The energy expansions for quantities closely related to the phase shifts may 
also be derived by generalizing Bethe’s technique. For this purpose we choose 


bi = Ge .coto,.. « Ne oe 6:31) 
Let u, Ug; v, Up be the two solutions (3.9) and (3.10) corresponding to the wave 
numbers k and ky, and F,_,* the complete solution of the auxiliary equation 
(i.e. (3.9) with P(r) = O(r) =0), viz. 
Bigeye COU0 * Es goa per le meee (3.32) 
It then follows, as in §2, that 


6lj9 + I)” 


[FF,! — FyF’],, =(# —k,?) | (uu) — F Fyldr + ar1 


To 


i O(Nou, —uopldr. 


In the latter derivation we have dropped the subscripts « and 7 where such 
omission causes no confusion. . 

Comparison of (3.33) with (2.10) indicates that the results are similar in form, 
though the expression (3.33) is of a more general character. The derivation of 
the phase shift and its expansion in the energy may be achieved in a manner 
outlined in §2. 


ALP PENG TD XA: 


ENERGY EXPANSIONS OF THE REGULAR AND IRREGULAR 
COULOMBIAN WAVE FUNCTION 
The functions S=[CA1) FOr, k) and N=C,R'F(r, k) are two linearly 
independent solutions of the equation 
a’G 
2 za 


7 Ki+1)jP+URrG=-kG. =a... se (A 1) 

‘Their expansions in the series 
S,=S%(r) + RES) + RASH) + AD 
N= NO(r) + RN Or) + RBNO(r+ oe. fo ee 


have been investigated by Yost, Wheeler and Breit (1936) for S, and recently 
by Breit and Hull (1950) for N,. The terms of these expansions have been 
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derived explicitly by other means for certain particular cases, and in our notation 
are, for all /, 


S§0(r) = RAT (214.2) (r/R) lara, ae (A3) 
N(r) =2R-4 1 /P(21-+ 2) (r/R)? Kora 3* 
for /=0, 'S,%(r) =(— 1/3) RE, } 
S.(r) =(1/18)(7/RY!2 R21, — (1/45) R271, Conk (A4) 
Ny(r) = (2/3)? Ke, 
N,(r) = (2/45) RPK, + (1/9)(7/R)!2RPK; ; | 
for J=1, S,(r) = —2P-R°[I, + (R/r)#2I5], } 
Sy(r) = (1/3)r*R[(r/R)PPT, + (4/5)(r/R)!L, + (8r/5R)Io, oa: (A5) 


N,(r) =(1/9R)r?2.Ky + (2/9)(7/R)2-K,, 
NM (r) =(1/54)r°[(7/ RPK + (12/5) Ka]. } 
In all of the above the argument of the Bessel function is 2(7/R)#?. 


It will be noticed that in the limit as Ro (i.e. the n—p case), the function S, 
reduces to the expansion of 


[(orkr/2)"2I 4 o(Ar)]2+D (1+ 3/2) /arb2keH4 =r? +[ 1 — 722 /2(2143) + ...] ...(A6) 


It is readily checked that the terms as defined by (A4) and (A5) conform to (A6) 
in this limit. On the other hand N, does not reduce directly to the corresponding 
irregular solution of the n—-p system in this limit, as may be surmised from 
inspection of the series defining F;(7, k). 

Blatt and Jackson (1949), in their consideration of the S state, have developed 
an alternative expansion for Ny which is slightly less restrictive than the 
expansion given in(A2). The essence of this expansion is to separate out of the 
irregular solution a part which reduces to the solution of the n—p system as 
R+>0o, Le. 


U(arkr/2)#2I_y_1)o(Rr)\(— 2-12 D+ 8) = (204 1) — F2r?/2(1 — 21) + . «J. 


Ne (A7) 

For /=0, we write iNomikg Reiley Sone « Tus Ue Baxi (A8) 

in which Xg=%o+7RU[In(7/R)—1/2]@). __...... (A9) 
Blatt and Jackson’s expansion is now written 

N=X + PXO+ RX + 2.2. + RBM g(a)[S(O+ PSO+ ...] 0.00. (A 10) 


in which the functions are already specified. Observing that the first few terms 
of the expansion in the energy of h, («) are 


Mg i= (4 )ROR? AA Z/IS)RER Se oe eye ce suse see (A 11) 
we have 
X= N., 
xX, ON Ve CR O)Sgee eo ee tue eh ote (A 12) 


X,) = N, — (R/3)S,© — (2R3/15)S,5. 
Similarly, for /=1, we write 
Ny =X, +(QRUR+IU4R)hg(a)S, ne wee (A 13) 


* K, is here the function defined by Watson, Theory of Bessel Functions (Cambridge, 1922), 
p. 79. For v integral, it differs from the function used by Breit by the factor (—1)”. 
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2 which 
= (1/37)b, + 7?(9R)-l{R? + 1/4.R?)(In(r/R)—11/6)—R7}0,. is. se (A 14) 
ae coefficients X,, of the expansion 
Ny = XV + RX + RAXY + . + (ORME + 1/4R2)ho(x){SoP + R28 + . . .} 
ee (A 15) 
are then XV=NY, Xy0 =N,® —(1/108R)S,, | . (A 16) 
XO = NO — (1/108R)LS,0 + (22R2/5) SP]. J 


APPEND ITEXSE 


DIFFERENTIAL SCATTERING CROSS SECTION FOR THE 
*P STATE AND THE *F, -COUPLING 
The following is the contribution to the triplet differential cross section for 
scattering due to the states °P), 3P, and °P, along with the coupling of the ?F, 
state to the °P, state. ‘The Coulomb scattering is taken into account completely. 
We write the cross section in the form 


(FdS_/dQ) = (PdS4/dQ), + 2dS5/dQ)xy + CdSo/dQ)ingg ss see (B 1) 
in which the subscripts C, N and int denote the purely Coulomb and nuclear 
scattering and the interference of the two, respectively. 


3 if 1 2 cos (« log tan? 6/2) 
CaS dQo= gra bs ah (saa 6/2 * cos*Oj2 sin®0/2cos*o2 )f (Bd) 


dS_/dQ)y = = Al / 2p, + (4/3) sin? 8) +4 sin? 5, + (20/3)d, sin? 8, 
20/3) int (0) sing) 4 nt — 0-8 
~ 6, sin? (8) — of) —8,) — (2/3) f sin? (f yf ~8,) — 4a, sin?® (8-8) (9) 
+ 2[ — pz —(1/9)f sin? (y® —y™ — 89) + d; sin? (32 — 8p) 
+ fe sin? (6, —6,) +d, sin? (y®) — y®) — 85) + 3cy sin? (y®) —y® —8,)] Y,1Y,1 * 


6(2 : : 
T9721 “al = peered — 206, sin® (yy —y)) —20b, sin? 3, 
1 


— 2c, sin? (y1) — 8) + 85) + 9a, sin? (y) — y®) +81) —9ay sin? (8, — 85) 
— 2c, sin? (By —8:) +4) V/6g sin’ a =/91+8,) | V,PY," 


1 4(21/6+w,+ a) her pean 
aia 4/21 | - pepe ay =; Co sin (A ay +o) 


ae €y sin? (5) —3,) 


“vie sin? (6, 6,)+ 37 sin® (70 —/48,) | YAY,! * 
4(16—18w,? — 
+55 5 | eee) yn 


2+w +o 
2 190\/6w — 1821/6a— 168 — 160,” 
4 Dt wnt op pa os 
a 6a, —71/6w2— 10w,?—8 : 
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In the above we have abbreviated as follows: 


1 2b, 2b, 
b= fgg Oh (5 aD) 


= O24 3a), aie ae, AS yh 


a = 9/6 + 7w, — w, — 34/6,” 
f= b,b.[/6(w,+@2)—1], g= Sag a ae 
_ 2r/6w," — 30, —w, — 61/6 
eo care ark 
and have for convenience written w, =w»%, wy = Ww", 5, =8%. 
0 1 5 S 
(dS6/4Q) in, = Y4° ye (127)"? {1 + 3 tot 3 A bos . 
> b, sin {a log cos? 6/2 +2(y9—y)} — sin {x log sin? 6/2 + 2(yO — y)} 
cos* e By sin? 0/2 


0 Nes 04 SNC RSE SER Wi he) 
ays 2, (28m) b, {-35 5 ee 


byte a 


sin (« log sin? 6/2 + 2(y— y®) + 8,)) 
z sin? 6/2 
in which, with 7=0, 1, 2, 


jpn eset Os 207) 
. cos? 6/2 
cos {a log sin? 6/2 + 6; + 2(y) — y)} 
a sin? 6/2 
sin {a log cos? 6/2 +2(y@—y)} — sin {« log sin? 6/2 + 2(y%) — y)} 
garaaieeiecoss0,2 50 HONG We 8zte ler Neintd2o ml ae Le 


In the derivation of the results (B3), (B4), use is made of the fact that the 
3P, state is uncoupled, 1.e. wy*=1, wyY>-— 0. The phase shift pertaining to this 
state was labelled 5). 

It is readily pute that in the limit of no coupling of the j7=2 states, 


i.e. w,=1, we>— 0, the form of the cross section formula reduces to the usual 
formula oe central teres 
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Average Numbers for the Development of the Nucleon 
Component of the Cosmic Radiation in Air 
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ABSTRACT. Comprehensive results are presented in graphic form for nucleons with 
energies equal to or greater than 500 mev., based on calculations using a new cross section 
for nucleon-nucleon collisions and in the light of recent experimental evidence. The 
results obtained are in good agreement with existing experimental observations and it is felt 
that the high-energy nucleon component of the cosmic radiation can now be satisfactorily 
accounted for by our theory. ‘The basic assumptions underlying the theory are discussed. 


§1. INTRODUCTION 

N a number of recent works (Messel and Ritson 1950, Janossy and Messel 

1951, Messel 1951 a,b) an attempt was made to account theoretically for 

the experimentally observed behaviour of the nucleon component of the 
cosmic radiation in air. ‘The theory of a nucleon cascade was given for both 
homogeneous nuclear matter (considering the nucleus to be infinite in size) 
and inhomogeneous matter (a finite absorber) with and without taking ionization 
losses into account and for various initial conditions. It was shown that the 
above theory was capable of accounting at least qualitatively for all the salient 
features of the high energy nucleon component of the cosmic radiation. Since 
the appearance of more accurate data on the primary spectrum of the cosmic 
radiation it has been possible to adopt a more accurate cross section for 
high-energy nucleon-nucleon collisions and thus to obtain a theory for the 
nucleon cascade which is both in quantitative and qualitative agreement with 
experiment. We shall present results below on the nucleon component which 
are in very good agreement with experiment. 

Though quantitative agreement is obtained, we do not in any way claim that 
the cross section used is the correct cross section for high energy nucleon-nucleon 
collisions. At the same time it appears unlikely that, when the true cross section 
does become known, its energy dependence would differ radically from the one 
used by us. 

Since there has been a good deal of misapprehension regarding the 
fundamental assumptions underlying our theory we shall endeavour to give 
these in detail and also to show that the results are valid on either the plural or- 
multiple meson production hypotheses. 


§2. ELASTIC AND INELASTIC SCATTERING 
Below primary nucleon energies of the order of 1,000 Mev. it appears that 
elastic scattering of nucleons is the dominant process, with the production of 
mesons due to inelastic scattering of the nucleons playing a minor role. For 
energies greater than 1,000 Mev. meson production due to inelastic scattering 
becomes important. In all the work which follows we shall only consider primary 
nucleons with energies equal to or greater than 1,000 Mev. and assume that the 
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scattering for this energy region is inelastic. Thus in a collision of a nucleon of 
energy Z, with another nucleon at rest we have Ey = E, + E, + E,, where E, and E, 
are the energies of the secondary and recoil nucleons and £,, is the energy loss 
of the primary into a meson or mesons. 

Now, following Heitler and Janossy (1949), we assume that the cross section 
for a nucleon-nucleon collision leading to secondary energies F,, E,+dE, and 
E,, E,+dE, is a homogeneous function of degree zero of the primary and 
secondary energies. ‘This assumption is probably equally valid on either scheme 
of meson production,. plural or multiple, for the energy region we are considering. 
Thus we have for the probability that in a nucleon-nucleon collision a primary 
nucleon of energy E, should give rise to secondary and recoil nucleons of energies 
and Es, 


WL fy, H,)\dE dE, =w(E,/Ey, E,/E,)E, “dE, dE, EF, +E, <2). ...... (1) 
From this it follows that the total probability for a collision is given by 
L-l=HiE, /E FF 
| [ w(t, —*)E,%dE,dE,=constant, ...... (1a) 
0/0 Ey £y 


which is independent of the primary and secondary energies. The function w 
is so chosen that the total probability is normalized to unity. 

We require no further assumptions in order to obtain the solution of the 
fundamental equations of a nucleon cascade. The actual form of the cross section 
appearing in (1) is discussed in a later section. 


§3. THE FUNDAMENTAL EQUATIONS 

Denote by S“” (E, 6)dE the average number of protons or neutrons (according 
to the suffix k) in the energy range F, E + dE at a depth Ad gm/cm? in inhomogeneous 
matter, due to a primary nucleon of energy Ey. A is the interaction mean free 
path of nucleons with nuclei in the material concerned. The suffix 7=1, 2 refers 
to the original primary particle and k=1, 2 to the secondary particles. The 
number 1 refers to protons and 2 to neutrons. The quantities N@”(E£,, E, 6) 
denote the corresponding average number of particles with energy greater 
than E. We denote the ionization loss of nucleons per interaction length by f. 
Also define V(E/E’)dE/E’ as the probability that a nucleon of energy E£’ falling 
on a nucleus gives rise to one or more nucleons such that one of the emerging 
nucleons has an energy in the interval FE, E+dE. ‘The fundamental equations 
are now given by 


z S% (BE, 0) +S (E, #) = : { {St (EB, 0) + S38 E’, 0)}V(E/E’)d EB’ /E' 
0 


+3, 48 a Sea O Ven Mar ehh e ekise, (2) 


The equations (2) and their exact solution were given by Messel (1951 c). Since 

in the work which follows we shall not be considering secondary energies less 

than 500Mev. we may to a good approximation drop the ionization term 

(Messel 1951 c). 

We take for the initial condition a proton integral power law spectrum 
ti 

NOE, E,8=0)= pee EE. 

1] E<E 


c 


(E, is the cut-off energy). 
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The average number of particles with energies greater than E is then given by 


aap Sue oay a yds 
(1, %) ie et (1, &) Nees 4 
NG. *(E., E, 6) mal. Ce) Pye (5,0) ES, O<%<y (4) 
The integration is carried out along a line parallel to the imaginary axis such 
that the condition 0<s)<y is fulfilled. In (4) 


Poh \(s, 6) =3{exp [—f(Dge,)0) 4(—1)* hexp (—0)} "esas (5) 
(oy Sieee eres Cae ae) eae (Gay 
ae l | (URE ual ed ee (55) 


We note that «,.. expresses the average energy loss into a meson or mesons. 
D, is the average number of collisions which a nucleon suffers on making a 
diametrical passage through a nucleus of atomic weight A. Its value will be 
discussed in the next section. 

For the number of nucleons (neutrons plus protons) with energies greater 
than £ at a depth A6 gm/cm?”, we have 


ea 7 oe ds 
NE, Ee) = =| “’. (=) exp {—f(D4,)6} cay OF Si eee (6) 


For E>E, the solution of (6) is obtained meres by noting that the integral 
has a pole at s=y, hence 


N( E28) =(E/EY expt —f(D72,)0),) She Vie eee (7) 
For E<E, (6) may be evaluated using a saddle-point integration (as may be (4)). 


$4. CHOICE OF CROSS SECTION 


The solution of equations (2) was given utilizing so far only the one 
assumption given by (1). In order to obtain numerical results from the 
formulae (4) and (6) we must choose explicit values for the probability function (1) 
and for D,. ‘Taking the nuclei to be spheres with radius 1:37 x 10-18 4% cm. 
we find D, =1-5A4?(R,/1-37 x 10-#)? where R;, is the range of nuclear forces 
and A is the atomic weight. We take for the range of nuclear forces R;,=h/z,,¢ 
with z,, equal to 286 electron masses. Thus for air D, =3:-7. 

We now proceed to choose our cross section utilizing the following three 
results, deduced from experiments : (a) the interaction mean free path A for 
high-energy nucleons in air is of the order of 75 gm/cm” (Sitte 1951); (6) the 
absorption length for high-energy nucleons in air is approximately 125 gm/cm?: 
this result has been confirmed many times (McMahon et al. 1950); (c¢) the 
integral proton primary power law exponent y appearing in (3) has the value 
y =0-90 to 1-1. In our previous work we used the value 1:7 which at the time 
the work was being done appeared to be the accepted value. The results of 
Van Allen and Singer (1950) and Winckler et al. (1950) have shown that the 
value 1:7 is too high and that a value close to 1-1 is valid for energies up to 
15,000 Mev. ‘There has also been mentioned the possibility that y increases for 
energies greater than 15,000 Mev. ‘To date there has been little or no evidence to 
confirm this ;* however, in order to take into account this eventuality throughout 
we shall adopt throughout our work the upper limit y=1-1. Our results depend 


m 


* See M. Schein and J. J. Lord, Phys. Rev., abstract of meeting of 26th April 1951. Evidence 
is given for y not increasing very rapidly for very high energies. 
g g g 
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fairly critically upon the choice of y, consequently it is not surprising that our 
former work, where the value y=1-7 was used, gave only qualitative agreement 
with experiment. 

Using (7) and (a), (6) and (c) above, we must have the following condition 
fulfilled : (3-7 «.,)/75 =1/125. From a table of the function f we find 


ie hg gage ORS es eee (8) 


Recall that «,,, expresses the average energy loss into a meson or mesons. 
In other words, if we were to choose y = 1:0 the average energy loss per collision 
into the meson component would be immediately defined. 

One may now choose a general form for the function w(e,, <9), (where ¢;= E;/Ep) 
and adjust it so that the condition (8) is satisfied. We omit the straightforward 
but laborious detail here and only give the results. It should be pointed out that 
if one wishes to obtain a reasonable function w/(e,, €,) then the choice is not as broad 
as one would expect, especially when one takes into consideration the fact pointed 
out in the previous paragraph. We find that 


W(€qr€q) =120ee,(1—ey—e€g)  =«=§-— = nneene (9) 


satisfies the conditions above and also yields «.)=0-33. This cross section is 
to be preferred on physical grounds to the type used previously in our work. 
We note that, for e, =0 or ¢, =0 or €, +€, =1, w(e,, €2) now vanishes, thus ensuring 
that our collisions are inelastic. The value of 33% is sufficiently high to be 
compatible with either a plural or a multiple meson production hypothesis. 
However, if one accepts the hypothesis of plural meson production then it is 
possible immediately to work out the meson spectrum for the energy region we 
are considering. Since quantitative agreement with experiment for the nucleon 
component has been obtained using the cross section following from (9), one 
would expect the same to hold true for the meson component—providing the 
plural hypothesis is correct. If the results obtained were incompatible with 
those observed experimentally, then a good deal of doubt would be thrown 
on the plural hypothesis. On the other hand we cannot hope to disprove a 
multiple meson producing hypothesis by the above method. One may always 
argue that the multiplicities to be expected in the energy region we are considering 
are too low to be critical. We hope to present results shortly along the lines 
mentioned above on this highly controversial point which at present we feel is far 
from settled. 


§5. RESULTS 
In Figures 1,2 and 3 we give the results for the vertical intensity of protons, 
neutrons and nucleons respectively above a given energy £ for various ratios E/E. 
Our results should not be trusted for F less than 500 Mev. We hope to present 
results shortly for the low energy nucleon component, making use of the high 
energy results. In Figure 4 we have plotted the ratio of neutrons to protons. 


(i) Latitude Effect 
The latitude effect at any given depth and for any given energy EF equal to or, 
greater than 500 Mev. may easily be obtained from our curves. One merely 
measures the difference between the two curves in question, subtracts it from the 
logarithm of the normalization factor for the latitudes concerned and takes the’ 
antilogarithm. 
PROC. PHYS. SOC. LXIV, 8—A 47 
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For instance, we wish to know the latitude effect for protons with E equal 
to or greater than 500 Mev. at an atmospheric depth of 312gm/cm? between 
latitudes 0° and 55° N. where F,=15,000 and 1,000 Mev. respectively: from 
Figure 1 we measure the difference between the curves for F,/E=30 and 


05 
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Figure 1. Log N G1) (Ec, E, 8), logarithm of the intensity of vertical protons with energy greater 
than E due to the primary proton spectrum given by (3), plotted against the atmospheric 
depth @ for various values of the ratio Ec/E. 
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Figure 2 Log N qa, 2s FE, 8), logarithm of the intensity of vertical neutrons with energy greater 
than E due to the primary proton spectrum given by (3), plotted against the atmospheric depth 
@ for various values of the ratio Ec/E. 


2 at 312gm/cm?. We find this value to be 0-84. This has now to be subtracted 
from the normalization factor, log (15 x 109/10°)!1, and the antilogarithm taken; 
the result is 2-9. Conversi (1950), working on protons under corresponding 
conditions (actually he considers a differential energy spectrum) finds a latitude 
effect of 3-2+0-5. McMahon et al. (1950), using an ionization chamber, find for 
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protons with a minimum energy of about 400 Mev. a latitude effect of 2-0 between 
latitudes 20° and 55° N. If we apply an inverse Gross transform to our results. 
we also find a value of 2-0. 


; 
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Figure 3. Log N(Ec, E, 8), logarithm of the intensity of vertical nucleons with energy greater 
than E due to the primary proton spectrum given by (3), plotted against the atmospheric 
depth 6 for various values of the ratio Ec/E. 
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Figure 4. A plot of the ratio of neutrons to protons against the atmospheric depth @ for various: 
ratios Ec/E. 


(11) Position of Maximum 
From (4) and (6) it is readily seen that the maximum value for protons, 
neutrons and nucleons is attained when f(D4«,)=0. This occurs for s~0-55. 
Thus, for example, in the case of nucleons @,,,, =54 log (E/E) gm/cm’, E<E,. 
This expression shows that for a fixed E, as one passes from the geomagnetic 
47-2 
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equator (large E,) to northern latitudes (smaller £,), the position of the cascade 
maximum shifts from a high to a small value of 6. From Figure 3 it is also 
evident that the height of the maximum decreases simultaneously. The same 
behaviour is naturally also shown by the proton and neutron components, 
separately. The above points have been verified experimentally both by 
Whyte (1951) and by Staker (1950) and our agreement is good. 


(iii) Absorption Length and Neutron—Proton Ratio 


A natural consequence of the simultaneous shift of the position and height 
of the cascade maximum for a fixed E, discussed in (ii), is the increase of 
absorption length for the vertical intensity, as one passes from northern latitudes 
to the geomagnetic equator. Thus, suppose we are interested in nucleons with 
E equal to or greater than 1,000Mev. The absorption length (Figure 3) at an 
atmospheric depth of 300 gm/cm? increases from a value of 126 gm/cm? at 57°N. 
to a value of 170 gm/cm? at 0°. The absorption length for neutrons is larger than 
that for protons. We wish to point out that if an inverse Gross transform is 
applied to our result then at northern latitudes the above trend for protons is 
actually reversed. Significant experimental results on the above points are still 
lacking—especially for the energy region we are considering—although the trend 
appears to be that given above. 

For the neutron-proton ratio our agreement with existing experimental 
evidence is good. For instance, at an atmospheric depth of 650gm/cm? the 
Bristol group * finds the value N/P=1-0 for #>1,500 Mev. and N/P=1:0+0-1 
for E>3,000 Mev. Our results give corresponding values of 1-03 and 1-1. At an 
atmospheric depth of 45gm/cm? the Bristol group finds N/P=~0-6+0-1 for 
E>1,500mev. and N/P=0-2 for E>3,000Mev. Our corresponding values 
are 0-45 and 0-15. In all cases we note the extremely rapid build-up of the 
neutron component. 
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ABSTRACT. Bursts in an ion chamber, accompanied by the discharge of one of more of 
18 surrounding Geiger counters separated from it by 15 cm. Pb, have been studied at the 
Jungfraujoch, at sea level and below ground at a depth of 60 m. water. With 15cm. Pb over 
the chamber the bursts at altitude were separated into those associated with air showers 
(accompanied bursts), and the remainder (local bursts). The accompanied penetrating 
bursts appear to be associated with the incidence of high density regions of air showers. 
The local penetrating burst size increases with the number of penetrating particles recorded 
consistent with their interpretation in terms of mixed showers involving the simultaneous 
production of charged and neutral mesons in energetic nuclear collisions of the N-component. 

Underground, about 10% of the bursts under 15 cm. Pb are of the mixed shower type. 
These underground bursts are attributed to nuclear excitations produced by fast ~-mesons, 
resulting again in charged and neutral meson emission. About one-third of the bursts 
observed below ground with no Pb over the chamber appear to be produced by single 
penetrating particles which release, by an unidentified process, 20 Mev. and more in the 
ion chamber. 


Si INTRODUCLION 

HE work described in this paper started in 1946 in this way. It was 

known that most of the particles in the atmospheric showers were 

electrons, and that these showers seemed to be examples of electron—photon 
cascades. On the other hand, Rogozinski (1944) had shown that they were not 
pure cascades since they contained a small admixture of particles more 
penetrating than electrons, presumed to be mesons. ‘The question of the origin 
of the penetrating particles is of some interest and one that has received a good 
deal of attention, particularly from Janossy, Greisen, and their collaborators. 
It occurred to us that the same question could be investigated in another way 
by selecting, with an ion chamber, large cascades produced locally in an absorber, 
and at the same time looking for penetrating particles associated with these 
ionization bursts. For the latter purpose we could use, below the ion chamber, 
an arrangement of Geiger counters sufficiently shielded to prevent their discharge 
by the electrons from the burst. A coincidence between an ionization burst 
and the discharge of the shielded counters we will call a penetrating burst. 

A preliminary investigation with a pilot set showed the feasibility of the 
experiment (George and Janossy 1946) and a more elaborate scheme was 
prepared. The first results obtained with the latter arrangement were described 
by George and Trent (1948). Later, an experiment along similar lines was 
proposed by Greisen (1948). 

A parallel investigation of penetrating bursts was also carried out by 
Rossi and his collaborators (Bridge and Rossi 1949) with a different 
experimental emphasis. Their results were in agreement with ours in suggesting 
that at high altitude the majority of the bursts recorded under thick absorbers 
are due to nuclear interactions of the N-component of the cosmic radiation. 
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Our observations fit into the picture that has now come to be accepted of 
high-energy nuclear collisions giving rise to the emission of charged and neutral 
mesons. On this view, the latter quickly decay into photons, giving rise to a 
soft cascade shower which registers, in our experiments, in the ion chamber. 
The former discharge the counters set to detect the associated penetrating 
particles. The use of an ion chamber has the advantage that it is. possible to 
investigate showers that are too large for quantitative study with a cloud chamber. 


§2. SCOPE OF THE PRESENT INVESTIGATION 

Observations were made first at the Scientific Station, Jungfraujoch, at an 
altitude of 3,457m. Bursts under 15cm. Pb only were recorded. ‘These were 
found to be frequently associated with groups of penetrating particles, and 
size-frequency distributions of these groups were determined separately for 
bursts associated and not associated with air showers. ‘The multiplicity of 
penetrating particles was found to increase with burst size for the unaccompanied 
bursts. 

A comparison of the frequency of bursts at the Jungfraujoch with that at 
sea level gives the fraction of bursts produced by the hard and the N-component 
at these two places. 

The observations were repeated in our underground laboratory in London, 
at a depth equivalent to 60 metres of water. Bursts observed here must be 
attributed to the penetrating component, i.e. the .-mesons. We found that 
the proportion of bursts associated with groups of penetrating particles is much 
reduced, but that this type of burst still occurs at these large depths. On 
analysis, the results indicate that probably about 90°% of the cascade bursts at 
this depth arise from bremsstrahlung and knock-on electrons produced by the 
energetic yz-mesons, while the remaining 10°, appear to arise from nuclear 
interactions of fast u-mesons resulting in charged and neutral 7-meson emission. 
The frequency of this latter type of burst is consistent with other recent results 
on the nuclear interactions of fast u-mesons (George and Evans 1950). 

Further observations enabled us to separate out experimentally those bursts 
produced by stars occurring in the gas or the walls of the ion chamber. 


$35 APPARATUS 
The main part of the apparatus is shown in Figure 1. It consists of a 
1$-litre duralumin ion chamber I filled with 99-9°% pure argon at 40 atmospheres, 


Figure 1. General arrangement. 


surrounded by lead screens S and T and several banks of Geiger counters A....D. 
‘he ends were closed by lead screens 12cm. thick. ‘The ion chamber was kindly 
lent to us by Dr. H. Carmichael, and has been described by Chou (1948). 
Negative ions were collected, the collection time being approximately 0-1 sec. 
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The burst pulses were amplified. and recorded by means of a cathode-ray 
oscilloscope and a continuously moving film camera. ‘The smallest size of 
burst measured corresponded to the passage of 20 minimum ionization particles 
through the chamber. 

Counter banks A and D were alternative arrangements not employed 
simultaneously. A hodoscope was used to record coincidences between the 
ionization bursts and the discharge of the Geiger counters. ‘The hodoscope 
possesses no novel features and need not be described. ‘The frequencies of 
accidental coincidences were determined experimentally and shown to be 
negligibly small. 
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§4. OBSERVATIONS AT THE JUNGFRAUJOCH 

At the Jungfraujoch the lead screen T was permanently in position, and 
its thickness was not varied. Counter bank A was not used, while counter bank D 
was. The requirement for operating the hodoscope was an ionization burst, 
plus a threefold coincidence in B plus one or more counters discharged in C, 
which may be written (I, B;, C.,). On the occurrence of such an event, the 
hodoscope indicated which other counters in C or counter groups in D were 
also discharged. 

The counter banks D were designed to indicate those events associated with 
air showers, and bursts in coincidence with counters D (I, D) will be referred 
to as ‘accompanied bursts’, while anticoincidences (I, —D) will be called 
“local bursts’. 

The apparatus was operated at the Jungfraujoch during September and 
October 1947. ‘The integral size—frequency distribution of the different types of 
burst that we classified are given in Figure 2. More detailed results are given 
separately for the local and accompanied bursts. 


Penetrating Bursts 


Accompanied 
Penetrating 
Bursts 


Frequency of Bursts of Size >ze (hr-') 
Frequency of Discharge of 7 Counters 
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Figure 2. Size-frequency dis- Figure 3(a). Localpenetrating Figure 3 (b). Accompanied 
tributions under 15 cm. Pb bursts. penetrating bursts. 
at the Jungfraujoch. 


4.1. Local Bursts 


From Figure 2 it is seen that a good proportion of the bursts are in fact 
penetrating bursts. A differential size-frequency distribution of the discharge 
of groups of counters in bank C is shown for the local bursts in Figure 3 (a). 
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The frequency of discharge of m screened counters in coincidence with local 
bursts F(z) shows a rapid decline with increasing m, and may be approximately 
represented by F(n) an. 

Another feature of interest concerning the local bursts is the fact that there 
is a marked correlation between n, the number of counters discharged in bank C, 
and the size of the associated ionization burst. ‘This is illustrated in Table 1, 


Table 1 
No. of counters n discharged in C 2 3 4 5 6 Ff 8,9 
Mean size of bursts in coincidence with 
nm counters (rays) 23°31 41 51 103s 06 aeaon 


from which the steady increase of burst size with n will be seen. 

One notable feature of the events associated with the discharge of more than 
one counter in bank C was the frequent occurrence of wide separations between 
the discharged counters. Details are given in Table 2. This suggested that in 
these events we were not dealing with knock-on electrons in bank C. 


Table 2 
No. of counters n Events with g intervening undischarged counters 
discharged in bank C q=0 a G2 q=—3 q=4 q25 

2 16 10 8 5 6 115 

3 4 6 8 8 6 11 

A 5 0 3 4 2 8 

5,6 1 1 5 a 1 5 

>6 0 3 4 3 3 2 


4.2. Accompanied Bursts 


The situation is different for the accompanied bursts. ‘The size-frequency 
distribution of the C counter bank discharges is shown in Figure 3(d). The 
distribution for these bursts no longer shows a rapid decrease with n, but has 
a flat maximum in the vicinity of m about 5 or 6. The results are consistent 
with a Poisson distribution with n=6. Furthermore, there is no marked 
correlation of burst size with m for the accompanied bursts, the mean burst size 
in coincidence with 7 counters in bank C is distributed over m at random. 


§5. OBSERVATIONS AT SEA LEVEL 

Our observations at sea level have been somewhat limited. It had 
been intended to take a series of observations at sea level after a run of 
about twelve months below ground, but interesting and unexpected results 
were found in the underground observations, and it was considered more 
profitable to continue these observations for longer than had been planned. 
A brief run of the full set was obtained before leaving for Switzerland, but the 
hodoscope results were too few to be significant. The results for bursts alone 
are of significance, and those which will be useful for later discussion are: the 
burst frequency with no lead absorber over the ion chamber, 2-7 +0-25 hr-}, 
and that under 15cm. Pb, 4-6 + 0-3 hr-1. 


§6. OBSERVATIONS BELOW GROUND 
The apparatus was moved to the laboratory on Holborn underground 
station, at a depth equivalent to 60 metres of water, and has been operated there 
since the summer of 1948. 
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6.1. T=15 cm. Pb 


A preliminary test with T=15cm.Pb showed that few, if any, of the 
underground bursts were accompanied bursts. In a run lasting 14 days, 
208 bursts were observed, none of which was in coincidence with counter 
discharges in bank D. Frequent checks of the satisfactory operation of the 
apparatus were made. ‘The use of counter bank D was therefore discontinued, 
and the hodoscope channels used to indicate the operation of the counters in 
bank A. The next observations were made with T’=15cm. lead with the same 
requirement for the production of a hodoscope master-pulse as in the high 
altitude measurements, viz. (I, B,, C.,). In arun lasting 1,550 hours, 846 bursts 
were observed (0-54 + 0-02 hr), of which 382 were penetrating bursts. Of these 
382, 336 occurred with the discharge of only one of the lower counters C. Of the 
remaining 46 events, 33 discharged 2 counters in bank C, 9 discharged 3, and 
4 discharged 4. 

The penetrating bursts underground show an increase of size with the 
number of counters discharged simultaneously in C, as already noted for the 
Jungfraujoch observations. ‘Table 3 gives the details. 


Table 3 

No. of counters n discharged in C 1 2 3 4 

Mean size of bursts in coincidence with n counters (rays) 40 42 52 66 

The distribution of the intervals between the discharged counters is given 
in Table 4. The frequent wide separation of the discharged counters was taken 
to indicate that in these events the C counters were not operated frequently by 
knock-on electrons. It was felt important to confirm this for the underground 
penetrating bursts, and this was achieved in the following way: about half-way 
through the run we placed vertical lead slats, 3mm. thick, between each of the 
counters C. These slats ran the entire length of the counter slot, and touched 
the top and bottom walls. This modification made no detectable difference, 
beyond statistical errors, to the penetrating burst type of event, although 
additional counter experiments, to be reported separately, showed a marked 
reduction in the frequency of knock-on showers. 


Table 4 
No. of counters n discharged in bank C 2 3 4 
Events in which all counters were adjacent 22 2 il 
Events with 1 intervening undischarged counter 8 4 3 
Events with 2 or more intervening 3 3 3 
Counters discharged in last row als ak! 14, 16, 18 sate 
7,49 OY 7,12, 14, 18 
3,18 3, BI S16) dls 
(The C counters were numbered serially 1... 18.) 


The integral size-frequency distribution of the underground bursts is shown 
in Figure 4, plotted logarithmically, and may be represented by an expression 
of the form F(N)=AN~S with S=1-6+0-15 for all bursts, and 1-75 +0-15 for 
penetrating bursts. The former result is consistent with the observations of 
Chou (1948) at the same station. 


6.2.7 =45 cm. Pb 


Also indicated on Figure 4 are the results for a run of 273 hours with 
T=45cm. Pb. Increasing the lead screen from 15 to 45 cm. was found to reduce 
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all the readings by about 5%. The proportion of penetrating bursts associated 
with one or more than one discharged counter in bank C did not differ 
significantly from the proportion observed with T=15cm. Pb. 

The counter bank A was only employed for a fraction of the two runs with 
T=15 and 45cm. Pb. For such time as it was employed, the A bank record of 
counter discharges always indicated the high degree of multiplicity expected 
from the large cascade bursts required to operate the hodoscope. 
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Figure 4. Integral size-frequency dis- Figure 5. 
tribution of underground bursts. 


6.3. T=Ocm. Pb 


It is well known that most of the smaller bursts observed in an unscreened 
ion chamber at sea level are caused by nuclear fragments-from ‘stars’ produced 
in the gas, walls or surrounding material of the ion chamber (Euler 1940, 
Carmichael 1948). ‘The hodoscope used above was intentionally made insensitive 
to this type of event. For the run with T=0 it was decided to modify the 
requirement for hodoscope master-pulses, making the system less selective for 
showers. ‘The master-pulse requirement decided on was an ion chamber pulse 
in coincidence with at least one counter in A, B and C, which may be written 
(I Ase Der, Cs 7): 

In a run of 580 hours, 89 bursts were observed, the exponent in the size 
distribution being again 1-6. Nineteen bursts were associated with operations 
of the hodoscope, and four of these events showed multiple counters discharged 
in banks A and B, and were classified as incident showers produced in the 
overlying earth. The breakdown of the remaining fifteen events is given in 


Table 5. 


Table 5 
No. of counters discharged in A 1 1 1 
”» ”> ” B 1 2 3 
» ”» ”» Cc 1 1 1 
No. of events i 2 1 


These events were all associated with the discharge of single counters in 
banks A and C. They could not be classified as showers, since the ion chamber 
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~walls, of thickness 1-2 cm. aluminium, were far too thin to allow the development 
of cascade showers from a single incident particle. ‘Their classification is 
-discussed further below. 


§7. DISCUSSION 
7.1. Contribution from Stars 


Before discussing the burst results in detail, we first estimate the contribution 
‘to the observed burst frequencies which may be attributed to nuclear fragments 
from stars. We base this estimate on the burst frequency, 2-7 hr1, observed 
-with the unscreened chamber at sea level. It is known (Euler 1940, Carmichael 
1948) that for the size-group in question these bursts are almost entirely due 
to star processes. From the known ratio (300:1) of the frequency of stars at 
:sea level to that at the same underground laboratory (George and Evans 1950) 
‘we get immediately the frequency of bursts due to nuclear fragments in the ion 
chamber below ground, 0-01hr-t. The required burst frequencies under Pb 
-due to stars at sea level and at the Jungfraujoch may be determined from the 
known decrease under Pb of the frequency of stars, and its increase with altitude 
(Bernardini et al. 1948, George and Jason 1949). The results are collected 
together in Table 6. 


‘Table 6 
Bursts per hour due to stars 
Jungfraujoch (15 cm. Pb) So 
Sea level: 0 Der) 
15 cm. Pb 1-6 
Underground: 0 0-01 
15 cm. Pb 0-01 


7.2. Bursts at the Jungfraujoch 
(a) Local bursts. 

The observed frequency of local bursts after subtraction of the star 
‘contribution is 14-Ohr-!. The stars are predominantly low energy events, and 
will only rarely be associated with particles capable of crossing the absorber 
S and discharging one of the counters C. In any case they would not be 
expected to give a three-fold coincidence in bank B. The observed frequency 


-of local penetrating bursts is 10-3hr-4. Since this is approximately equal to 
the frequency of all local bursts with stars subtracted, we may conclude that the 


majority of those bursts under 15cm. Pb not produced by nuclear disintegration 
fragments are produced by electron cascades which contain penetrating 
particles increasing with the size of the electron cascade (Table 1). These 
observations are consistent with the assumption that they arise from the 


‘simultaneous production of charged and neutral mesons in energetic nuclear 


collisions of the N-component of the cosmic radiation. On this view, the neutral 
mesons decay rapidly, with a lifetime less than 2 x 10-1 sec., into photons (Kaplon 
et al. 1949, Carlson et al. 1950, Lord et al. 1950), giving rise to the cascade shower 
recorded as the burst in our ion chamber. The charged mesons of sufficient 
energy produced simultaneously that escape nuclear absorption in the screen S 
are those that are recorded in the counter bank C. 

Sufficient data now exist to put these ideas on a more quantitative basis. 
Thus, assume that a nuclear excitation of energy FE is produced. From the 
data of Powell (1950, Figure 29), we obtain—by extrapolation if necessary—the 
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mean number of charged mesons n(7*) emitted, and hence, from Carlson et al- 
(1950), the mean number of neutral mesons, n(7°)=4n(7*). We also have 
available from Camerini et al. (1950), and Carlson et al. the energy spectrum 
of the charged and neutral mesons. This enables us to calculate, in a 
straightforward manner, the maximum size of the cascade shower produced by 
the neutral mesons and the number of charged mesons able to penetrate the 
lead screen S. In calculating the penetration of the charged mesons we have 
assumed they suffer nuclear absorption with a probability determined by the 
geometric nuclear cross section. For the purpose of our approximate calculation 
it is sufficient to assume that the bursts are produced by cascade showers near 
their maximum, and to neglect fluctuations. By repeating for various values. 
of E, it is possible to determine the number of charged mesons penetrating S: 
as a function of the size of burst in the ion chamber. The results, so calculated, 
are given in the full curve in Figure 5, and are seen to be in reasonable agreement 
with the observations. 

The fact that at high altitude most of the bursts under thick lead are produced. 
by the N-component has been well established by the work of Bridge and 
Rossi (1949). The remaining bursts can be understood in terms of electromagnetic. 
interactions of x-mesons resulting in knock-on electrons and bremsstrahlung. 
Taking the absorption mean free path of the N-component producing these 
mixed showers to be 120 gm/cm?2, we obtain, from a comparison of the sea ievel 
and Jungfraujoch results, the following contributions to the observed burst 
frequency under 15cm. Pb at the Jungfraujoch: stars, 51°% ; mixed showers, 40% ; 
u-mesons (knock-on electrons and bremsstrahlung), 9°%. 

Our results and conclusions concerning the penetrating bursts at altitude 
appear to be in agreement with those of Hudson (1949). 


(b) Accompanied bursts. 


Here we confine ourselves to a demonstration that a reasonable interpretation. 
of the accompanied penetrating bursts lies in the assumption that they are: 
produced by the high-density regions of air showers incident on the apparatus. 
A more detailed analysis is precluded at present by the very complexity of these: 
events, possessing, as seems probable, a high density of all the known fundamental. 
particles. 

Cocconi et al. (1949) have used an arrangement of counters comparable in 
size to Our arrangement to select high-density regions in air showers. ‘They 
found the frequency of these events at Echo Lake, 3,500m. altitude, to be 
1-4hr-!, which may be compared with our observed frequency of accompanied. 
penetrating bursts at the same altitude, 0-8 hr-t. 

Further, McCusker (1950) has shown that of the penetrating particles in an 
extensive shower, about one-quarter show the strong interaction characteristic: 
of the N-component. As a rough approximation, we assume this same ratio: 
holds in the high-density region, and that the total number of penetrating particles 
striking our apparatus is proportional to the energy of the primary particle causing 
the shower. ‘Then we would expect the frequency of events in Figure 3(b) to 
increase with decreasing number of counters 7 simply on account of the power 
spectrum of the primary radiation. We would not, however, expect this increase 
to be maintained indefinitely, since the requirement of an ionization burst 
entails a low energy cut-off. The fact that the cut-off sets in for 7 about 5 or 6, 
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indicating the incidence of between one and two nucleons per event, does not 
seem unreasonable. Thus the frequency of the events, and the number—frequency 
distribution of the associated counter discharges in the bank C, do not seem 


_ inconsistent with the assumption that they are produced by the high-density 


regions of air showers. In view of the complexity of the events, it does not seem 
profitable to attempt to push the analysis further. 


7.3. Bursts at Sea Level 


The same analysis used for the bursts at altitude gives the following 
results for the contributions to the bursts observed at sea level under 
15cm. Pb: stars, 35°4; mixed showers, 15°; u-mesons (knock-on electrons 
and bremsstrahlung), 50%. ‘These observations are consistent with the cloud- 
chamber observations of Fretter (1948) at sea level. They suggest that 
arguments for -meson spin based on burst observations under thick lead 
(Christy and Kusaka 1941) are rendered a little uncertain unless the contributions 
to the observed burst frequency from components other than the u-mesons of 
the cosmic radiation are taken into account. Actually, the main conclusion 
of Christy and Kusaka was that the value of the (u) meson spin was not unity. 
This conclusion is not modified when we allow for the fact that only about 
one-half of the bursts at sea level are due to uw-mesons. 


7.4. Bursts Underground 


(a) Bursts under thick lead, T=15 and 45 cm. 


We shall assume that all the bursts observed below ground are produced by 
the energetic p.-mesons of the hard component. The very small reduction in 
burst frequency on increasing J from 15 to 45cm. Pb is consistent with this 
assumption. 

The results in §6.1 show that approximately 50°, of the bursts are 
accompanied by the discharge of a single counter in the bank C. These events 
may be simply interpreted as those in which the fast ».-meson causing the event 
passes through one of the lower counters. The fact that the u-meson is only 
detected in about 50% of the cases may be attributed to the geometry of the 
counter arrangements. 

Of more interest is the fact that 12-5°% of the penetrating bursts discharged 
more than one counter in bank C. While we may attribute some of these events 
to the production of knock-on electrons in the screen S, the frequency and 
spatial distribution of the discharged counters (‘Table 4) shows that they may 
not all be explained in this way. Observations on the frequency of knock-on 
electrons were made with the same apparatus, the frequency observed being 2:5%%. 
Subtracting this from the 12-5°% mentioned above, we conclude that 10°% of 
the bursts observed underground are associated with penetrating showers 
containing two or more penetrating particles. ‘These results do not agree 
with those of Hudson (1949). Working in a tunnel under 2,000 gm/cm?, Hudson 
reports that none of the bursts were associated with penetrating showers. 

In Figure 5 we have also plotted the mean burst size associated with the 
discharge of m of the screened counters C for the underground observations. 
The results are seen to lie close to the calculated curve, and to the results obtained 
at the Jungfraujoch, which is consistent with the viewpoint that the penetrating 
bursts observed above and below ground are essentially of the same nature, 
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and arise from the emission of charged and neutral 7-mesons from energetic 
nuclear collisions. 

We have calculated the cross sections for the production of bursts below 
ground for the following interactions of s-mesons: knock-on electrons, 
bremsstrahlung and neutral meson production, and the results are shown in 
Figure 6. In all cases we have averaged the cross sections over the u-meson. 
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Figure 6. Cross sections for bursts produced underground by p-mesons. 


spectrum, assumed to be of the form N(£)dEadk(E+£,)-*, where Ey is the 
energy lost in ionization in the overlying material, equal to 1-2x10!ev. The 
values of the cross sections for knock-on electrons and bremsstrahlung 
production by mono-energetic u-mesons of spin 4 were taken from Christy and 
Kusaka (1941). The figures for neutral meson emission were obtained in the 
following way. We assume that the excitation is produced by the electromagnetic 
interaction of fast z-mesons in passing through atomic nuclei (cf. George and. 
Evans 1950). ‘The mechanism can be understood in terms of the Williams— 
Weizsicker photon spectrum of a fast-charged particle of energy E, 
2a dv ie 

N(v)dv= airs log p> 
where v is the photon frequency and « is e?/he=1/137. Ifo is the cross section 
for star production by photons, assumed constant in the following, then the 


cross section for the production of a star of energy in the interval W, dW is,. 
when averaged over the y-meson spectrum, 


2x dW W+E E 
o(W) dw = = —— ND ee ee: 
Cee a 7 | lo8{ W } Wee, | ) 


The value of o was fixed by normalizing to the frequencies of stars observed in 
nuclear plates exposed below ground. The resulting value of o was found to. 
be 2 x 10-*8 cm? per nucleon, a value close to that derived by Miller (1951) from 
observations of stars formed in nuclear emulsions exposed directly to the y-ray 
beam of the 340 Mey. electron synchrotron at Berkeley. At large excitations 
most of the energy is liberated in the form of charged and neutral mesons 
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(Powell 1950). The required cross section for the emission of neutral mesons of 
energy W, dW was then taken to be ®(3W)d(3 W) (equation (1)), since Carlson 
et al. report that only one-third of the energy liberated as mesons is in the form 
of neutral mesons. The results of the calculation performed in this way are 
indicated in Figure 6. 

In order to express the observed burst frequencies in the same units, we 
assumed that the experimental cross section for the production of an energetic 
electron or photon of energy E, o(£) varied in the same way as the total theoretical 
cross section (Figure 6) with & when averaged over the u-meson spectrum, but 
was not necessarily equal to it numerically. This then enabled us to perform 
the integral giving the number of bursts of size greater than N rays/cm™ sect: 


F(N)=ILA ik | VGG)PGN dE a, (2) 


where J=frequency of u-mesons cm™sec}, L=number of atoms gm™!cm?, 
Az=atomic number, x=absorber thickness in gm.cm™?; P(E, N,x) is the 
probability that an energetic photon or electron of energy EF gives a shower of 
more than N particles after a thickness x of absorber. P(E, N,x) was obtained 
from the results of Bhabha and Chakrabarty (1943), assuming the number of 
electrons to follow a Poisson distribution. Equating the integral (2) to the 
observed burst frequency enabled us to fix the numerical values of (Z). The 
observed burst frequency agrees well with the sum of the various calculated 
contributions. The calculated frequency of bursts of the mixed shower type 
arising from nuclear interactions of u-mesons is 6% of that arising from the 
production of knock-on electrons and bremsstrahlung, and this is reasonably 
consistent with the observed proportion, 10%. 

In calculating the frequency of the mixed shower bursts, we have assumed 
that the characteristics of meson emission do not depend on the mechanism of 
nuclear excitation. There is some evidence from recent observations at Berkeley 
that nuclear excitation by photons results in the emission of a higher proportion 
of neutral mesons than is found for excitation by nucleons. If the proportion of 
neutral mesons is increased to, say, 50°%, as indicated by these experiments, then 
this would improve the agreement between theory and observation by increasing 
the calculated proportion of mixed showers in the underground bursts, and would 
account for the underground points lying above the curve in Figure 5. 

The observed proportions of the various events giving rise to bursts under 
thick lead below ground are: stars, 2°% ; mixed showers, 10°{ ; knock-on electrons. 
and bremsstrahlung of »-mesons, 88%. 


(6) Bursts with lead removed, 1.e. T=0. 

The observed frequency of bursts with T=0 is 0-15hr-!, and the upper 
limit to the possible contribution from stars is 0-01 hr-!; hence we conclude that, 
contrary to the conditions at sea level, with only air above the ion chamber, the 
bursts observed underground with T=0 are not due to stars. Although there 
was no lead above the ion chamber, there was a considerable thickness of earth, 
3 m., above the apparatus, and the majority of the bursts were probably produce 
by electron cascades incident from the roof. This explanation cannot apply 
to those 15 events in coincidence with the hodoscope counters in Table 5, since 
only one counter in bank A was discharged. The frequency of these events is 
0-025 + 0:006 hr-1, and since the hodoscope efficiency is approximately 50% 
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(§7.4(a)), we may infer that the actual frequency of these events is 0-05 + 0-01 hr. 
The hodoscope records clearly indicate that most of these events are associated 
with the passage of a single particle through the banks of counters and, from 
their frequency, we know that not more than 20° of these events can be due to 
stars produced by fast yz-mesons. A few of the hodoscope records may have been 
produced by accidental coincidences, but the majority indicate that bursts of 
more than 20Mev. of ionization in the chamber were produced by single 
penetrating particles by a process so far unidentified. 


§8. CONCLUSIONS 

From our observations we have drawn the following conclusions: 

(i) 90% of the local bursts observed under 15cm. Pb in our ion chamber 
at the Jungfraujoch are produced by the N-component of the cosmic radiation, 
50% being produced by heavily ionizing fragments from nuclear disintegrations, 
and 40% from electron cascades containing penetrating particles (mixed 
showers); the remaining 10° are produced by p-mesons. 

(ii) The number of penetrating particles in a local penetrating burst under 
15cm. Pb increases with the burst size, the proportion being about 1 per 
50 electrons at the shower maximum. 

(iii) ‘The observations of penetrating bursts at the Jungfraujoch are consistent 
with the assumption that they are produced by energetic nuclear disintegrations 
resulting in the simultaneous emission of charged and neutral mesons, the latter 
decaying to form cascade showers. 

(iv) About 3% of the bursts observed under 15cm. Pb at the Jungfraujoch 
are produced by the incidence of high-density regions of energetic air showers 
on the apparatus. 

(v) Under 15cm. Pb at sea level only about 50% of the bursts are produced 
by fast u-mesons. 

(vi) At a depth equivalent to 60m. of water about 10% of the bursts under 
15cm. Pb are of the mixed shower type, and presumably arise from energetic 
nuclear disintegrations emitting charged and neutral mesons produced by fast 
p-mesons. 

(vii) The frequency of these mixed-shower events below ground is consistent 
with observations of nuclear disintegrations observed in nuclear emulsions at 
the same depth. 

(viii) The remaining bursts under lead are produced by the knock-on 
electrons and bremsstrahlung produced by fast u-mesons. 

(ix) With no lead immediately above the ion chamber, two-thirds of the 
bursts are due to cascade showers produced in the earth above the apparatus. 

(x) About one-third of the bursts observed underground with no lead above 
the apparatus are produced by single penetrating particles crossing the ion 
chamber and releasing 20 Mev. or more of ionization in the chamber. The 
mechanism of this energy release has not been identified. 
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Observations of the Latitude Effect of Penetrating Showers 


By J. W. SFURGESS 
Birkbeck College, University of London 


Communicated by E. P. George; MS. received 13th February 1951 


ABSTRACT. The latitude effect of the local penetrating showers of cosmic radiation 
produced in paraffin has been measured at sea level and found to be in the ratio of 1:8+ 0-4 
between London and Dakar, French West Africa. ‘The extensive penetrating showers 
showed no latitude effect, which is consistent with the greater energy of these showers. 
The background frequency of local showers, i.e. those recorded with no absorber above the 
counter system, was found to be insensitive to latitude and reasons for this are discussed. 


Sie INLRO DUCTION 
Roe the latitude dependence of the nucleon component of the cosmic 


radiation has attracted some attention. Simpson and Uretz (1949) and 

Yuan (1949) have studied the latitude effect for low energy neutrons, while 

observations on disintegration ‘stars’ observed in photographic plates have been 
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reported by Salant et al. (1950) and Lattimore (1950). Results for the more 
energetic type of event, the penetrating showers, observed with systems of 
lead-screened Geiger counters have been given by Walsh and Piccioni (1950) 
at altitude and by Appapillai and Mailvaganam (1950) at sea level. In addition 
a series of theoretical papers has been published by the Dublin School of Cosmic 
Physics (see Messel and Ritson (1950) for a list of references). 

In the present paper, results are given of observations of the penetrating 
showers produced in paraffin at sea level at Dakar, 21-6° N and London 54:2°N 
geomagnetic latitude. 

Paraffin was chosen in order to minimize the response of the apparatus to 
the extensive showers (cf. Broadbent and Janossy 1948) which, being events of 
considerable energy, were expected to show negligible latitude effect. 

Since it seemed probable that the latitude effect would be larger the smaller 
the energy required to discharge the counter system, it was decided to reduce to 
a minimum the lead screening while still maintaining the system insensitive to 
cascade electron showers. A thickness of 15 cm. Pb between extreme counters 
was used, and separate experiments are reported which show that this screening 
was adequate. 


§2. APPARATUS 

The apparatus is shown in Figure 1. It is similar to that used by Broadbent 
and Janossy (1947), except for the reduction in the lead thickness from 30 cm. to 
15 cm., as stated above. Sevenfold coincidences were recorded between the 
discharges of the three counter groups in the top bank, and the two groups in. 


heey cm 


E 


LIL 
soos Sane 


Figure 1. The apparatus. 


each of the middle and the bottom banks of counters. Butler et al. (1950) have 
shown that the mean energy of the showers detected by the Broadbent—JAnossy 
arrangement is 10!°ev., whereas the latitude cut-off energy at Dakar is 
1-15 x 10'ev. in the vertical direction. The mean energy of the showers 
detected in the present experiment may be taken as 5 x 10% ev., well below the 
cut-off energy. 

The unscreened tray E was used in the usual way to detect events associated 
with extensive air showers, and was connected in coincidence with the main 
shower set P. ‘The resolving time of the coincidence circuits was 10-® sec. and 
the accidental rate was negligible. The paraffin absorber was cast in blocks 
45 cm. x 45 cm. x 15 cm. 
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§3. RESULTS 
(i) Extensive Showers 
Table 1 gives the frequencies of extensive showers (coincidences (P, E)) at 
London and Dakar, for various absorber thicknesses. It is seen that there is 
no significant difference between the results at the two latitudes. 


Table 1 
Paraffin absorber thickness (cm.) 0 15 30 
Extensive penetrating { London 0:17+ 0:02 0:18+ 0:02 Om Sat.0 02 
showers per hour Dakar 0:17+0-02 0:16+ 0:04 0-17+0-03 


(11) Local Showers 


Table 2 gives the frequencies of occurrence of local showers (anticoincidences 
(P, —E)) corrected for differences in barometric pressure (Janossy and Rochester 


1944). 


Table 2 
Paraffin absorber thickness (cm.) 0 15 30 
Local penetrating showers [London 0-40+40-03 0-73 + 0:04 0:85+ 0-04 
per hour | Dakar 0:-41+ 0:04 0:51+ 0-07 0:66+ 0-05 


The two straight lines of best fit relating the shower frequency to paraffin 
thickness were found by the usual procedure of minimizing the value of ?. 
The latitude effect determined from the ratio of the slopes of these two lines is 
1-8 + 0-4. 

(iii) Local Showers produced in Pb at London 

In order to verify that the coincidences were mostly caused by penetrating 
and not by electron cascade showers, the frequency of showers under various 
Pb absorbers was determined at London. The results, corrected for barometric 
fluctuations and the effect of the geometry of the top absorber (cf. George and 
Jason 1950), are plotted in Figure 2. Fitting an inverse exponential curve to 
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Figure 2. Corrected lead transition curve. 


the results, in the manner indicated by these authors, gave the value of the collision 
length in Pb of the particles generating the showers as L,=150+70 gm/cm?, 
consistent with recent results. 

Thus, from the transition curves obtained in paraffin and in lead, it appears 
that in spite of the reduction in the amount of Pb absorber between the counters 
of the set P, the events recorded were still of the penetrating shower type, and the 
latitude effect measured is that for events generally accepted as penetrating 
showers. 

48-2 
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§4. DISCUSSION 
(i) Extensive Showers 


The lack of latitude effect for the extensive showers is not unexpected, since 
the mean energy of the extensive showers recorded must have been in the region 
of 104 ev., well beyond the geomagnetic cut-off energies. 


(ii) The Background Frequency of Showers with no Material 
above the Counter System 


Table 2 shows that the background frequency of penetrating showers, with 
no material above the top tray, is the same at both latitudes. One might expect 
these background rates to show the same ratio as the measured latitude effect. 
As possible explanations of this latitude-insensitive background, the following 
mechanisms have been considered: (a) low density extensive air showers, whose 
penetrating particles discharge the main set but miss the unscreened counter 
tray E, or (b) the discharge of the top trays by back-ejected shower particles from 
an event in the lead below, such events being of a higher energy than the normal — 
showers recorded, and hence less latitude-sensitive. 

The contribution of extensive air showers has been estimated numerically 
in the following manner. For showers of density greater than 10m the density 
spectrum and proportion of penetrating particles have been taken from Cocconi 
et al. (1949). For lower density showers, where few measurements have been 
made the density spectrum and proportion of penetrating particles have been 
calculated from the primary spectrum and the known cross section for meson 
production. It was found that below densities of 1m~? the proportion increases 
markedly beyond the 2°%% observed for denser showers, because these low density 
and therefore low energy showers are well beyond their cascade maxima and 
the associated mesons are being absorbed much less rapidly than the cascade 
electrons. However, only 1% of the observed background rate can be accounted 
for by this means. 

Franzinetti and Perkins (1951) have investigated the backward ejected 
relativistic particles from ‘stars’ in nuclear emulsions, and have shown that the 
ratio of back-ejected shower particles to the total number of shower particles is 
independent of the size of the star and is about 0-10. ‘Thus the number of back- 
ejected shower particles increases with energy. 

To discharge the top tray of the set P three particles are required, but the 
incident particle is charged in about half the events so that on the average 2:5 back 
particles are required. The average number of shower particles in such an 
event is approximately 2-5/0-10=25 which corresponds to a mean energy of 
about 2:5 x 10!ev. This may be compared with the estimate of 5 x 10% ev. for 
a shower generated in the top absorber. ‘The nucleons causing the apparatus 
to discharge will have energies greater than these values. If the energy spectrum 
of these high energy nucleons is taken as N(>E)« E~1 (cf. Barton et al. 1951), 
the ratio of the frequencies of nucleons above these two energies will be 516 or 
about 12:1. 

Table 2 shows that the frequency of showers produced in 30 cm. paraffin 
is comparable with the background rate of showers. As the collision length in 
paraffin is about 80 gm/cm? (George and Jason 1950), the chance of a collision 
in 30 cm. paraffin is about 0-25. But the efficiency of recording showers in 
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: paraffin is about one-third from geometric considerations, whereas for the back- 
ground showers produced in lead it is nearer unity. Therefore the chance of 
recording those showers produced in lead with sufficient back-ejected particles 
is probably so much greater that it compensates for the lower frequency of 
nucleons of sufficient energy (>2-5 x 10'ev.), and the shower frequency from 
30 cm. paraffin is approximately equal to the background rate. 

It seems probable therefore that the small latitude effect in the background 
rate may be understood in terms of the higher energy required to produce the 
background coincidences. 


(iii) Local Showers 


The latitude effect observed was 1-8 + 0-4, a value somewhat larger than other 
recent determinations discussed below. A nucleon of energy 5 x 10®ev. which 
will discharge the penetrating shower set must have had an energy of 7 x 10% ev. 
at the top of the atmosphere. The cut-off energy at Dakar is 1-15 x 10ev., and 
if the penetrating showers arise from primaries only, the index y of the primary 
spectrum in this energy region will be given by (11-5/7)”’=1-8+0-4, whence 
y=1:2+0-4. This figure is in good agreement with the results of Van Allen 
and Singer (1950) who obtained a value of 1-1 for the primary spectrum in the 
energy range 3 to 15 x10%ev. from rocket experiments. At a residual depth 
of 15 gm/cm?, Winckler, Stix, Dwight and Sabin (1950) obtained a value of 
0-9 in the same energy range. If secondary nucleons are also considered to 
discharge the set, the energy of the primary nucleons from which these arise 
must be greater and the index will be greater. 

Messel (private communication) has extended the original work of Heitler 
and Janossy (1949) by including the effect of nucleon recoils and ionization loss, 
thus giving, in effect, a theory of the nucleon cascade. He gives a value of 1-7 
- for the latitude effect for nucleons at sea level with energy greater than 5 x 10° ev., 
using y=1-1. 


(iv) Comparison with other Observations 


The only other determination of the latitude effect at sea level for penetrating 
showers using a counter system is that of Appapillai and Mailvaganam (1950) 
who obtained a negative result. ‘There are several possible reasons for this. 
First, Appapillai and Mailvaganam compared the frequency of showers deter- 
mined by them at Colombo with the frequency observed at Manchester by 
Broadbent and Janossy, using similar but not identical apparatus. Secondly, 
as mentioned above, the average energy required to discharge this particular 
equipment is about 101° ev., and a smaller latitude effect would be expected in 
any case. ‘Thirdly, a small difference existed in the two sets of apparatus which 
would tend to reduce the observed effect. ‘The top absorber used by Broadbent 
and Janossy was 45 cm. square, whereas that used by Appapillai and Mailvaganam 
was 50 cm. square. Thus, from the null effect reported by Appapillai and 
Mailvaganam, one may conclude that the latitude effect was probably of the 
order of (50/45)? or 1-2, a not unreasonable value for showers of 10! ev. energy. 
In the present experiment, the same apparatus was used at the two latitudes, 
though it had necessarily to be dismantled for transport between them. 

A latitude effect of 1-1 was observed at sea level and at altitudes of 25,000 
and 30,000 feet by Walsh and Piccioni (1950) between 21° and 55°N. These 
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authors estimate the mean energy of the penetrating showers discharging their 
apparatus at 2 x 10! ev., which is above the latitude cut-off energy. 

Observations of the nuclear component made with photographic plates and 
with neutron detectors have in general shown a bigger latitude effect than that 
reported here; this is reasonable in view of the smaller energy of these events. 


§5. CONCLUSIONS 


The conclusions resulting from this investigation are as follows: (1) the 
latitude effect of local penetrating showers of mean energy about 5 x 10° ev. is 
1-8+0-4 between 20°N and 54°N geomagnetic latitudes. From this the index 
of the integral energy spectrum for primary nucleons for energies between 
5 x 10% ev. and 1-15 x 10 ev. is found to be 1:2+0-4. (ii) The range in lead of 
the particles producing the showers investigated is 150+ 70 gm/cm?. (iii) The 
local showers observed with no material above the top counter tray show a 
negligible latitude effect. This implies that the mean energy of the particles 
producing the showers under this condition is higher than with material above 
the top counter tray, consistent with the assumption that in the former case the 
top tray is discharged by back-ejected particles. (iv) The extensive showers 
show no detectable latitude effect. This is consistent with the large value, of 
the order of 104° ev., of the energy of these showers. 
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LETTERS TO THE EDITOR 


Investigation of Superconductivity in Lead Compounds 


Meissner, Franz and Westerhoff (1933) have made electrical resistivity measurements 
on lead sulphide down to 1:3° k. and concluded that the very pure substance is a semi- 
conductor, whilst the occasional appearance of superconductivity in specimens showing a 
metallic characteristic must be due to small amounts of excess lead. Dunaev and 
Maslakovitz (1947) have worked with specimens containing small amounts of lead impurity 
(carrier concentration of the order 101° per cm*) and find a metallic behaviour below 
450° c. but no superconductivity down to 2:15° K. 

Quite recently Darby, Hatton and Rollin (1950) reported superconductivity in PbS, 
PbSe, and PbTe all at about 5° k., as measured by a magnetic method. The materials 
were prepared by fusion together of the constituents in vacuo; no X-ray investigation 
was made, apparently, to look for excess lead. The figure 5° K. is most interesting, being 
relatively far removed from the normal lead transition temperature of 7:3°K. Lark- 
Horovitz, in a discussion of these results pointed out that it is now generally accepted 
that the chemically very pure lead compounds behave as semiconductors. If they should 
become superconducting at low temperatures this would be the only known case of such 
behaviour and would be of great interest from the theoretical point of view. On the other 
hand, x-ray investigations of lead compounds in the Purdue laboratories have shown that 
frequently lead is found if the material is electrically n-type (excess conduction). 

We have investigated the problem using the ballistic-throw magnetic method of 
detection; this enables one to distinguish most readily between a bulk effect and an 
impurity effect. The specimen under investigation is enclosed in the bulb of a gas 
thermometer, which is surrounded by a vacuum case immersed in a bath of liquid helium. 
It was hoped to reproduce the results of Darby et al. and then to show that the 
specimens were impure, but we were unable to obtain a 5° transition and the results are 
essentially a confirmation of Meissner’s work. 

Pure lead sulphide may be prepared in the form of a fine powder by precipitation from 
solutions of, for example, lead acetate and sodium sulphide. Care must be taken in the 
method of fusing together lead and sulphur, for the reaction is strongly exothermic and some 
of the sulphur is vaporized. This difficulty may be overcome by heating the mixture in 
a graphite crucible inside an evacuated, sealed-off quartz tube. After the reaction has 
taken place the temperature of the furnace is raised to the melting-point of the sulphide 
(1,120° c.) and held there for a time before cooling again. Ingots may also be prepared 
by using the precipitated powder instead of the lead sulphur mixture. 

We have made measurements on ingots prepared by fusion and on the precipitated 
powder both finely divided and sintered,* In most cases no transition could be observed 
down to 1:3° k. and the x-ray photographs contained only PbS lines.t In a few specimens 
a very small (partial volume) effect was found at 7:3°«., the normal lead transition 
temperature, and the x-ray patterns for these specimens included a very weak extra line 
presumably due to excess lead. 

Crystals of natural lead sulphide (galena) and of the pure commercial product also 
gave negative results. A fusion sample of lead telluride became superconducting at 
7:3°., the magnitude of the effect indicating that it was confined to a small fraction of 
the total volume of the specimen, and the presence of lead was confirmed by x-ray. A 
single crystal of lead telluride (kindly supplied by Mr. R. A. Smith of Telecommunications 
Research Establishment, Malvern) showed no transition down to 1:3° k. 

These experiments, in essentially reproducing Meissner’s results—and extending them 
to lead telluride—strongly suggest that these lead compounds are not superconductors. 
We feel that the results of Darby and his co-workers may be explained by the presence of 
excess lead in their specimens (a possibility admitted by the authors) since their transition 
temperatures and critical field curves are nearly the same for the three different compounds. 
A few per cent of excess lead could easily form a superconducting network of lead veins 
that would give an apparently complete volume effect (cf. Hudson 1950) and the stated 


* We are indebted to L, M. Roth for assistance in the preparation of these specimens. _ 
+ We are indebted to I. G. Geib for all the x-ray investigations connected with these experiments. 
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procedure of rapidly cooling the melt to room temperature might result in internal strains 
sufficient to alter the transition temperature appreciably. We have, however, tried to 
obtain such a strain effect without success. 

These experiments were carried out at the suggestion of K. Lark-Horovitz. Since 
their completion, Dr. Rollin has sent us an account (to be published later) of more recent 
experiments with very pure specimens which are in agreement with the above results. 


Department of Physics, R. P. Hupson. 
Purdue University, 
Lafayette, Indiana. 
22nd May 1951. 
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Effect of Temperature on the Height of Potential Barriers and 
on the Breakdown Voltage of Contact Rectifiers 


Theoretical calculations (Kimball 1935) have indicated a widening of the energy gap 
of covalent lattices, such as diamond and silicon, with decreasing lattice parameter. This 
is confirmed by Taylor (1950) who finds the intrinsic resistivity of germanium to increase 
by 10-2 10->/atm. With the values of volume compressibility and of the coefficient of 
thermal expansion this corresponds to a change in the gap of —0°87x10-*ev/°c. This is 
confirmed by Shockley and Bardeen (1950), who obtain a value of —0:9x10-* ev/°c. 
from the change of mobility of electrons and holes in germanium with temperature. On the 
other hand, Fan (1950) finds a value of —4-0 x 10-4 ev/°c. from the long-wave limit of the 
optical absorption band; this high value may be caused by increased line width (Radkowsky 
1948). 

Apart from these effects in the semiconductor itself, it might be expected that the height 
of the potential barrier arising at the interface when in contact with a metal would also be 
affected by a change in gap width. Such changes in barrier height would affect certain 
characteristics of barrier cells very sensitively and would thus make the changes highly 
susceptible to measurement. Some such effects have recently been observed here on 
rectifiers and photocells (Billig and Ridout 1951). Measurements of the contact resistance 
of standard selenium, germanium and silicon rectifiers and of selenium photocells over a 
large range of temperatures indicate that the height of the barrier in all these cells remains 
constant above room temperature, but gradually decreases as the temperature is lowered. 
Extrapolation beyond the lowest temperature reached so far (78°K.) seems to indicate a 
complete disappearance of the barrier around absolute zero temperature. Such a drop in 
barrier height is supported by an analysis of the inverse current-voltage characteristic of 
contact rectifiers at low temperatures. As has recently been shown (Billig 1951), electrical 
breakdown in rectifiers is mainly due to the onset of thermal instability and is critically 
dependent on the height of the potential barrier as well as on the cooling. As these two 
effects operate in opposite directions, a lowering of the temperature may lead to either an 
increase or a decrease in the breakdown voltage. This is actually borne out by experience : 
selenium rectifiers, for instance, when cooled considerably below room temperature break 
down at a lower voltage, whilst in germanium point contact rectifiers the turnover voltage 
actually increases at low temperatures, at least down to liquid air temperature (Benzer 1949, 
confirmed by unpublished measurements of the present author). 

Quantitatively, the results obtained could not so far be correlated with the effects reported 
above in the semiconductor material itself: taking Taylor’s value of —0-87 x 10-4 ev/°c., 
a drop of 100°c. in the temperature should increase the gap width in the germanium by 
0-0087 ev., whilst, according to our measurements, cooling the germanium rectifier by 
100°c. depresses the height of the barrier by 0:125 v. 


Letters to the Editor 753 


According to the simple model of Schottky (1942), the height of the potential barrier— 
the ‘ diffusion voltage ’—should be equal to the difference in work functions of the contacting 
metal and semiconductor. Very little information is available on the change with 
temperature of the work function, particularly in semiconductors; it is, therefore, impossible 
at this stage to decide in what way, according to Schottky’s model, the observed change in 
barrier height should be correlated with the change in width of the energy gap. Further 
measurements, particularly a study of the effect of large hydrostatic pressures on the 
characteristics of barrier cells, should provide useful evidence for or against that simple 
model of the potential barrier. Such experiments are being planned and will be reported 
in due course. 
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K Capture in I 


The nucleus 181 lies between two stable isobars: 1°T’e and !8Xe. It has long been 
known to decay by f-emission to #?8Xe; more recently Reynolds (1950) has shown by mass 
spectrometric analysis of an old neutron-irradiated iodine sample that 5% of #81 nuclei 
decay to *Te. This alternative decay process must consist principally of K capture, 
since less than one in a thousand positrons are emitted (Barber 1947). 

We have observed x-rays emitted by 1281 by means of a Nal (T activated) crystal and 
photomultiplier. They consist of Te K-radiation and are not emitted in coincidence 
with the #8I B-ray. By comparing with the x-rays from 4In we have estimated the 
K capture branching ratio. 

Sources of 1°8I were prepared by the Szilard—Chalmers method from neutron irradiated 
ethyl iodide. Owing to the fact that the K absorption edge of Ag lies below the x-ray 
energies involved, the use of silver iodide precipitate was avoided in the extraction, and 
active iodine was prepared as cuprous iodide. 

Beta-rays were screened off with 1 cm. of plastic, and the sources were placed in front 
of a thin slice of sodium iodide crystal mounted on an E.M.I. photomultiplier tube. An 
analysis of the pulse height spectrum with a single channel kicksorter showed a peak at an 
energy 1:15 times that corresponding to the peak found with "4In. This may be compared 
with the ratio of 1:13: 1 to be expected for Te and In K x-rays. Further experiments 
with absorber foils identified these x-rays as predominantly characteristic of Te. 

To support this conclusion an experiment was devised to distinguish x-rays which were 
in coincidence with f-particles from those which were not. An !°8I source was sandwiched 
between two crystals of naphthalene—anthracene and mounted on a photomultiplier tube. 
This formed a B-counter of 47 solid angle. The 47 counter and source were then brought 
in contact with the x-ray counter. A coincidence circuit generated a pulse which could be 
used either to block or to pass counts through the kicksorter. In this way the kicksorter 
could be used to analyse the spectrum in coincidence with f-decay, or not in coincidence 
with f-decay, or simply the total spectrum due to the source. 

Curves obtained with and without f-decay anticoincidence are shown in the Figure 
to the same scale. The 14In x-ray peak is added as a calibration mark. Beyond the peak, 
anticoincidence reduces the count rate to 16°, this remainder presumably representing 
the inefficiency of the 47 counter; in view of the thickness of the source and the consequent 
B-absorption and loss of light, this is a plausible figure. The reduction within the x-ray 
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peak is due to the removal of bremsstrahlung and of small pulses produced by the 428 kev. 
y-ray line of !8I, both of these being in coincidence with b-decay. Observation of the 
X-ray spectrum in coincidence with fB-decay showed the same type of spectrum as that 
obtained from the bremsstrahlung of a =2P (pure 8) source. No significant quantity of 
characteristic x-radiation could be detected. 

With the 47 counter and anticoincidence circuit we hope to proceed with the examination 
of other nuclei; it is possible that in some cases bremsstrahlung and y-ray emission have 
hitherto obscured any X-rays arising from K capture. 
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In order to estimate the proportion of 1*8I nuclei decaying by K capture, a comparison 
of B-ray emission and of x-ray emission was made between sources of !4In and 81. A 
highly converted isomeric level of 44In decays with a 48-day half-life to the ground state 
of 44In, and this in turn decays 97°% by B-ray emission to Sn and 3% by K capture to 
14Cd. K capture and K shell conversion of the isomeric level lead to a fraction (45 + 4) x 10-? 
of K shell vacancies per f-ray emitted (Boehm and Preiswerk 1949). Since both f-ray 
and x-ray energies lie near one another in these two cases, a comparison is particularly 
easy to make. The correction for Auger effect need not be made since a similar fraction of 
K shell vacancies will be thus filled in both types of atom. A branching ratio, K capture 
to B-decay, of (6-3+0-7) x 10-? was found. 

The authors would like to thank Dr. P. F. D. Shaw for his assistance with the 
Szilard—Chalmers separation, and Lord Cherwell for putting at their disposal the facilities 
of this Laboratory. 


The Clarendon Laboratory, W. B. Mims. 


Oxford. H. HALBan. 
21st May 1951. 
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Magnetic Analysis of «-Particle Groups from (d, a) Reactions 


A heavy-particle magnetic focusing spectrometer has been used to measure the energy 
spectrum of «-particles emitted during the bombardment of light elements with the deuteron 
beam from the Liverpool 37 inch cyclotron. The spectrometer is similar to those described 
by Ringo (1940) and by Burcham and Freeman (1949), and uses a 50° sector field to give 
focusing in the horizontal plane. Particles emitted within + 14° of the normal to the field 
boundary enter the field and those with a 25 cm. radius of curvature in the field are focused at 
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an adjustable slit. The target and slit are symmetrically placed with respect to the 50° 
sector on a line through the intersection of the effective field boundaries. A proportional 
counter is used as a detector. Thus different types of particle of equal momenta may be 
distinguished and «-particle energies up to 4 Mev. may be measured. ‘The vertical separ- 
ation of the pole faces is 1 cm., giving a $° angle of collection in this plane, and a total solid 
angle for collection of about 1/2400. 

The thin targets used (~1 mm. air equivalent) were deposited on strips of gold foil 
1-2 cm. wide supported vertically with their faces making an angle of 30° with the deuteron 
beam. In the target chamber the beam was about 2 cm. wide and 3cm. high, so that a 
target area of about 3 cm? was illuminated. The analyser was arranged to measure particles 
emitted at a mean angle of 90°. A proportional counter was used to monitor the beam 
passing through the target by counting elastically scattered deuterons in a small solid angle 
at 45° to the deuteron beam. Absorbers could be placed in an air gap between the cyclotron | 
extension tube and the target chamber to reduce the incident beam energy. The energy 
resolution was limited to 10°% by the energy spread of the incident deuterons, the spread in 
angle between incident and emitted particles, and the width and thickness of the targets. 

The current to the high resistance windings of the spectrometer magnet was supplied 
from a 1,200 v. rectifier set and was stabilized electronically. "To measure the magnetic 
field, the charge generated in a flip coil was compared with that stored in a condenser charged 
to a voltage known by comparison with that of a standard cell. As the spectrometei was in 
the stray field of the cyclotron magnet it was calibrated in situ by slowing down a-particles 
from Thorium C in an air absorption cell. 

The «-particle spectra observed from thin deposits of carbon and zinc oxide are shown in 
Figures 1 and 2. The higher energy peak (2°81 Mev.) in the carbon spectrum corresponds 
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to the ground state of B, giving an incident beam energy of 6:77 Mev. using the Q value 
—1-44 mev. for this reaction. This beam energy has been used to calculate the O values 
and energy levels given in the Table. The accuracy of the excitation energies is estimated 
to be closer than + 0:08 mev. 


Reaction 220 (daca) 8 16QO(d, «)14N 16O(d, «)4N 16O(d, «)14N 
«-particle energy (Mev.) 2°31 2:28 2°83 3-70 
QO (mev.) —2-15 —2:64 —1-95 —0-89 


Excitation energy of residual 


nucleus (Mev.) 0-71 5-70 5-01 3-95 


The «-particle group corresponding to the known level in 1°B at 0:71 Mev. (Hornyak et al. 
1950) has also been observed in photographic plates by Berlman (1950). A range of energy 
levels of 9B up to 1:3 Mev. would be covered by our measurements, but they show no 
a-particle groups corresponding to levels at 0-41 and 1:02 Mev. suggested to explain y-ray 
lines from 1°B* (Hornyak et al. 1950). 

There is a possibility that the 2-28 and 2:83 Mev. groups from oxygen might be due to 
carbon contamination although measurements with freshly prepared zinc oxide targets 
indicate that this is unlikely. The «-groups corresponding to known levels of 4N at 3-95 
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and 5 Mev. (Hornyak et al. 1950) have also been reported by Berlman (1950), while the 
2:28 Mev. group would correspond to the suggested level at 5-5 Mev. (Hornyak et al. 1950). 
A range of energy levels of 14N from 3:4 to 6:1 Mev. would be covered by the spectrum shown 
and this range has been extended to 2:0 Mev. by using a reduced deuteron beam energy of 
5:31 mev. No further groups were discovered. There is no sign of a-particle groups 
corresponding to the known levels at 2:3 and 3:4 mev. or the doubtful level at 6-1 Mev. 
(Hornyak ez al. 1950). 

Using the reduced beam energy of 5-31 Mev. we have also observed a 4°32 MeV. a-particle 
group from magnesium. As this element has three isotopes, this could correspond to 
excited states of 2:90, 8:06 and 3-05 mev. for ?2Na, ?*Na, ?4Na respectively. Since 4Mg 
is considerably the most abundant isotope, ?2Na is the most probable residual nucleus. 

A more detailed account of this work will be published elsewhere. 


George Holt Physical Laboratory, A. ASHMORE. 
University of Liverpool. J. F. RAFFLE. 
3rd May 1951. 
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The Reaction “C(an)'*O 


We have established the existence of the reaction C(an)!®O by bombarding a thick 
target of separated °C with alpha-particles of 1:05 mev. The neutrons were detected 
at 90° to the alpha-particle beam by studying the pulse-size distribution due to the recoils 
they produced in an ionization chamber containing deuterium at 38 atmospheres. The 
energy scale was provided by photoprotons ejected from the deuterium by gamma-rays 
of 6:14 Mev. ‘The energy of the neutrons from the *C was found to be 2:68+0:15 Mev.; 
this determination was based on about 20,000 deuteron recoils. The thin target energy 
expected from the atomic masses is 2:92 Mev. 

The excitation function rises slowly throughout the range of alpha-particle energy 
investigated—up to 1,200 kev. This probably corresponds to the formation of 17O in 
the known broad state around 7:7 Mev. found by Nuckolls et al. (1946) in the total neutron 
cross section of ®O. 

We should like to thank Dr. W. D. Allen and Sir John Cockcroft for the supply of the 
18C from the Atomic Energy Research Establishment, Harwell. 


JONES. 


Cavendish Laboratory, G. A. 
D. H. WILKINSON. 


Cambridge. 
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Angular Distributions of the Two Long-Range Neutron Groups 
of the Reaction ‘°O(d, n)'’"F 


The angular distributions of charged particles emitted by certain light elements subject 
to deuteron bombardment have recently been published by several authors. Previous to 
this communication there have been no reports on the angular distributions of separated 
neutron groups. Using the photographic plate technique, we have investigated the neutron 
spectra and angular distributions of a number of light elements bombarded by 8 mev. 
deuterons. In this note we present the angular distributions of the two long-range neutron 
groups from the reaction O(d,n)!’F. A detailed account of the experiment will be 
published elsewhere. 
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A lead peroxide target of thickness 0:15 Mev. deposited on a thick gold backing was 
bombarded by the 8 Mev. deuteron beam of the Liverpool 37-in. cyclotron. ‘Twenty 
photographic plates, of emulsion thickness 4001, were orientated at 5° intervals about the 
target, and simultaneously irradiated. From measurements on the pjates ranging from 
angles of —5° to 55° with respect to the incident deuteron direction, we were able to 
establish the existence of two long range neutron groups. These correspond to the 
ground state and a level of the residual nucleus 2’F at 0:53+ 0:06 mev. The Q-value for this 
1eaction was found to be —1:51+0-05 mev., and from this the mass of !7F was calculated 
to be 17-:007398 + 0:00005 mass units using the mass of the neutron and deuteron as given 
by Bethe (1947). ; 

The angular distributions of the neutron groups were determined by measuring the 
total number of tracks belonging to these two groups in a known volume of emulsion. 
After correcting for the variation with energy of the n—p scattering cross section, the results 
were transferred to the centre-of-mass system and plotted as shown in the Figure. The 
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theoretical curves were calculated by H. Newns on the basis of the theory of Bhatia and 
Huang (unpublished). Also, Butler (private communication) obtained similar curves by 
applying his theory (Butler 1950). In both cases the agreement was achieved only by 
assuming the angular momentum / of the absorbed proton to be 2 in the case of the ground 
state and zero for the first excited state. Since the ground state of 1®O has spin zero and 
even parity, then the ground state of 1*F must have spin 3/2 or 5/2 and even parity, and 
the first excited state must have spin 4 and even parity. 

We conclude by pointing out the similarity between the mirror nuclei ’O and lF. 
The observed level (0:53 Mev.) in 1’F presumably corresponds to the 0°876 mev. level 
in”’O. The experiment of Burrows, Gibson and Rotblat (1950) on the angular distribution 
of the protons from the reaction 1*O(d, p)!”O led them to exactly similar conclusions regarding 
the spins and parities of the ground and first excited state of "O. Recently the spin of the 
ground state of 1’O has been reported by Alder and Yu (1951) using the nuclear induction 
method to be 5/2; it therefore appears likely that the spin of the ground state of 17F may 
also be 5/2. 


George Holt Physics Laboratories, F. A. EL-BepEwt!. 
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Quadrupole Moments of Copper 63 and 65 


Hyperfine structure in the paramagnetic resonance spectrum of the copper Tutton 
salts has been extensively studied in this laboratory since its discovery by the late Dr. R. P. 
Penrose (1949). Ingram (1949) found that in addition to magnetic interaction between 
the electrons and the nucleus, another effect was present which was interpreted as due to 
an electric quadrupole interaction. Abragam and Pryce (1951) have shown how both this 
and the anomalous magnetic interaction can be accounted for theoretically. A brief 
review of the experimental results is presented in this letter. 

The observed spectra have been interpreted in terms of the usual Hamiltonian 
H =B{g, AS, + 8)(A2S_+HySy)}+ ASzel, + BSgl y+ Syly) + QU? 41([+1)}—ypyH_I. 
By the use of crystals where most of the copper ions are replaced by diamagnetic zinc ions, 
the half-width at half intensity of the lines can be reduced to about six gauss. Although 
the spectrum is by no means fully resolved (cf. Ingram 1949), its principal features are in 


accord with the predicted behaviour (Bleaney 1951), and the values of Table 1 are obtained 
for the constants in the Hamiltonian. 


Table 1 
Salt git gs alee ea ies es 
Copper potassium sulphate 2°45 2:14 103 35 inl 
Copper ammonium sulphate 2°47 DAS 1132 20 9 
Copper rubidium sulphate 2°45 2:12 120 9 11 


| A|,|B| and| QJ values are to be multiplied by 10-* cm=}. 


One of the principal drawbacks of the lack of resolution is that the quadrupole moments 
of the two copper isotopes (63, 65: abundances 69%, 31%) could not be determined 
accurately, though they appeared nearly equal. In addition, the relative signs of A, B 
and Q could not be found. 

By growing crystals of diluted copper potassium sulphate from heavy water we have 
succeeded in reducing the half-width of the lines to 2 gauss. This reduction is achieved 
because of the smaller magnetic moment of the deuterons in the water of crystallization, 
which is located immediately around the paramagnetic ion. A typical spectrum at 3 cm. 
wavelength is shown in Figure 1. In Figure 2 some of the lines which are normally 
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Figure 1. Hyperfine structure of the electronic Figure 2. The (4, 3/2) <> (—3, 4) and 
4<> —# transition showing the strong transi- (4, 4) <> (—3, 3/2) transitions at 
tions (4, m) <> (—4, m) in which the nuclear high resolution and intensity. The 
magnetic quantum number m does not change. lines due to the less abundant 
The m=3/2 and m=—3/2 lines of the less isotope 65 are displaced to higher 
abundant isotope 65 are clearly resolved at each fields owing to its larger nuclear 
end of the spectrum. Weak transitions (see magnetic moment. 

Figure 2) in which m changes by 1 are indicated 
by arrows. 


forbidden (if Q=0) are shown, the smaller lines being due to the isotope 65. These are 
shifted outwards owing to its larger nuclear magnetic moment. From a number of such 
spectra we have determined the ratio of the quadrupole terms O,3/Qg¢5. as 1:08+ 0-02, all 
the values lying within this range. This is in agreement with the value 1:081+0-003 
obtained by Becker and Kriiger (1951) from the quadrupole splitting of the nuclear resonance 
spectrum of copper in a diamagnetic salt. The ratio of the magnetic hyperfine splittings 
parallel to the tetragonal axis of the copper ion is 1 : 1:069+ 0-003, in agreement with the 
nuclear resonance ratio 1 : 1:0711+0-0002 (Pound 1948). 
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The relative signs of A, B and Q have been found through the direct effect of the 
external magnetic field on the nuclear magnetic moment, which produces a small asymmetry 
in the spectrum (see Bleaney 1951). To obtain a more pronounced effect, measurements 
were made at 1:2 cm. wavelength (higher field). "These showed that B and Q had the same 
sign for both isotopes, and were opposite in sign to A. This is in agreement with the 
theory of Abragam and Pryce and we may adopt their actual signs for these coefficients. 
The results on the heavy water crystal are summarized in Table 2; the values of the 
coefficients are more accurate than those in Table 1, as the magnetic field was calibrated by 
means of a proton resonance. 


Table 2. Coefficients for CuK, (SO,4). . 6D,O 


Isotope Zu can A B Q 
63 PED AND —96:9 +34-0 +10-341 
65 2:42. —F12 —103+5 +36:3 +9-5+1 


A, B and Q values are to be multiplied by 1074 cm7}. 


There is some variation in the values of g, and B in the plane normal to the tetragonal 
axis, showing that there is present a rhombic component of anisotropy. The values given 
refer to the crystallographic ab-plane, which coincides with the magnetic K,Kz plane in this 
salt. No similar variation in Q could be detected, and the error given is such that all 
measurements lay within this range. 

Application of the theory of Abragam and Pryce shows that the values of the nuclear quad- 
rupole moments are —0-12,+0-01 x 10-*4 cm? for isotope 63 and —0-11,+ 0-01 x 10-4 cm? 
for isotope 65. Previous values were —0:1+0-1 (Schuler and Schmidt 1936), later 
re-calculated as —0:26+0-10, —0-15+0-10 (Brix 1949) for the isotopes 63 abd 65 
respectively. 

The crystals used in the heavy water experiments were grown from two different 
solutions containing copper and zinc in the ratios 1: 1,000 and 1 : 200 respectively. No 
difference in line width could be observed, indicating that the contribution to the half-width 
from other copper ions is less than 1 gauss. At 90° k. a small increase in line width, compared 
with 20°. is attributed to sign—lattice relaxation, and the value hence obtained for the 
relaxation time rt is 5X 10-® second. Allowing for experimental error, this is not signi- 
ficantly different from the values of 9 x 10-8 sec., for copper potassium sulphate (Broer and 
Kemperman 1947) and 6 x 10-8 sec., for copper ammonium sulphate (Gorter 1947). The 
latter values are for the concentrated salts; it appears therefore, that in this case there is 
no great change in y on dilution, unlike the results of Benzie (1951) at helium temperatures. 


Clarendon Laboratory, B. BLEANEY. 
Oxford. K. D. Bowers. 
13th June 1951. D. J. E. INGRAM. 
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REVIEWS OF BOOKS 


Molecular Spectra and Molecular Structure: Vol. 1. Spectra of Diatomic Molecules, 
by GERHARD HERZBERG. Pp. xv+658. 2nd Edition. (New York: Van 
Nostrand; London: Macmillan, 1950.) $9.75; 73s. 6d. ; 


The first edition of Professor Herzberg’s book on the spectra and structure of diatomic 
molecules appeared in 1939 and is well known as an excellent text and reference book on the 
subject. The new edition has been revised, extended and brought up to date, not simply 
by the addition of new paragraphs, but by a most thorough process of alterations and 
improvements affecting almost every page. 

In spite of the considerable increase in the material, the book has fortunately retained its 
character as a textbook. The author has achieved this by making a clear distinction between 
paragraphs dealing with fundamentals and paragraphs dealing with less essential details, and 
by marking this distinction by means of different types of print. An outstanding feature of 
this book is the extreme thoroughness and care with which it has been written, and the 
attention given to every detail. 

New sections on microwave spectra, hyperfine structure in molecular spectra and inten- 
sities of electronic transitions have been added. ‘The list of references has been brought up 
to date, and an appendix covering 80 pages contains a complete table of molecular constants 
for excited as well as for normal states. 

The new edition of this work by one of the first authorities on the subject will be warmly 
welcomed by many physicists and chemists, and will prove invaluable to all those whose 
work is connected with molecular spectra and molecular structure. H. K. 


Negative lons, by H. S. W. Massey. (Cambridge Monographs on Physics.) 
Pp.x+136. 2nd Edition. (Cambridge: University Press, 1950.) 12s. 6d. 


The first edition of this excellent monograph was reviewed in the Proceedings of the 
Physical Society in 1938. This second edition is better bound with hard covers, has a 
bibliography, and is nearly 30° larger, the increase being due mainly to the complete 
re-writing of the last chapter on modern developments in the study of negative ions. 

The first four chapters closely follow those of the first edition apart from small 
modifications due to later work. Chapter I discusses the theory and structure of negative 
atomic ions and describes experimental determinations of electron affinities. (There are 
two small misprints in the last paragraph: 102% being printed for 10-12, and O7 being 
printed for O-~.) Chapter II discusses the nature of molecular ions in general, and of 
the simple cases H,~, O,~, halogen ions and OH~ in particular. Chapter III deals with 
modes of formation of negative ions: it discusses fully the various processes of electron 
attachment. Calculations are given of free-electron radiative capture by atoms of hydrogen 
and of oxygen, and capture by Na, Hg, Cl is also considered : then follows a discussion of 
the possibility of experimental confirmation of these theoretical results. Attachment in 
three-body collisions is next examined, and attachment by molecules (a form of three- 
body problem) is discussed. ‘The two groups of experiments on electron attachment 
are then described : (i) the homogeneous electron beam method of Tate and Lozier and the 
mass spectrograph of Hagstrum and Tate at low pressures; (ii) the electron group methods 
such as the filter method of Loeb and Cravath and the diffusion method of Bailey at higher 
gas pressures. A good review of the main results is given. This important chapter 
closes with a survey of other modes of ion formation by capture of bound electrons from 
neutral atoms or metals. The next chapter discusses processes of electron detachment 
in the light of the principle of detailed balancing. This discussion is mainly theoretical, 
and reflects the lack of experimental data in this field : in fact, one of the features of this 
book is the emphasis it lays on the need for more experimental data on these collisional 
processes, since many of the views now held are, of necessity, conjectural. 
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Detachment is more difficult to study experimentally than attachment, as the main 
desideratum is a high value of the ratio of negative ion to electron concentration, whereas 
for attachment the more readily obtained low values of this ratio are required. 

The last chapter occupies more than a quarter of the book, and deals with applications 
of negative ions and interpretation of collisional processes, whether in glow discharge tubes 
or in the terrestrial or stellar atmospheres. Some possible effects of negative ions in glow 
discharges are outlined qualitatively. Many new data on the earth’s upper atmosphere 
are included, and a careful summary of its properties is given in tabular form. The 
Chapman theory of ionized layer formation is outlined. The influence of negative ions in 
solar and stellar atmospheres is discussed, and it is shown in some detail how negative ions 
of hydrogen practically determine the spectral distribution of solar continuous radiation 
over the observable frequency range, 

This book is very readable, and can be highly recommended for its treatment of the 
subject, ; F. LLEWELLYN JONES. 
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ASTRACTS FOR SECTION B 


Diffraction from a Bent Crystal, by M. BLACKMAN. 


ABSTRACT. The diffraction of plane waves by a bent crystal is considered. It is shown 
that the peak of the diffracted ‘ intensity ’ can be represented in reciprocal space by an arc, 
the limits of the arc being the total angle of bend of the crystal. The ‘ intensity ’ function 
in reciprocal space has been studied numerically in a particular case, and it is shown that 
both ‘intensity ’ and shape change with order. 


Diffraction from a Curved Linear Lattice, by M. BLACKMAN. 


ABSTRACT. The diffraction of plane waves by a curved linear lattice is studied by 
calculating the ‘ intensity ’ in reciprocal space. It is shown that the regions of high intensity 
are arc-shaped in character. Further, the intensity associated with a particular order of 
diffraction (considered as a function of the cylindrical radial coordinate in a plane parallel 
to that containing the curved lattice) can be represented by the square of a Bessel function of 
high index in the neighbourhood of the main peak, but dies off sharply on that side of the 
peak furthest from the origin. ‘This ‘ cut-off’ is characteristic of the total angle of bend 
as long as this angle of bend is less than 180°. The zero order diffraction peak is of a 


different type from the others, but is similar to the peak characteristic of a straight chain ~ 


except for angles of bend greater than 180°. 


The Morphology of Zinc Oxide Smoke Particles, by J. M. Cow.ey, A. L. G. REES 
and J. A. SPINK. 


ABSTRACT. Single-crystal electron diffraction patterns with doubled spots obtained from 
zinc oxide smoke particles originate from pairs of thin sheets which are commonly seen in 
electron micrographs of this material joining up the members of a group of four spines. 
The plane of each sheet is parallel to a (1010)-type lattice plane. The sheets are continuous 
with one of the spines and spaced at 120° intervals around this spine. The other three 
spines, and additional spines which grow from corners of the sheets, are related to this one 
spine and the sheets by twinning on (1122) planes. Intensity measurements of single- 
crystal patterns allow the shape-transform and thickness of the sheets to be deduced. 
Sharp lines and extra spots indicate the occasional presence of lattice imperfections in the 
form of parallel faults. Electron micrographs show wide deviations from the idealized group 
of spines and sheets. 


The Coagulation and Deposition of a Stirred Aerosol, by J. CorNER and E, D. 
PENDLEBURY. 


ABSTRACT. Langstroth and Gillespie have determined experimentally the rate of change 
of paiticle number in ammonium chloride smokes in still and moving air, and have analysed 
their results into deposition on the walls and coagulation inside the aerosol cloud. It is 
shown here that the empirical constants describing these processes can be explained on 
theoretical grounds, both in order of magnitude and as regards their variation with the rate 
of stirring. 


On Oscillations in Electron Streams, by R. Q. Twiss. 


ABSTRACT. A critical discussion is given of the conventional procedure for determining 
whether a particular electron stream is unstable. It is shown that in some cases where this 
analysis suggests instability the stream is stable, while in others, where instability does occur, 
neither the frequency nor the rate of build up of oscillations are as predicted. The physical 
causes of instability in an electron stream are also considered. 
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Dielectric Constant and Electron Density in a Gas Discharge, by R. E. B. MakINson, 
P. C, THONEMANN, R. B. KinG and J. V. Ramsay. 


ABSTRACT. Measurements have been made of the impedance at frequencies of the 
order 1,000 Mc/s. of a coaxial transmission line terminating in a small electrode immersed in 
the plasma of a low pressure discharge in mercury vapour. ‘The frequencies f chosen were 
in the neighbourhood of that for which the real part €, of the dielectric constant of the space 
near the electrode, as given by ¢,=1—ne?/(7mf?) where n is the local electron density, was 
near zero. The effect of making a small change in the size of the electrode was also studied. 
The electron density in the neighbourhood of the electrode was measured by Langmuir 
probes. A minimum in reactance and a maximum in resistance (when plotted against 
discharge current) were found to indicate «,=O near the electrode, in agreement with a 
qualitative argument given. ‘The possible use of such impedance measurements as a means 
of measuring electron density is discussed. 


Conductivity, Hall Effect and Thermo-electric Power of Selenium Single Crystals, 
by K. W. PLESSNER. 


ABSTRACT. Single crystals of selenium have been grown from the vapour phase and 
their conductivity and thermo-electric power were measured over a range of temperatures. 
The Hall effect was determined for one crystal. The conductivity obeys a relatively simple 
exponential law o=o, exp (—eE/RT), with eH=0-13 ev. and the room temperature conduc- 
tivity between 10-5 and 0:5x10-®ohm-!cm~}. Both thermo-electric power and Hall 
effect rise slightly with increasing temperature, indicating a fall in the current carrier con- 
centration. Values calculated from thermo-power and Hall effect agree reasonably well, and 
give about 1014 positive holes per cm®. Because of an experimentally established correlation 
between thermo-power and conductivity, differences in conductivity between various 
crystals are ascribed to different hole concentrations. Deviations from Ohm’s law, the 
rapid increase in mobility with temperature and its low absolute value (about 0-14 cm?/v.sec. 
at room temperature) all point to the existence of potential barriers within the crystals. 
It is shown that a range of barriers of different heights could explain the observed behaviour. 
The conductivity is anisotropic (ratio about 10 : 1), whereas the thermo-power is not. 


Conductivity and Hall Effect of Micro-crystalline Selenium containing Iodine 
Impurities, by K. W. PLESSNER. 


ABSTRACT. ‘The conductivity and Hall effect of micro-crystalline selenium have been 
measured as a function of iodine content under two different conditions of crystallization, 
over a range of temperatures. ‘The conductivity increases with iodine concentration, a 
maximum being reached at a concentration of about 0:03 atomic per cent. The Hall effect 
shows this increase to be due to an increasing mobility, the carrier concentration remaining 
nearly constant. The effect of carrying out the first crystallization of the amorphous Se 
at 175° c. instead of 110° c. is to reduce the conductivity by a factor of ten. This is due to 
reductions in both mobility and current carrier concentration. With increasing temperature 
the mobility rises and the carrier concentration falls, the product of the two resulting in a 
‘rising conductivity with samples crystallized at 175° and a slightly falling conductivity with 
samples crystallized at 110° c. The results are discussed in terms of the grain structure of 
the material. Microscopic examination has shown that the grain size is different for different 
crystallization temperatures and it is suggested that larger crystals contain a smaller density 
of acceptor centres. The fall in carrier concentration with increasing temperature is thought 
to be due to the relief of stress set up due to anisotropic contraction when the material is 
cooled from the annealing temperature. The low absolute value of the mobility and its 
rise with temperature are ascribed to intergranular high resistance barriers. 
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Some Conduction Properties of the Oxides of Cadmium and Nickel, by C. A. 
HoGartTH. 


ABSTRACT. The electrical conductivity and thermo-electric power of specimens of 
cadmium oxide and nickel oxide are measured as functions of temperature and of oxygen 
pressure, and the results are found to agree in general form with theories previously proposed 
by C. Wagner and by the present author. The results are used to verify Wagner’s concept 
of a differential thermo-electric circuit, and a band model for CdO is proposed. 


The Conductivity-Temperature Characteristics of Lead Sulphide: The Influence of 
Oxygen and of the Rate of Heating, by W. EHRENBERG and J. Hirscu. 


ABSTRACT. In general, the conductivity-temperature curves of impurity semiconductors 
should be affected by the rate of heating. An arrangement is described by which such curves 
can be recorded with rates of heating up to 5,000 centigrade degrees per second. Extensive 
tests made on lead sulphide films in an oxygen atmosphere lead to conclusions about the 
stages of oxidation. 


Piezoelectric Behaviour of Partially Plated Square Plates Vibrating in Contour 
Modes, by R. BECHMANN and P. L. Parsons. 


ABSTRACT. In extension of a previous paper, the piezoelectric behaviour of partially 
plated square plates vibrating in contour modes is used to check the theoretical solution for 
the motion. For the longitudinal mode 2 recently described by Ekstein a particularly 
sensitive test of the theory is possible. A critical electrode arrangement should cause the 
piezoelectric excitation to vanish. ‘This is confirmed by experiment. 


The Astigmatic Correction of Telescopes, by R. S. LONGHUuRST. 


ABSTRACT. It is shown theoretically that if a telescope has positive field curvature 
the best possible image on the Gaussian plane is obtained by making the tangential field flat. 
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Magnification =58-5 x 


(a) Maze pattern found after electrolytic (d) (111) plane 
polishing but before annealing. Field 26:5 oersteds. 


(6) Line pattern on (111) plane. (e) Complicated closure pattern on side of 
Field 16-3 oersteds along [110] direction— crystal, i.e. (112) plane. 
Field parallel to [110] direction> 


(c) (111) plane. (f) Pattern on (112) plane showing branching 
Field 25-0 oersteds— at grain boundary. 


Figure 1. 
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